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1.0  ABSTRACT 


This  1-year  SEED  project  advanced  the  field  of  passive  sampling  on  two  major  fronts:  1) 
Developed  an  engineering  approach  to  extend  in-situ  passive  sampling  to  high  molecular  weight 
compounds  in  sediment  porewater,  and  2)  Explored  an  equilibrium  passive  sampling  approach 
for  MeHg  in  sediment  porewater.  Both  proof-of-concept  developments  are  novel  and  were 
successful  in  paving  the  way  for  further  research  to  operationalize  the  extension  of  passive 
sampling  to  these  heretofore  difficult  to  measure  analytes  in  sediment  porewater. 

Passive  sampling  for  the  measurement  of  freely  dissolved  concentrations  of  organic  pollutants  in 
sediment  pore  water  has  emerged  as  a  very  promising  approach,  but  in  situ  measurements  are 
complicated  by  slow  mass  transfer  of  strongly  hydrophobic  compounds.  The  primary  resistance 
to  mass  transfer  arises  in  the  sediment  side  where  a  static  boundary  layer  develops  in  the  vicinity 
of  the  polymeric  passive  sampling  material.  The  slow  mass  transfer  results  in  under-equilibrated 
passive  sampler  measurements  that  need  to  be  corrected  for  equilibrium,  typically  by 
extrapolation  of  the  loss  kinetics  of  performance  reference  compounds.  Such  corrections  are 
prone  to  large  errors,  especially  when  deviation  from  equilibrium  is  large.  In  this  research  we 
address  the  challenge  of  slow  mass  transfer  by  disrupting  the  external  aqueous  boundary  layer 
around  an  in-situ  passive  sampler  using  periodic  mechanical  vibration.  We  report  an  engineering 
innovation  of  adapting  low-cost  motors  used  for  producing  haptic  feedback  in  cell  phones  for  the 
use  in  disrupting  aqueous  boundary  layer  in  a  passive  sampler  deployed  in  sediments.  We 
demonstrate  through  laboratory  experiments  and  numerical  modeling  that  short  periodic  shaking 
of  a  passive  sampler  deployed  in  static  sediment  greatly  enhances  the  rate  of  mass  transfer  and 
reduces  the  difference  in  the  extent  of  equilibrium  achieved  compared  to  a  well  stirred  laboratory 
equilibrium.  The  improvement  over  static  sediment  deployment  is  especially  evident  for  the  high 
molecular  weight  compounds  such  as  benzo(a)pyrene.  We  also  demonstrate  this  method  for 
strongly  hydrophobic  chlorinated  organics  using  PCB  congeners  in  the  log  Kow  range  of  6-8. 
Deployment  of  the  vibrated  passive  sampler  in  laboratory  mesocosms  of  field-collected 
sediments  shows  that  the  porewater  concentrations  of  up  to  octachlorobiphenyls  can  be  measured 
accurately  even  with  a  7-day  deployment  of  the  sampler  in  situ.  In  contrast,  a  static  deployment 
of  passive  samplers  for  28  days  results  in  a  measurement  that  is  4-fold  higher  than  the  true 
equilibrium  value.  Further  work  is  needed  to  operationalize  the  vibrating  passive  sampler 
concept  and  include: 

1)  Further  optimization  of  vibration  frequency  through  laboratory  experiments. 

2)  Testing  on  strongly  hydrophobic  compounds  such  as  dioxins  and  furans  in  sediments. 

3)  Confirmation  through  laboratory  experiments  that  the  concept  of  using  freely  dissolved 
concentrations  in  sediment  porewater  to  predict  biouptake  can  be  extended  to  high  molecular 
weight  PCBs  and  dioxins  and  furans. 

4)  Field  testing  of  the  sampler  platform  along  with  traditional  deployment  of  passive  samplers. 
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This  research  also  laid  the  groundwork  for  the  development  of  a  passive  equilibrium  sampler  for 
MeHg  in  sediment.  Several  polymer  composites  were  developed  and  tested  based  on  the  concept 
of  equilibrium  partitioning  of  MeHg  between  bioavailable  species  in  porewater  and  a  well 
characterized  sorbent  held  within  a  polymer  film.  Among  the  materials  tested,  several 
combinations  showed  the  greatest  promise  based  on  sorption  capacity  and  linearity  of  sorption 
behavior:  activated  carbon  embedded  in  agarose  or  polyvinylidinefluoride  polymer,  cystine  in 
polyethylene  terephthalate,  and  thiols  embedded  in  agarose  or  polyvinylidinefluoride.  While,  the 
presence  of  dissolved  organic  matter  reduced  uptake  in  the  polymers,  this  is  expected  based  on 
speciation  of  MeHg.  Going  forward  we  should  be  able  to  design  the  polymers  to  respond  to  the 
bioavailable  fraction  of  MeHg.  We  believe  that  initial  results  presented  in  this  report  show  great 
promise  for  pursuing  the  concept  of  equilibrium  passive  sampling  for  Hg  and  MeHg. 

Further  work  is  needed  to  advance  the  technology  development  for  passive  sampling  of  Hg  and 
MeHg  and  are  listed  below: 

1)  Complete  testing  of  a  suite  of  selected  polymers  in  a  range  of  MeHg  complexes  prior  to 
testing  against  benthic  organicsms. 

2)  Empirically  determine  which  passive  sampler  sorbents  and  pore  size  predict  MeHg  uptake  by 
benthic  organisms. 

3)  Develop  enriched  stable  isotope  spikes  as  performance  reference  compounds  in  these  samplers 
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1 . 1  List  of  acronyms 


Cfreei  freely  dissolved  concentration 
PRCs:  performance  reference  compounds 
PCBs:  polychlorinated  biphenyls 
PAHs:  poly  aromatic  hydrocarbons 
PE:  polyethylene 

DGT:  diffusive  gradient  in  a  thin  film 

FCVs:  final  chronic  values 

Kpew:  polyethylene-water  partitioning  coefficient 

LP:  half  thickness  of  polyethylene  strip 

Lw:  thickness  of  sediment/porewater 

S:  sediment  concentration 

T:  time 

Kd  :  sediment-water  partition  coefficient 
k:  first  -order  desorption  rate  constant 
Cw:  chemical  concentration  in  water 
Cpe:  chemical  concentration  in  PE 
Dpe:  chemical  diffusivity  in  PE 
p:  sediment  bulk  density 
e:  porosity 

Dw:  diffusivity  in  water 
So:  initial  chemical  concentration  in  sediment 
CwO:  initial  chemical  concentration  in  water 
Ksed-  pe:  sediment-polyethylene  partitioning  coefficient 
ko:  overall  mass  transfer  rate  constant 
5:  diffusion  distance 

PCDD/Fs:  polychorinate  dibenzo-p-dioxins/dibenzofurans 
kf:  fast  desorption  rate  constant 
ks:  slow  desorption  rate  constant 

CpE,eq:  chemical  equilibrium  concentration  in  polyethylene 
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2.0  OBJECTIVES 


Accurate  bioavailability  measurements  are  needed  for  improved  site  risk  assessments,  proper 
selection  of  remedy,  and  post-remediation  monitoring.  While  freely  dissolved  porewater 
concentrations  of  organic  contaminants  and  aqueous  concentrations  of  bioavailable  forms  of 
trace  metals  in  sediments  have  been  demonstrated  to  be  the  ideal  metrics  for  assessing 
bioavailability,  accurate  measurements  of  the  low  environmentally  relevant  concentrations  has 
been  a  challenge.  For  hydrophobic  organics,  the  primary  resistance  to  mass  transfer  is  in  the 
sediment  side  where  a  static  boundary  layer  develops  in  the  vicinity  of  the  polymeric  passive 
sampling  material.  The  slow  mass  transfer  results  in  under-equilibrated  passive  sampler 
measurements  that  need  to  be  corrected  for  equilibrium,  typically  by  extrapolation  of  the  loss 
kinetics  of  performance  reference  compounds.  Such  corrections  are  prone  to  large  errors, 
especially  when  deviation  from  equilibrium  is  large.  In  addition,  there  has  been  little 
development  of  equilibrium  passive  sampling  for  metals.  Key  objectives  of  this  research 
therefore  are  to: 

1)  Explore  a  novel  approach  of  enhancing  mass  transfer  by  introducing  mechanical  vibration  to 
disrupt  the  aqueous  boundary  layer  around  passive  samplers  deployed  in  situ, 

2)  Test  this  approach  for  strongly  hydrophobic  compounds,  and 

3)  Evaluate  a  range  of  polymer  types  for  use  as  passive  equilibrium  samplers  for  methylmercury. 

The  ultimate  goal  is  to  develop  an  in  situ,  actively  shaken  deployment  platform  that  can 
accommodate  multiple  types  of  passive  samplers  (for  strongly  hydrophobic  organics  and  metals). 
This  research  directly  addresses  the  SEED  SON  through  new  innovation  in  passive  sampling 
measurements  for  pollutant  bioavailability  in  sediments.  Specifically,  the  research  addresses  two 
key  needs  described  in  the  SEED  SON: 

i)  Development  of  a  multi-purpose  passive  sampling  device  capable  of  collecting  data  on  several 

contaminants  of  interest.  We  developed  an  actively  shaken  platform  for  the  deployment  of 
equilibrium  passive  samplers  in  the  field.  The  active  shaking  innovation  will  allow  achievement 
of  a  state  of  pollutant  uptake  in  the  passive  samplers  that  is  closer  to  thermodynamic  equilibrium 
and  will  reduce  the  high  uncertainties  associated  with  corrections  using  performance  reference 
compounds  when  the  sampler  is  far  from  equilibrium. 

ii)  Development  of  passive  sampling  strategies  for  very  hydrophobic  chemicals  such  as 
dioxins/furans  and  metals  such  as  Hg.  Our  strategy  for  dioxins/furans  is  to  deal  with  the  issue  of 
slow  mass  transfer  expected  of  these  compounds  by  actively  enhancing  in  situ  mass  transfer  rates 
using  mechanical  disruption  of  the  aqueous  boundary  layer  around  the  passive  sampler.  This  has 
never  been  tried  before  and  is  a  promising  innovation  in  the  field  of  passive  sampling  that 
addresses  the  high  uncertainty  associated  with  performance  reference  compounds  when  the 
extent  of  equilibrium  achieved  is  low.  For  mercury,  we  focused  on  the  species  most  relevant 
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from  a  risk  standpoint,  methylmercury  (MeHg),  and  to  identify  a  polymer  partitioning  approach 
by  developing  and  testing  a  range  of  polymeric  material  types.  We  also  anticipate  that  future 
work  can  expand  this  approach  to  include  other  metals  of  interest  such  as  As,  Zn,  Pb,  Cd,  and 
inorganic  Hg  through  the  right  choice  of  ligands  embedded  in  the  polymer. 


3.0  BACKGROUND 


3.1  Relevance  of  porewater  measurements  for  predicting  bioaccumulation  of 
dioxins/furans  and  MeHg. 

Several  studies  have  shown  that  freely  dissolved  concentrations  (Cfree)  of  organics  in  sediment 
porewater  can  be  related  to  toxicity  and  bioaccumulation  (Mayer  et  al.  2014;  Peijnenburg  et  al. 
2014).  Our  recent  work  with  soil  contaminated  with  dioxins  and  furans  demonstrated  that  uptake 
by  earthworms  is  correlated  to  freely  dissolved  concentrations  measured  using  passive 
equilibrium  sampling  (Figure  3.1a;  Fagervold  et  al.  2010  Supporting  Information).  Similarly,  we 
showed  that  uptake  of  MeHg  by  sediment-dwelling  worms  is  correlated  to  filterable  MeHg 
concentrations  in  sediment  porewater  (Figure  3.1b;  Gilmour  et  al.  2013).  In  both  of  these  cases, 
the  measurements  of  porewater  concentrations  were  made  in  ex  situ  laboratory  studies.  There  is  a 
great  need  to  be  able  to  perform  concentration  measurements  in  sediment  porewater  and  surface 
waters  in  the  field  to  capture  actual  exposure  concentrations  for  native  biota. 


a.  dioxins/furans 


b.  methylmercury 
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Figure  3.1.  Correlation  of  dioxin/furan  uptake  in  earthworms  with  uptake  in  passive  samplers  (a; 
Fagervold  et  al.  2010)  and  correlation  of  MeHg  uptake  in  aquatic  worms  to  porewater  MeHg 
measured  directly  (b;  Gilmour  et  al.  2013). 
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3.2  Challenges  for  in  situ  passive  sampling. 

The  in  situ  approach  is  used  when  it  is  critical  to  capture  conditions  in  the  field.  In  this  approach, 
the  polymer  is  inserted  directly  into  sediments  or  suspended  in  the  water  column  above  the 
sediment  in  the  field  and  left  in  place  for  sufficient  duration  to  allow  the  derivation  of  Cfree  (Oen 
et  al.,  2011;  Fernandez  et  al.  2009a, b;  Beckingham  and  Ghosh  2013;  Lampert  et  al.  2013). 
However,  the  ability  to  attain  equilibrium  and  demonstrate  that  equilibrium  has  been  achieved  is 
often  more  difficult  for  the  in  situ  approach  as  compared  to  the  ex  situ  approach.  To  overcome 
these  challenges,  Huckins  et  al.  (1993,  2002)  suggested  the  use  of  performance  reference 
compounds  (PRCs)  to  calculate  Cfree  from  non-equilibrium  PSM  measurements  (Cp).  PRCs  are 
analytically  noninterfering  chemicals  that  are  embedded  in  the  passive  sampler  prior  to 
environmental  exposure  (Huckins  et  al.  2002).  Examples  of  surrogate  chemicals  include  stable 
isotope-labeled  or  deuterated  forms  of  the  analytes  of  interest,  substances  with  a  log  Kow  that  is 
similar  to  that  of  the  target  analytes  (Huckins  et  al.  2002;  Fernandez  et  al.  2009a),  or  rare  PCB 
congeners  (Tomaszewski  and  Luthy  2008).  The  depletion  rate  of  a  PRC  during  sampler 
deployment  reflects  the  uptake  rates  of  a  target  analyte,  assuming  isotropic  exchange  kinetics 
occur  (Figure  3.2).  Because  of  the  differences  in  the  compound  properties  for  the  PRC  and  the 
target  analyte,  correction  is  needed  to  calculate  the  fractional  approach  to  equilibrium  for  the 
target  analyte  (C(o/C(ss))  from  the  fractional  PRC  dissipation  (I-CprciVCprcco))  at  time  t.  In 
addition,  PRC  correction  becomes  difficult  if  sorption  in  the  surrounding  media  is  concentration- 
dependent  or  compound-dependent.  Several  approaches  to  calibration  using  PRC  data  have  been 
suggested  (Fernandez  et  al.  2009b;  Huckins  et  al.  2006;  Tomaszewski  and  Futhy  2008;  Reible 
and  Fotufo  2012).  In  all  of  these  approaches,  the  uncertainties  introduced  by  the  PRC  correction 
are  larger  when  the  extent  of  equilibrium  is  low,  which  is  the  case  for  strongly  hydrophobic 
compounds  such  as  dioxins  and  furans. 


correction  for 

compound  properties  PRC  uptake  (mirrored) 


Figure  3.2.  PRC  dissipation 
and  compound  uptake  kinetics 
generally  assumed  for  the 
performance  reference 
compound  (PRC)  approach. 
C(t)  and  C(ss)  refer  to  target 
analyte  concentrations  in  the 
passive  sampler  at  time  t  and 
steady  state,  respectively; 
CPRC(t)  and  CPRC(0)  refer  to 
PRC  concentrations  in  the 
passive  sampler  at  time  t  and  0, 
respectively.  (Ghosh  et  al. 
2014) 
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3.3  Measurement  of  pore  water  dioxins/furans. 

The  field  of  passive  sampling  has  been  sufficiently  advanced  for  organic  compounds  in  the  low- 
to  mid-range  of  hydrophobicity,  including  PCBs  up  to  pentachlorobiphenyls  and  PAHs  at  least 
up  to  chrysene.  The  USEPA  has  published  a  guidance  document  for  the  monitoring  of  organic 
contaminants  at  Superfund  sediment  sites  which  discusses  the  use  of  passive  sampling  (USEPA 
2012).  A  recently  published  practical  guidance  (Ghosh  et  al.  2014)  developed  as  part  of  a 
SET AC  Pellston  workshop  on  passive  sampling  provides  the  current  consensus  for  the  use  of 
passive  sampling  for  PAHs  and  PCBs.  By  contrast,  very  little  information  exists  on  the  use  of 
passive  sampling  for  dioxins  and  furans.  Our  previous  work  with  dioxin-contaminated  soils 
demonstrated  that  freely  dissolved  equilibrium  concentrations  in  soil  correlated  well  with  uptake 
in  earthworms  (Figure  3.1a;  Fagervold  et  al.  2010).  In  these  studies,  we  used  a  well-mixed  batch 
system  to  achieve  close  to  an  equilibrium  state.  We  anticipate  that  the  major  obstacle  to 
performing  in  situ  passive  sampling  measurements  of  dioxins  and  furans  in  porewater  will  be  the 
slow  external  mass  transfer  for  these  strongly  hydrophobic  compounds  (log  Kow  range  from  6.4- 
8.8). 

3.4  Measurement  of  porewater  methylmercury. 

Accurate  information  regarding  the  amounts  and  distribution  of  the  various  species  of  mercury  in 
sediment  and  other  aquatic  environmental  matrices  is  critical  to  a  realistic  assessment  of  risk. 
Methylmercury,  the  mercury  species  of  greatest  toxicological  concern,  is  formed  primarily  by 
anaerobic  microorganisms,  which  take  up  and  methylate  inorganic  mercury  species.  In  sediment 
environments,  MeHg  is  taken  up  by  benthic  and  epibenthic  organisms  and  undergoes  substantial 
biomagnification,  ultimately  leading  to  toxic  risk  for  humans  (Mason  and  Lawrence  1999). 
Porewater  concentration  measurements  are  most  useful  to  risk  assessors  when  they  can  be  used 
to  predict  MeHg  uptake  and  bioaccumulation  potential  in  these  lower  trophic  level  organisms 
residing  at  the  base  of  the  aquatic  food  chain,  but  to  date  no  passive  sampling  scheme  has 
achieved  wide  acceptance  for  this  purpose. 

Neither  the  free  MeHg  cation  nor  simple  MeHg  salts  (e.g.  MeHgCl)  are  dominant  forms  of 
MeHg  in  sediment  pore  waters.  Rather,  MeHg  is  mainly  complexed  with  colloidal  natural  organic 
matter  (Bergamaschi  et  al.  2012);  which  can  enhance  partitioning  from  sediments  into  the  total 
porewater  phase  (Hammerschmidt  et  al.  2004;  Hollweg  et  al.  2009)  but  can  limit  uptake  by 
benthic  infauna  (Luengen  et  al.  2012;  Lawrence  and  Mason  2001).  Methylmercury  thiols  may  be 
particularly  bioavailable  (Leaner  and  Mason  2001),  but  overall  bioavailability  of  MeHg  in 
sediment  porewaters  remains  poorly  understood.  The  bioavailability  of  MeHg  to  organisms  is 
highly  dependent  on  phase  and  complexation.  Like  PCBs,  the  sediment:water  partition 
coefficient  for  MeHg  can  be  a  strong  predictor  of  bioavailability  to  benthic  fauna  (Wang  et  al. 
1998;  Williams  et  al.  2010).  Further,  we  found  that  the  concentration  of  filter-passing  MeHg  is 
well  correlated  with  uptake  by  benthic  organisms  in  laboratory  tests  (Gilmour  et  al.  2013). 
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Currently,  the  most  widely  used  passive  sampler  for  MeHg  is  the  diffusive  gradient  in  a  thin  film 
(DGT)  device.  Briefly,  analytes  diffuse  through  a  gel  of  known  thickness  and  the  molecules 
which  traverse  the  entire  gel  within  the  sampling  time  are  tightly  bound  to  a  high-affinity 
stationary  phase  which  is  collected  and  extracted  for  mercury  analysis  (Divis  et  al.  2005a).  DGT 
has  been  developed  and  refined  for  roughly  twenty  years,  but  doubts  remain  as  to  its  ability  to 
generate  reliable  estimates  of  bioavailability,  particularly  in  sediment.  For  one,  the  question  of 
which  species  are  measured  (that  is,  which  ones  are  sufficiently  “labile”  to  penetrate  the  gel  and 
reach  the  stationary  phase  during  the  deployment  time)  is  not  well  characterized  (FIsu-Kim  et  al. 
2013).  In  addition,  the  fact  that  DGT  operates  in  a  kinetic  (as  opposed  to  equilibrium)  uptake 
regime  presents  an  important  challenge  unique  to  sediment  environments.  Concerns  persist 
regarding  porewater  depletion  and  resupply  of  MeHg  in  the  vicinity  of  the  DGT  device  (Clarisse 
et  al.  201 1).  Passive  polymer  sampler  surfaces  may  offer  an  empirical  approach  to  monitoring 
bioavailability  by  mimicking  the  partitioning  of  MeHg  complexes  in  porewaters  to  cell  surfaces. 
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4.0  ORGANISATION  OF  THE  REPORT 


This  research  project  entailed  three  key  research  objectives  as  described  below.  Results  and 
discussion  associated  with  each  of  the  objectives  have  been  written  separately,  with  the  ultimate 
goal  of  being  able  to  submit  each  for  separate  peer-reviewed  publications.  Also,  keeping  them 
separate  allows  more  fussed  discussion  of  each  objective  separately.  The  first  document 
associated  with  objective  1  has  been  published  in  the  journal  Environmental  Science  and 
Technology  (Jalalizadeh  and  Ghosh,  Environ.  Sci.  Technol.  2016,  50,  8741-8749).  The  other 
two  are  still  being  worked  into  final  manuscripts. 

1)  Explore  a  novel  approach  of  enhancing  mass  transfer  by  introducing  mechanical  vibration  to 
disrupt  the  aqueous  boundary  layer  around  passive  samplers  deployed  in  situ.  (Chapter  5) 

2)  Test  this  approach  for  strongly  hydrophobic  compounds  like  dioxins  and  furans.  (Chapter  6) 

3)  Evaluate  a  range  of  polymer  types  for  use  as  passive  equilibrium  samplers  for  methylmercury. 
(Chapter  7) 
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5.0  RESEARCH  OBJECTIVE  #1:  Explore  a  novel  approach  of 
enhancing  mass  transfer  by  introducing  mechanical  vibration  to  disrupt 
the  aqueous  boundary  layer  around  passive  samplers  deployed  in  situ. 


5.1  Introduction 

The  freely  dissolved  concentration  of  hydrophobic  pollutants  in  sediment  porewater  is  a  critical 
measurement  that  is  useful  in  assessing  fate,  transport,  and  bioavailability  of  these  compounds.7 
Accurate  measurement  of  low  aqueous  concentrations  of  hydrophobic  compounds  is  challenging 
due  to  the  association  with  colloidal  and  dissolved  organic  matter  in  porewater.  This  has  led  to 
the  development  of  passive  sampling  approaches  using  well  characterized  polymeric  materials. 
When  the  polymeric  material  is  able  to  reach  equilibrium  with  the  sediment  porewater,  such  as  in 
a  well-stirred  laboratory  measurement,  the  estimation  of  freely  dissolved  porewater 
concentration  (Cfree)  becomes  a  trivial  exercise  based  on  the  measured  partition  constant  of  the 
polymer.2-3  However,  in  several  situations,  an  in  situ  measurement  in  sediment  is  desired.  Such 
measurements  have  been  challenged  by  the  difficulty  in  reaching  equilibrium  between  porewater 
and  the  passive  sampler  as  mass  transfer  through  the  depletion  layer  outside  the  commonly  used 
polymers  (with  commonly  used  thicknesses)  becomes  limiting  in  the  absence  of  active  mixing.4 
It  has  been  shown  that  for  strongly  hydrophobic  compounds  equilibrium  may  not  be  achieved  in 
the  field  even  after  one  year  in  51  pm  polyethylene  (PE).5  Several  researchers  have  adopted  the 
use  of  performance  reference  compounds  (PRCs)  dosed  in  the  polymer  to  assess  the  kinetics  of 
mass  transfer  and  correct  for  non-equilibrium.6'74  While  corrections  based  on  PRC  loss  work 
reasonably  well  for  compounds  with  low  to  mid-range  hydrophobicity,  the  corrections  become 
increasingly  erroneous  for  strongly  hydrophobic  compounds  when  the  departure  from 
equilibrium  increases.74  Several  approaches  for  calibration  using  PRC  data  have  been 
suggested.6'70  In  all  of  these  approaches,  the  uncertainties  introduced  by  the  PRC  correction  are 
larger  when  the  extent  of  equilibrium  is  low,  which  is  the  case  for  strongly  hydrophobic 
compounds  in  the  field. 

A  primary  uncertainty  in  the  PRC  correction  arises  from  the  fact  that  nearly  always  it  is  the 
sediment  side  mass  transfer  in  the  immediate  vicinity  of  the  passive  sampler  that  controls  overall 
kinetics  4  The  retarded  diffusion  in  the  sediment  side  is  impacted  by  the  site-specific  sorption 
characteristics  of  the  sediment,  which  can  vary  across  orders  of  magnitude.  For  example 
Hawthorne  et  al.75  reported  a  3-4  orders  of  magnitude  range  for  site-specific  Koc  values  for 
polycyclic  aromatic  hydrocarbons  (PAHs).  Thus,  to  be  able  to  correct  for  non-equilibrium  and 
estimate  in  situ  porewater  concentrations  we  need  to  first  have  an  estimate  of  site-specific 
partitioning  of  the  analytes  of  interest.  The  loss  kinetics  of  a  few  PRC  compounds  into  sediment 
is  used  to  infer  the  desorption  behavior  of  a  large  range  of  analyte  compounds  from  sediment.8 

Thinner  polymeric  materials  can  be  used  to  increase  the  surface  area  to  volume  ratio  and  reduce 
the  depletion  per  unit  area6  and  reduce  deployment  time.  However,  even  with  some  of  the 
thinnest  polymers  practically  deployable  in  the  field  (e.g.  25  pm  thick  PE),  sediment-side  mass 
transfer  limitation  can  be  significant.4  Making  the  polymers  too  thin  makes  them  prone  to 
damage  during  deployment  in  sediment,  reduces  the  total  mass  of  polymer  sampling  material 


16 


(impacting  detection  limits),  and  also  poses  a  physical  challenge  of  insertion  in  sediments  if  the 
polymer  surface  area  is  very  large. 

To  address  these  challenges,  we  take  a  very  different  approach  of  manipulating  the  external 
depletion  layer  in  the  sediment  side  of  a  passive  sampler.  We  demonstrate  for  the  first  time  a 
novel  approach  of  in  situ  passive  sampling  that  overcomes  the  slow  approach  to  equilibrium  for 
hydrophobic  organic  compounds  in  static  sediments.  Our  approach  involves  mechanical 
disruption  of  the  depletion  layer  outside  the  polymer  surface  using  periodic  vibration  performed 
in  situ.  We  adapted  small  vibration  motors  used  for  haptic  feedback  in  cell  phones  and  other 
electronic  devices  to  create  vibrations  in  in  situ  passive  sampling  frames  programmed  to  trigger 
at  pre-determined  time  intervals.  The  concept  of  passive  sampling  still  applies  to  the  vibrated 
deployment  because  the  disruptions  introduced  by  the  periodic  vibration  is  not  unlike  natural 
sediment  processes  induced  by  burying  and  sediment  ingesting  worms  that  also  achieve  chemical 
equilibrium  through  small-scale  physical  disruptions  in  sediments.  Laboratory  experiments  were 
performed  using  PAH-impacted  field  sediments  to  compare  the  approach  to  equilibrium  under 
static,  well-mixed,  and  differently  shaken  passive  sampling  modes.  We  also  performed  numerical 
modeling  of  the  mass  transfer  process  to  mechanistically  explain  our  observations  and  optimize 
the  duration  and  periodicity  of  the  vibrations  in  order  to  minimize  the  energy  requirement  to 
drive  the  vibration  motors. 


5.2  Materials  and  methods 

Materials.  Low  density  PE  sheets  (25pm  thickness),  manufactured  by  Poly- America  (Grand 
Prairie,  TX,  USA)  were  purchased  from  the  Home  Depot.  PAH  and  deuterated  PAH  stock 
solutions  were  purchased  from  Fisher  Scientific  (Pittsburg,  PA,  USA).  Cylindrical  vibrating 
motors  with  a  diameter  of  9  mm  and  length  of  25  mm  were  purchased  from  Precision 
Microdrives  (London,  UK).  The  motors  operated  at  a  rated  voltage  of  3V,  operating  current  of 
130mA,  and  vibrating  speed  of  13,500  rpm.  Pulse-pause  timers  (model  60H)  were  purchased 
from  Velleman  Inc.  (Fort  Worth,  TX,  USA).  Prior  to  use,  PE  sheets  were  soaked  twice  in 
hexane/  acetone  (50/50)  and  left  on  a  shaker  for  24  h  each  time  to  remove  oligomers  and  any 
target  and  non-target  contaminants.  Clean  PE  sheets  were  then  cut  into  6mm  x  2cm  strips  (2.8 
mg)  and  were  soaked  in  a  PRC  solution  (80:20  methanol: water  with  pyrene-dlO  and 
phenanthrene-dlO)  and  allowed  to  equilibrate  for  15  days  on  an  orbital  shaker. 

After  impregnation,  all  strips  were  removed  from  the  PRC  solution  and  rinsed  with  DI  water  to 
remove  residual  methanol  on  the  surface.  Two  strips  were  extracted  immediately  in  a  1:1  hexane 
and  acetone  mixture  (3  x  24  h)  to  determine  the  initial  PRC  concentration  in  the  PE  strips. 


In  situ  shaken  passive  sampler  design.  Two  motors  were  connected  in  parallel  to  a  timer  and  a 
power  supply  (2  rechargeable  batteries  1.2  V,  700  mAh  each)  (Figure  5.1).  The  timers  were 
powered  by  a  12  V  power  supply  and  were  programmed  to  control  motor  vibration  duration  and 
frequency:  5  seconds  pulse  and  2  minutes  pause  (high  frequency),  2  seconds  pulse  and  5  minutes 
pause  (low  frequency).  The  small  PE  sheets  described  above  were  enclosed  in  stainless  steel 
mesh  (14  mesh  from  TWP  Inc,  Berkeley,  CA)  and  were  attached  like  radial  fins  on  cylindrical 
motors  (four  fins  on  each  motor)  as  shown  in  Figure  5.1  and  SI  Figures  S5.1-5.3. 
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Figure  5.1-  Wire  diagram  of  a  vibrating  passive  sampler.  PE  samplers,  enclosed  in  stainless  steel  mesh, 
were  attached  to  the  motors.  Two  motors  were  connected  in  parallel  to  a  timer  and  a  power  supply  (2 
rechargeable  batteries  1.2  V  each).  The  timers  were  powered  by  a  12  V  battery. 


PAH  contaminated  river  sediment  samples  from  the  vicinity  of  former  manufactured  gas  plants 
studied  previously  (identified  as  HD-3  and  HD-5  in  Khalil  et  al/6)  were  used  in  the  present 
work.  The  two  sediments  were  combined  to  obtain  a  sufficient  quantity  for  the  present 
experiments.  The  sediment  was  mixed  with  DI  water  containing  200  mg/L  sodium  azide  to 
make  a  slurry  with  weight  ratio  of  1:2  (dry  sediment/water).  Water  addition  was  performed  to 
make  the  sediment  more  fluid  and  allow  good  mixing  in  the  jars/  For  the  well-mixed  exposure, 
250  mL  of  the  sediment-water  slurry  and  eight  of  the  prepared  PE  samplers  were  transferred  to  a 
wide  mouth  jar  with  a  Teflon-lined  cap.  The  jar  was  placed  on  a  rotary  agitator  and  tumbled  at 
28  rpm.  The  rest  of  the  sediment  slurry  was  placed  in  a  large  glass  tray  (25cmx35cmx6cm).  Four 
motors  were  placed  inside  the  sediment  in  the  tray  with  10  cm  separation  to  prevent  influencing 
each  other.  Two  motors  vibrated  at  the  high  frequency  mode  and  two  at  the  low  frequency  mode. 
Eight  additional  PE  samplers,  enclosed  in  stainless  steel  mesh  without  motors,  were  placed 
inside  the  same  tray  in  a  static  mode  at  least  10  cm  away  from  any  vibrating  motor  to  simulate  a 
static  deployment.  The  tray  was  then  covered  with  aluminum  foil.  Two  PE  strips  were  removed 
from  the  well  mixed,  static,  and  vibrating  systems  and  analyzed  for  native  PAHs  and  PRC 
compounds  after  7,  14,  28  and  56  days. 


18 


PAH  extraction  and  analysis.  Upon  removal  from  the  sediment,  PE  strips  were  rinsed  with  water 
and  wiped  with  laboratory  tissue  to  remove  water  and  adhering  particles.  Prior  to  extraction 
anthracene-dlO  surrogate  was  added  to  assess  the  effectiveness  of  sample  processing,  and 
extracts  with  lower  than  80%  surrogate  recoveries  were  discarded.  Samplers  were  extracted  with 
a  1:1  hexane  and  acetone  mixture  (3  x  24  h,  with  sequential  extracts  pooled).  The  final  extraction 
volumes  were  blown  down  to  lmL  using  a  stream  of  nitrogen  gas.  PAHs  from  sediment  were 
extracted  by  sonication  (EPA  method  3550B)  and  cleaned  using  activated  silica  gel  (EPA 
method  3630C).  Four  internal  standards  were  added  to  the  final  extracts  before  analysis  (1- 
fluoronapthalene,  p-terphenyl-dl4,  benzo(a)pyrene-dl2,  and  dibenz(a,h)anthracene-dl4).  The  16 
USEPA  priority  pollutant  PAHs  were  analyzed  in  an  Agilent  6890  gas  chromatograph  coupled  to 
an  Agilent  5973N  MS  detector  as  described  in  Khalil  et  alA6 


Modeling  uptake  of  analytes  from  sediment  porewater 

A  numerical  modeling  approach  was  used  to  describe  the  mass  transfer  process  involved  in  the 
static,  vibrating,  and  well-mixed  deployments  of  passive  samplers.  The  modeling  structure  for 
static  and  vibrating  deployments  are  similar  and  are  described  first  with  respective  boundary 
conditions.  Finally  the  well-mixed  model  is  described  where  mass  transfer  limitation  is  only  in 
the  polymer  side. 

Static  system.  Two  modeling  approaches  were  implemented  to  simulate  mass  transfer  of  PAHs 
from  sediment  particles  into  porewater,  and  from  porewater  into  polymer  during  static 
deployment.  The  first  approach  is  based  on  the  one-dimensional  diffusion  model  presented  by 
Fernandez  et  al.,s  but  solved  numerically.  In  this  model  mass  transfer  within  the  sediment 
particle  is  assumed  to  be  fast  and  instantaneous  equilibrium  between  sediment  particles  and 
porewater  is  assumed  (local  equilibrium  model).  Mass  transfer  within  the  polymer  and  in 
porewater  is  described  by  Fick’s  2nd  Law  of  diffusion.  Instantaneous  equilibrium  is  assumed  at 
the  polymer  surface  boundary  with  porewater  as  done  in  previous  work.4  5  This  model  is 
explained  in  more  detail  in  the  SI. 

The  second  approach  is  based  on  the  assumption  that  mass  transfer  from  sediment  particles  into 
porewater  is  described  by  first  order  kinetics  A /A  The  model  was  solved  using  two  different 
hypothesis:  1)  All  PAHs  are  associated  with  the  slow  desorbing  pool  in  sediment  and  desorption 
is  characterized  by  the  slow  desorption  rate  constant.  2)  All  PAHs  are  associated  with  the  fast 
desorbing  pool  in  sediment  and  desorption  is  characterized  by  the  fast  desorption  rate  constant. 
First  order  desorption  rate  constants  were  obtained  from  a  study  by  Ghosh  et  al./7,  where 
sediment  from  a  similar  manufactured  gas  plant  impacted  site  was  used  (see  SI). 

In  a  system  containing  a  PE  strip  with  the  thickness  of  21P  and  sediment/porewater  with  thickness 
of  lw  on  both  sides,  sediment  concentration  (S)  changes  following  first  order  kinetics  A /,s‘ 

^  =  k(KdCw  -  S)  lp  <  x  <  lp+lw  and  -  lp  <  x  <  -lp-lw  (1) 
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Where  t  is  time  (s),  S  is  chemical  concentration  in  sediment  (ng/g),  Kd  is  sediment-water 
partition  coefficient  (cm3/g),  k  is  first  order  desorption  rate  constant  (s'1)  and  Cw  is  chemical 
concentration  in  water  (ng/cm3) 


For  a  PE  strip  with  concentration  Cpe  and  at  point  x  and  time  t: 


dCPE  _  n  d2CpE 

~dT  ~  UpE  ~d^~ 


—lp  <  x  <  lp 


(2) 


Where  Cpe  is  chemical  concentration  in  PE  (ng/cm3)  and  Dpe  is  chemical  diffusivity  in  PE 
(cm2/s) 


The  transport  equation  in  porewater  with  concentration  of  Cw  at  point  x  and  time  t  will  be  as 
follows: 

^f  =  D^r-(0^  lp  <  x  <  lp+lw  and  —  lp  <  x  <  —lp—lw  (3) 

Where  p  is  sediment  bulk  density  (g/cm3)  and  s  is  porosity  (cm3/  cm3) 

D  is  the  diffusivity  in  water  (Dw)  after  correction  for  tortuosity: 


Substituting  Equation  (1)  in  Equation  (3),  the  transport  equation  in  porewater  can  be  re-written 
as: 


d2Cw 
dt  dx2 


0) k(KdCw-S ) 


(5) 


Initial  conditions 

The  polymer  was  initially  clean  and  porewater  was  assumed  to  be  in  equilibrium  with  sediment. 
CPE  —  0  —Ip  <  x  <  Ip  (6) 

Cwo=so/Kd  lp  <  x  <  lp+lw  and  -  lp  <  x  < -lp-lw  (7) 

Where  So  is  initial  chemical  concentration  in  sediment  (ng/g)  and  Cwo  is  initial  chemical 
concentration  in  water  (ng/cm3) 

Boundary  conditions 

Continuity  of  flux  and  equilibrium  condition  was  assumed  at  the  PE-water  boundary  as  done  in 
previous  work.8 

n  dCpE  =  Ddc^ 
rCj  dxPE  dxw 

CpE  —  KpEw^w 


x  =  lp  and  x  =  —lp,  t  >  0  (8) 

x  —  lp  and  x  —  —lp  (9) 
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Where  Kpew  is  the  PE-water  partition  coefficient  (cm3/cm3) 

Due  to  symmetry  the  flux  will  be  zero  at  the  center  of  the  PE  sheet.  Porewater  concentration  is 
equal  to  the  initial  concentration  far  away  from  the  polymer  at  x=lP+lw  and  does  not  change  over 


time: 

dcPE  _  q 
dx 

o 

II 

X 

(10) 

dCw  =dS  = 
dt  dt 

X  =  l  +1 
*A'  Lp  1  LW 

and  x  —  —lp—l 

Vibrating  system.  In  the  vibrating  system,  when  the  motor  is  in  pause  mode,  the  mass  transfer  is 
similar  to  the  static  mode  and  sediment  and  porewater  concentration  in  the  vicinity  of  the 
polymer  depletes  with  time.  Every  time  the  motor  vibrates,  the  sediment  and  porewater  in  the 
vicinity  of  the  polymer  is  mixed  up.  Our  visual  observation  of  the  vibration  system  indicated  that 
the  extent  of  fluidized  mixing  was  at  least  1  cm  from  the  surface  of  the  device  with  pressure 
waves  from  the  vibration  extending  to  2-3  cm  from  the  device  (SI  Figure  S5.3).  We  assumed  that 
this  mixing  is  enough  to  increase  sediment  concentration  to  the  initial  concentration  in  sediment 
(So).  Porewater  concentration  right  after  each  vibration  pulse  will  also  increase  to  the  initial 
concentration: 


cw  Cwq  ,  S  Sq  Ip  x  Zp-fZ^y  and  Ip  x  Ip  Z^,  (13) 

Well-mixed  system.  Diffusion  in  a  polymer  with  thickness  of  21P  follows  Fick’s  second  law: 

~^r  ~  dpe  dx 2  ~lp  <  x  <  Ip  (13) 

The  passive  sampler  was  assumed  to  be  initially  clean.  Since  the  system  is  perfectly  mixed,  the 
porewater  concentration  remains  constant  and  equal  to  the  initial  value  (So/Kd)  during  the 
deployment  time.  The  boundary  condition  is  defined  as: 

CPE  —  KPEwCw0  x  —  Ip  and  x  —  —Ip  (14) 

The  model  equations  were  solved  in  Matlab  using  an  explicit,  finite-difference  numerical 
modeling  technique.39  To  check  the  numerical  model,  the  fractional  uptake  of  chrysene  into  PE 
was  simulated  and  the  solution  was  found  to  be  identical  to  the  analytical  solution  based  on 
Fernandez  et  al.s  Details  of  parameter  estimation,  numerical  solutions  and  Matlab  codes  are 
provided  in  the  SI. 


5.3  Results  and  discussion 

Equilibrium  PAH  concentration  in  sediment  and  porewater.  The  concentrations  of  the  16  EPA 
priority  pollutant  PAHs  in  sediment  are  shown  in  Figure  5.2.  The  four  most  abundant  PAHs  were 
phenanthrene,  fluoranthene,  pyrene,  and  benzo(a)pyrene.  The  total  concentration  of  PAHs  in  the 
sediment  was  128  pg/g.  The  2-4  ring  PAHs  comprised  60%  of  the  total  in  sediment.  The 
equilibrium  concentrations  in  PE  (CpE,eq)  were  determined  by  extracting  the  PE  samplers  after  56 
days  of  deployment  in  the  well-mixed  PE-sediment  system.  Another  time  point  measurement  of 
PE  concentration  after  77  days  of  deployment  confirmed  that  equilibrium  had  reached  in  PE  for 
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all  PAHs  in  56  days  (SI  Table  S5.3).  Both  PRCs  (pyrene-dlO  and  phenanthrene-dlO)  were 
depleted  completely  in  the  well-mixed  exposure  in  30  days  (SI  Figure  S5.5).  The  freely 
dissolved  porewater  concentrations  were  calculated  from  equilibrium  concentration  of  PAHs  in 
PE  and  reported  values  for  Kpew5  (Equation  9). 

As  shown  in  Figure  5.2,  acenaphthene  and  phenanthrene  are  the  dominant  PAHs  in  porewater. 
PAHs  up  to  chrysene  contributed  to  99%  of  the  total  porewater  concentration  at  equilibrium. 
Toxic  units  were  estimated  for  each  PAH  by  dividing  porewater  concentrations  by  final  chronic 
values  (FCVs).20  The  total  toxic  units  of  16  PAHs  measured  in  the  sediment  porewater  was 
approximately  9.33  indicating  that  PAHs  in  this  sediment  likely  pose  narcosis  toxicity  to  benthic 
invertebrates.20 


Figure  5.2-  PAH  concentration  in  porewater  (pg/L)  and  sediment  (pg/g).  Sediment  concentrations  are 
mean  values  from  triplicate  measurements.  Porewater  concentrations  were  obtained  using  eq  9  and  CpE,eq, 
which  are  mean  values  from  duplicate  measurements  of  PE  concentration  after  56  days  of  exposure  in  the 
well-mixed  system. 


PAH  uptake  in  static  deployment.  As  expected,  the  uptake  of  all  PAHs  was  slowest  in  the  static 
exposure.  In  fact,  none  of  the  PAHs  from  fluoranthene  to  indeno(l,2,3,-cd)pyrene  reached 
equilibrium  in  56  days  of  contact.  The  fractional  uptake  of  benz(a)anthracene, 
benzo(k)fluoranthene  and  indeno(l,2,3,-cd)pyrene  were  only  35%,  11%,  and  8%,  respectively 
after  28  days.  Previous  studies  have  reported  similar  slow  uptake  in  a  static  exposure/7'25 
especially  for  larger  molecular  weight  compounds.  For  example,  Fernandez  et  al.77  calculated  the 
fractional  equilibration  of  PCBs  into  25  pm  PE  by  measuring  the  fractional  loss  of  13C-labeled 
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PCBs  after  44  days  of  field  deployment.  The  average  fractional  equilibration  of  penta-  and 
hexachlorobiphenyls  from  different  locations  in  the  field  were  only  37%  and  33%,  respectively. 


Effect  of  vibration  on  polymer  uptake  rate.  Periodic  vibration  of  the  PE  sampler  resulted  in 
faster  uptake  compared  to  static  deployment  for  all  PAHs  measured.  As  shown  in  Figure  5.3,  the 
uptake  of  chrysene  after  56  days  was  only  40%  of  equilibrium  in  static  exposure  compared  to 
100%  for  the  vibrating  system.  Even  after  7  days  of  exposure,  the  vibrating  system  reaches  63% 
of  equilibrium  for  chrysene  compared  to  20%  in  the  static  system.  Comparison  between  the 
static  and  vibrating  system  for  several  other  PAH  compounds  is  shown  in  Figure  5.4  and  SI 
Figure  S5.4.  Both  PRCs  were  lost  more  than  90%  at  the  first  time  point  of  measurement  for  the 
vibrating  systems  and  well-mixed  systems  (SI  Figure  S5.5),  which  generally  agreed  with  the 
uptake  of  the  analogous  PAHs.  As  PRC  performance  was  not  the  focus  of  this  study,  using  a  full 
range  of  high  molecular  weight  PRCs  and  performing  shorter  duration  measurements  of  PRC 
loss  are  recommended  in  future  work  in  order  to  provide  a  more  robust  analysis  of  PRC 
performance  in  the  vibrating  system. 

For  PAH  compounds  less  hydrophobic  than  chrysene  (e.g.  fluoranthene  and  pyrene),  the  static 
system  showed  reasonable  uptake  of  close  to  70%  in  56  days  while  the  vibrating  system  reached 
equilibrium  during  that  exposure.  When  the  passive  sampler  is  close  to  equilibrium,  the 
correction  for  non-equilibrium  is  relatively  accurate.  However,  as  the  kinetics  slow  down  and 
only  a  small  fraction  of  equilibrium  is  achieved  in  the  sampler,  the  correction  for  non¬ 
equilibrium  becomes  error-prone.9  This  is  the  case  for  the  PAH  compounds  that  are  more 
hydrophobic  than  chrysene.  The  improvement  over  static  system  was  more  evident  for  larger 
molecular  weight  PAHs.  For  example,  the  measured  fractional  uptake  of  benzo(a)pyrene  was 
improved  from  6%  to  55%  in  7  days  and  from  23%  to  90%  after  56  days.  The  fact  that  vibration 
was  more  effective  on  uptake  rate  of  heavier  PAHs  can  be  explained  based  on  the  nature  of  the 
overall  mass  transfer  resistance  in  a  passive  sampler6 


KPEwDPE 


Where  ko  is  the  overall  mass  transfer  rate  constant  and  8  is  the  diffusion  distance. 


The  first  term  on  the  right  hand  side  of  Equation  15,  describes  the  sediment-side  mass  transfer 
resistance.  The  large  Kpew  values  of  the  more  hydrophobic  compounds  makes  the  second  term 
smaller  than  the  first,  thereby  making  the  first  term  or  the  water  side  mass  transfer  more 
dominant.  The  ratio  of  diffusivities  in  water  and  polymer  (Dw/Dpe)  also  increases  as  the 
compound  becomes  more  hydrophobic.  However,  Kpew  tends  to  be  the  dominant  factor 
influencing  which  side  controls  mass  transfer  since  Kpew  changes  over  a  wider  range  of  values 
compared  to  the  range  where  the  diffusivity  ratios  vary  for  different  compounds.27  As  a  result, 
disrupting  the  depletion  layer  (decreasing  apparent  S)  by  means  of  vibration  will  be  more 
effective  on  the  overall  mass  transfer  rate  of  larger  molecular  weight  PAHs.  Our  experimental 
data  indicated  that  concentration  of  larger  molecular  weight  PAHs  (log  Kow>5.5)  reached  more 
than  50%  of  their  equilibrium  concentration  in  PE  within  14  days  of  deployment  in  the  vibrating 
system  (SI  Figure  S5.4).  Thus,  field  deployment  of  passive  samplers  with  periodic  vibration  will 
encounter  less  challenges  of  non-equilibrium  correction  for  larger  molecular  weight  compounds. 
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Results  from  the  vibration  work  demonstrated  a  great  improvement  of  PE  uptake  after  short 
exposure  times  (1-2  weeks).  The  faster  exchange  will  also  allow  the  use  of  high  molecular 
weight  PRCs  to  accurately  correct  for  the  remaining  non-equilibrium  in  a  vibrating  system 
without  the  need  for  long  exposure  times. 

Comparison  of  three  modeling  approaches  in  static,  vibration  and  well-mixed  deployments. 
Figure  5.3-a  and  5.3-b  demonstrate  modeling  results  for  chrysene  in  static  and  vibrating  systems 
based  on  local  equilibrium,  slow  desorption,  and  fast  desorption  models.  As  shown  in  Figure  5.3- 
a,  in  static  deployment,  overall  mass  transfer  into  the  polymer  is  slow  and  the  predictions  based 
on  local  equilibrium  assumption  and  the  fast  desorption  rate  model  are  indistinguishable.  Also, 
both  predictions  are  close  to  the  observed  uptake  of  chrysene  in  the  static  deployment.  However, 
the  slow  rate  of  desorption  model  predicts  uptake  that  is  slower  than  the  observed  values.  Thus, 
it  appears  that  for  modeling  static  deployments,  due  to  the  rate  limiting  mass  transfer  through  the 
depletion  layer,  it  is  adequate  to  assume  local  equilibrium  between  sediment  and  water. 

For  the  vibrating  system,  assuming  that  the  sediment  is  effectively  mixed  after  vibration,  the 
model  predictions  based  on  local  equilibrium  and  fast  rate  of  desorption  deviated  as  shown  in 
Figure  5.3-b.  While  the  fast  rate  of  desorption  model  prediction  is  close  to  the  measured  uptake 
of  chrysene,  the  prediction  based  on  local  equilibrium  greatly  over-predicts  uptake.  The  same 
results  were  observed  for  pyrene  as  indicated  in  SI  Figure  S5.6.  Thus,  it  appears  that  as  the 
depletion  layer  is  disrupted  by  vibration  and  mass  transfer  is  enhanced,  desorption  from  sediment 
becomes  limiting  and  local  equilibrium  between  sediment  and  water  can  no  longer  be  assumed. 
The  fast  desorption  rate  adequately  predicts  uptake  while  the  slow  desorption  rate  still  under¬ 
predicts  uptake  in  the  polymer.  Past  work  has  demonstrated  that  PAH  desorption  from  MGP- 
impacted  sediments  is  characterized  by  slow  and  fast  desorbing  PAH  fractions. 17  However,  the 
overall  kinetics  for  the  vibration  system  appear  to  be  driven  by  the  fast  desorption  kinetics  and 
not  the  slow  desorption  likely  because  the  sediments  are  not  being  depleted  enough  to  reach  the 
slow  desorption  regime. 

For  the  well-mixed  system,  a  local-equilibrium  model  was  implemented.  As  shown  in  Figure 
5.3-c,  the  model  agrees  reasonably  with  the  observed  fast  uptake  of  chrysene  and  the  attainment 
of  equilibrium  in  the  matter  of  a  few  days.  The  well-mixed  system  brings  a  large  volume  of 
sediment  to  contribute  to  the  required  uptake  by  the  passive  sampler  and  hence  the  kinetics  are 
fast  enough  to  appear  close  to  instantaneous  equilibrium. 
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Figure  5.3-  Fractional  uptake  of  chrysene  in  PE  passive  sampler  in  static  (a),  vibrating  (b),  and  well- 
mixed  (c)  systems.  Experimental  data  have  been  shown  by  circles,  and  model  simulations  have  been 
shown  by  lines.  Note  that  for  the  static  system  the  local  equilibrium  and  fast  desorption  models  overlap. 


Modeling  uptake  of  pyrene,  chrysene,  and  benzol ajpyrene.  Experimental  and  modeling  results 
for  pyrene  (log  Kow  =  4.9),  chrysene  (log  Kow  =  5.7)  and  benzo(a)pyrene  (log  Kow  =  6.1)  in  static, 
vibrating,  and  well-mixed  systems  are  shown  in  Figure  5.4.  Based  on  the  discussion  above,  only 
the  fast  desorption  model  was  used  for  predicting  uptake  in  static  and  vibrating  systems,  and  a 
local  equilibrium  model  was  used  for  predicting  the  well-mixed  system.  As  shown  in  Figure  5.4, 
the  model  predictions  reasonably  agreed  with  the  uptake  profiles  of  the  3  PAHs  in  all  three 
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modes  of  exposure.  For  all  three  PAHs,  although  there  were  small  differences  between  the 
experimental  results  from  2  min  and  5  min  pause  times  (not  statistically  significant  except  for 
pyrene  at  day  14),  the  model  predicted  that  there  should  not  be  a  significant  difference  between 
the  two  (see  SI  Tables  S5.4,  S5.5,  and  S5.6).  In  the  model  we  assumed  that  both  pulse  durations 
are  long  enough  to  mix  up  the  sediment  to  initial  bulk  concentration.  However,  in  reality  short 
pulse  time  (2  seconds)  in  the  low  frequency  system  may  not  be  long  enough  to  satisfy  the  model 
assumption. 


Figure  5.4-  Fractional  uptake  of  pyrene  (a),  chrysene  (b),  and  benzo(a)pyrene  (c)  in  PE  passive  sampler  in 
four  differently  exposed  systems.  Experimental  data  have  been  shown  by  symbols,  and  fast  desorption 
model  simulations  have  been  shown  by  lines. 
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To  better  illustrate  how  the  periodic  vibration  impacts  development  of  the  sediment-side 
depletion  layer,  the  concentration  profiles  in  the  sediment  and  half  width  of  the  PE  sampler  was 
plotted  for  chrysene  as  a  function  of  exposure  time.  As  shown  in  Figure  5.5  left  panel  for  the 
static  deployment,  chrysene  in  the  sediment  side  is  depleted  well  into  300  pm  distance  from  the 
PE  surface  after  20  days.  The  concentration  values  for  PE  plotted  in  Figure  5.5  have  been 
multiplied  with  Ksed-PE  to  make  the  values  directly  comparable  to  the  sediment  concentration.  As 
indicated  in  Figure  5.5,  left  panel,  the  concentration  in  PE  reaches  equilibrium  with  the  depleted 
concentration  in  the  sediment  adjacent  to  the  surface  in  2  days  or  less  and  the  mass  transfer 
limitation  is  moved  to  the  sediment- side  depleted  layer.  Also,  after  the  first  few  hours  of 
exposure,  there  is  no  concentration  gradient  within  the  polymer  indicating  all  mass  transfer 
resistance  is  in  the  sediment  side.  In  stark  comparison,  for  the  periodic  vibration  deployment, 
since  the  sediment  side  is  mixed  up  at  frequent  intervals,  the  deep  depletion  layer  in  the  sediment 
side  is  not  able  to  develop  and  the  mass  transfer  resistance  is  limited  to  a  few  microns  near  the 
polymer  surface.  In  fact,  as  seen  in  the  2-day  simulation  for  the  concentration  gradient  within  PE 
in  the  vibrating  system,  the  model  suggests  that  there  is  still  some  polymer-side  resistance.  As  a 
result,  the  concentration  in  PE  rises  much  more  rapidly  than  in  the  static  case. 
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Figure  5.5-  Model  simulations  of  chrysene  concentration  profile  within  sediment  and  PE  for  static  and 
periodic  vibration  deployments.  The  values  shown  on  the  Y-axis  are  chrysene  concentration  in  sediment 
and  concentration  in  PE  multiplied  by  sediment-PE  partition  constant  (Cpe  x  Kseci-  pe). 
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For  the  well-mixed  system,  pyrene  uptake  is  fast  and  modeled  reasonably  well  by  the  well-mixed 
local-equilibrium  model  (Figure  5.4-a).  However,  the  model  is  not  accurate  for  early  times  for 
the  more  hydrophobic  compounds,  especially  benzo(a)pyrene  (Figure  5.4-c).  This  could  be  due 
to  inaccuracy  of  reported  diffusivity  in  PE,  slow  desorption  from  sediment,  or  inefficient  mixing 
of  the  system.  Inefficient  mixing  slows  the  uptake  of  more  hydrophobic  compounds  since 
sediment-side  mass  transfer  is  dominant  for  high  molecular  weight  PAHs  (as  described  earlier) 
and  disruption  of  depletion  layer  by  efficient  mixing  is  more  important.  Previous  studies22  have 
observed  slower  uptake  of  larger  PCBs  in  well  mixed  sediment- slurry  systems  and  attributed  the 
slow  kinetics  to  inefficient  mixing  and  slow  desorption. 

Optimization  of  vibration  frequency.  Results  of  the  vibration  model  for  PE  uptake  of  chrysene 
were  used  to  test  the  effect  of  different  vibration  frequencies  and  to  optimize  power  requirement. 
All  vibration  frequencies  had  pulse  duration  of  5  seconds  and  varied  in  pause  time.  The  modeling 
results  were  based  on  the  fast  desorption  model  and  pause  times  were  set  at:  2  min,  5  min,  6 
hour,  1  day,  and  5  days  (Figure  5.6).  Small  differences  (6%  or  less)  were  observed  between  the 
2-min  and  6-hour  pause  models  (see  SI  Table  S5.7).  Even  the  1-day  pause  model  predicted  an 
uptake  profile  which  was  not  greatly  different  from  the  2  min-pause  model.  The  5-day  pause 
profile  tracked  the  unmixed  profile  for  the  first  5  days,  then  jumped  up  to  a  higher  uptake  profile 
as  the  mixing  altered  the  boundary  condition.  Although  the  5  day-pause  model  does  not  show 
significant  improvement  over  the  static  system  in  short  deployment  times  (<7  days),  the  model 
still  predicts  75%  fractional  uptake  after  28  days.  In  comparison  the  fractional  uptake  of 
chrysene  is  only  40%  in  static  system  for  the  same  exposure  duration.  These  results  indicate  that 
significant  improvement  over  the  static  system  is  possible  with  lower  frequencies  of  vibration 
(large  pause  times).  Increasing  pause  time  is  desirable  as  there  is  no  energy  consumption  during 
the  pause  time  when  the  motor  is  not  vibrating.  In  order  to  optimize  pause  time,  energy 
consumption  was  estimated  for  deploying  a  vibration  system  for  15  days  and  a  range  of  pause 
times.  The  fractional  uptake  of  three  PAHs  (pyrene,  chrysene  and  benzo(a)pyrene)  in  PE  after  15 
days  was  also  determined  with  the  model  for  each  pause  time  and  plotted  versus  energy 
consumption  (SI  Figure  S5.7).  Only  30m Wh  is  required  to  enhance  the  fractional  uptake  of 
chrysene  to  80%  in  15  days  of  deployment.  In  comparison,  a  typical  AA  size  NiMH  rechargeable 
battery  can  provide  about  1000  mWh  of  energy.  However,  further  enhancement  of  fractional 
uptake  in  PE  from  80%  to  82%  increases  the  energy  requirement  to  2340mWh.  This  is  also 
evident  from  Figure  5.6  where  we  see  little  enhancement  of  uptake  when  the  pause  time  is 
reduced  below  6  hours. 
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Figure  5.6-  Model  simulation  of  chrysene  uptake  in  vibrating  system  with  different  pause  times  of 
vibration.  Experimental  data  have  been  shown  by  symbols,  and  fast  desorption  model  simulations  have 
been  shown  by  lines. 


5.4  Implications 

A  recent  review  article  on  passive  sampling  by  Booij  et  a\.23  concluded  that  options  to  reduce 
time  for  equilibrium  are  limited  to  manipulation  of  area/volume  ratio  of  sampler,  choice  of 
sampler  material,  and  flow  rate.  Past  work  has  led  to  significant  optimization  of  the  physical 
aspects  of  the  passive  samplers,  yet  attainment  of  equilibrium  remains  difficult  for  strongly 
hydrophobic  compounds  especially  for  in  situ  sediment  pore  water  measurements.  Also,  there  is 
little  that  can  be  done  in  the  sediment  environment  to  enhance  porewater  velocity.  This  study 
advances  the  practice  of  passive  sampling  by  addressing  a  key  bottleneck  through  the  novel 
introduction  of  periodic  vibration  in  the  sampling  platform  to  disrupt  the  depletion  layer  that 
develops  in  the  sediment  side  and  slows  mass  transfer.  Through  mass  transfer  modeling,  we  also 
provide  a  mechanistic  interpretation  of  how  periodic  vibration  enhances  the  approach  to 
equilibrium.  While  PRC  corrections  have  allowed  extension  of  passive  sampling  to  compounds 
that  do  not  achieve  equilibrium  during  a  reasonable  period  of  deployment,  the  introduction  of 
vibration  greatly  enhances  approach  to  equilibrium,  reduces  deployment  times,  and  extends  the 
use  of  passive  sampling  in  conjunction  with  PRCs  to  strongly  hydrophobic  compounds  as  shown 
for  benzo(a)pyrene  in  this  study.  Further  development  is  necessary  to  build  robust  prototypes 
and  perform  tests  in  the  field.  The  size  of  the  devices  and  motors  may  also  need  to  be  increased 
to  enable  deployment  of  larger  passive  samplers  to  target  low  concentrations  in  sediment. 
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6.0  RESEARCH  OBJECTIVE  #2:  Test  the  in- situ  vibration  approach 
for  strongly  hydrophobic  compounds  like  dioxins  and  furans. 


6.1  Introduction 

In  the  previous  chapter,  we  proposed  using  a  vibrating  platform  for  the  deployment  of  passive 
sampling  devices  in  sediments.  We  demonstrated  through  laboratory  measurements  and 
numerical  modeling  that  the  proposed  platform  greatly  enhanced  the  rate  of  mass  transfer  of 
PAHs  into  PE  passive  samplers.  In  this  chapter  we  present  an  extension  of  the  work 
demonstrating  effectiveness  for  strongly  hydrophobic  chlorinated  organics  in  sediments.  As 
shown  in  Figure  6.1,  most  polychlorinated  dibenzofurans  (PCDFs)  and  polychlorinated 
dibenzodioxins  (PCDDs)  have  the  same  range  of  Kows  (6.5<  Fog  Kow<8)  as  hexa-,  hepta-,  and 
octachloro-  PCBs.  Thus,  for  this  SEED  project,  the  most  abundant  hydrophobic  PCB  congeners 
were  used  as  analogs  for  dioxins  and  furans.  A  sediment  sample  from  a  DoD  site  containing 
highly  chlorinated  PCBs  was  used  in  laboratory  mesocosms  to  test  the  effectiveness  of  the 
vibrating  platform  on  uptake  rate  of  hydrophobic  PCB  congeners  into  PE  passive  samplers.  In 
this  section,  the  experimental  procedure  and  results  of  PCB  uptake  in  PE  under  static,  vibrating, 
and  well  mixed  conditions  are  presented. 
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Figure  6.1.  Comparison  of  Log  Kows  of  PCBs  (obtained  from  Hawker  and  Connell,  1988)  and  PCDD/Fs 
(obtained  from  Sacan  et  al.,  2005) 
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6.2  Materials  and  methods 

Materials.  Low  density  polyethylene  (LDPE)  sheets  (25pm  thickness),  manufactured  by  Poly- 
America  (Grand  Prairie,  TX,  USA)  were  purchased  from  the  Home  Depot  (local  store  in 
Baltimore,  MD).  PCB  stock  solutions  were  purchased  from  Fisher  Scientific  (Pittsburg,  PA, 
USA).  Cylindrical  vibrating  motors  with  a  diameter  of  24  mm  and  length  of  31  mm  were 
purchased  from  Precision  Microdrives  (London,  United  Kingdom).  The  motors  operated  at  the 
rated  voltage  of  3  V  and  operation  current  of  190  mA.  Their  rated  vibrating  speed  and  typical 
normalized  amplitude  of  vibration  were  5000  rpm  and  10  G,  respectively.  Pulse-pause  timers 
(model  60H)  used  to  control  the  vibration  motors  were  purchased  from  Velleman  Inc.  (Fort 
Worth,  TX).  Prior  to  use,  PE  sheets  were  soaked  twice  in  hexane/  acetone  (50/50)  and  left  on  a 
shaker  for  24  h  each  time  to  remove  oligomers  and  any  target  and  non-target  contaminants. 

Clean  PE  sheets  were  cut  into  20  mg  (2cmx5cm)  and  40  mg  (4cmx5cm)  strips.  PE  strips  were 
then  soaked  in  a  80:20  methanol/water  solution  containing  four  PRCs  (2,4,5-trichlorobiphenyl 
(PCB  29),  2,3',4,6-tetrachlorobiphenyl  (PCB  69),  2,2',4,4',6,6'-hexachlorobiphenyl  (PCB  155), 
and  2,3,3',4,5,5',6-heptachlorobiphenyl  (PCB  192)).  The  PRC  solution  was  allowed  to  equilibrate 
for  15  days  on  an  orbital  shaker.  After  impregnation  with  PRCs,  the  strips  were  transferred  into  a 
DI  water  jar  and  left  on  an  orbital  shaker  overnight  to  remove  the  methanol  from  PE  strips.  Upon 
removal  from  water,  three  strips  were  extracted  immediately  in  hexane  (3  x  24  h)  to  determine 
the  initial  PRC  concentration  on  PE  strips. 

Source  of  Sediments.  PCB  contaminated  sediment  samples  from  Site  102  Abraham’s  Creek  were 
used  in  this  study.  The  site  is  located  in  the  Marine  Corps  Combat  Development  Command 
(MCCDC)  military  training  facility  in  Quantico,  Virginia,  35  miles  south  of  Washington,  D.C. 

In  Situ  Shaken  Passive  Sampler  Design.  A  platform  made  of  copper  was  fabricated  with 
assistance  from  the  UMBC  machine  shop  (Figure  S6.4).  The  platform  consisted  of  a  copper  pipe 
section,  inside  which  the  motor  was  located.  Four  fins  made  of  copper  plates  and  meshes  were 
attached  radially  to  the  copper  pipe.  The  copper  wire  mesh  fins  were  designed  as  pockets  that 
could  hold  PE  sheets  sandwiched  within.  A  20  mg  PE  strip  (prepared  as  discussed  above)  was 
located  inside  each  pocket  and  the  open  edge  was  sewn  with  copper  wire. 

Two  platforms  were  connected  in  parallel  to  a  timer  and  a  power  supply  (2  rechargeable  batteries 

1.2  V,  700  mAh  each).  The  timers  were  powered  by  a  12  V  power  supply  and  were  programmed 
to  control  motor  vibration  duration  and  frequency.  Two  experiments  were  conducted  in  order  to 
measure  the  uptake  rate  of  PCBs  into  PE  under  different  deployment  modes.  In  the  first 
experiment  (experiment  1),  PE  sheets  were  deployed  in  sediment  under  static,  well-mixed,  and 
high  frequency  vibrating  modes.  In  the  second  experiment  (experiment  2),  PE  sheets  were 
deployed  under  static  and  low  frequency  vibrating  modes.  The  pulse  duration  of  vibrating  motors 
was  5  s  and  the  pause  period  was  2  min  and  5  d  in  the  high  and  low  frequency  vibrating  systems, 
respectively.  Details  of  experimental  set  up  are  as  follows: 

Experiment  1-  The  sediment  was  mixed  with  DI  water  containing  200  mg/L  sodium  azide  to 
make  a  slurry  with  weight  ratio  of  1:2  (dry  sediment/water).  For  the  well-mixed  exposure,  three 
wide  mouth  jars  with  Teflon-lined  lids  were  prepared.  The  sediment-water  slurry  (400  mL)  and 
four  of  the  20  mg  PE  samplers  were  transferred  into  each  jar.  The  jars  were  placed  on  a  rotary 
agitator  and  tumbled  with  the  speed  of  28  rpm.  The  rest  of  the  sediment  slurry  was  placed  in  two 
large  glass  trays  (60  cmxlO  cmx  10  cm).  Three  motors  were  placed  inside  the  sediment  in  the 
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trays  (two  of  the  motors  were  placed  in  one  tray  with  10  cm  separation  to  prevent  influencing 
each  other).  The  timer  was  programmed  to  vibrate  the  motors  with  the  frequency  of  2  min  pause 
and  5  s  pulse.  Twelve  additional  PE  samplers  with  a  mass  of  40  mg  each  were  enclosed  in 
stainless  steel  mesh  without  motors  and  were  placed  inside  the  same  trays  in  a  static  mode  at 
least  10  cm  away  from  the  vibrating  motors  and  5  cm  apart  from  each  other  to  simulate  a  static 
system.  The  trays  were  then  covered  with  aluminum  foil.  PE  strips  were  removed  from  the  well- 
mixed,  static,  and  vibrating  systems  and  analyzed  for  PCBs  and  PRC  compounds  after  7,  14,  28 
and  56  days.  At  each  time  point,  three  strips  were  collected  from  each  system  for  extraction  and 
analysis. 

Experiment  2-  PE  deployment  under  static  and  vibrating  modes  was  repeated  in  a  slurry  made 
with  fresh  portion  from  the  same  sediment  as  discussed.  This  time  the  motors  vibrated  with  the 
frequency  of  5  s  pulse  and  5  d  pause  by  manual  connection  to  the  power  supply.  PE  strips  were 
collected  in  triplicates  from  static  and  vibrating  systems  and  analyzed  for  PCBs  after  7,  14,  28 
and  56  days. 

PCB  Extraction  and  Analysis.  Upon  removal  from  the  sediment,  PE  strips  were  rinsed  with  water 
and  wiped  with  laboratory  tissue  to  remove  water  and  adhering  particles.  Prior  to  extraction 
PCB  congeners  BZ#14  and  65  were  added  as  surrogate  standards  to  assess  the  effectiveness  of 
sample  processing,  and  results  with  lower  than  80%  surrogate  recoveries  were  discarded. 
Samplers  were  extracted  with  hexane  (3  x  24  h,  with  sequential  extracts  pooled).  The  final 
extraction  volumes  were  blown  down  to  lmL  using  a  stream  of  nitrogen  gas.  Sediment  PCBs 
were  extracted  by  sonication  (EPA  method  3550B).  Final  extracts  from  passive  samplers  and 
sediment  were  cleaned  using  deactivated  silica  gel  (EPA  method  3630C).  Two  internal  standards 
were  added  to  the  final  extracts  before  analysis  (PCB  congeners  BZ#30  and  204).  PCBs  were 
analyzed  in  an  Agilent  6890N  Gas  Chromatograph  with  an  electron  capture  detector  (ECD) 
based  on  analytical  method  described  in  Beckingham  and  Ghosh.12 

PCB  Desorption  Study.  The  desorption  rates  of  PCBs  from  sediment  particles  into  water  were 
measured  based  on  procedures  in  Ghosh  et  al.13  Briefly,  Tenax  beads  (0.5  g)  and  sediment 
sample  (5  g)  were  added  to  12-mL  glass  vials.  DI  water  (10  mL)  containing  sodium  azide  (1  g/1) 
was  added  to  each  vial.  The  vials  were  palced  on  a  rotary  agitator  and  tumbled  with  the  speed  of 
28  rpm.  At  sampling  times  of  2  h,  4  h,  8  h,  12h,  and  48  h  the  sediment  was  allowed  to  settle  and 
Tenax  beads  to  float  up.  Tenax  beads  were  scooped  out  from  each  vial  and  a  fresh  portion  was 
added.  The  Tanax  beads  were  dried  with  anhydrous  sodium  sulfate  and  extracted  by  adding  three 
volumes  of  15  mL  50/50  (v/v)  hexane-acetone  in  a  20-mL  vial  and  shaking  horizontally  on  a 
rotary  shaker  for  24  h.  Prior  to  extraction  PCB  congeners  BZ#14  and  65  were  added  as  surrogate 
standards  to  assess  the  effectiveness  of  sample  processing,  and  extracts  with  lower  than  80% 
surrogate  recoveries  were  discarded.  The  final  extraction  volumes  were  blown  down  to  lmL 
using  a  stream  of  nitrogen  gas  and  were  cleaned  using  deactivated  silica  gel  (EPA  method 
3630C).  The  extracts  were  then  analyzed  for  PCBs. 

A  two-compartment  model  (Ghosh  et  al.13)  was  used  to  describe  the  desorption  kinetics  of  PCBs 
from  sediment  particles  into  porewater  (see  Supporting  Information  for  more  details).  The  model 
was  fitted  to  the  normalized  desorption  data  for  PCB  128,  183,  and  194  in  order  to  find  the 
desorption  rate  constants  (Supporting  Information  Figure  S6.5). 
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PRC  Correction  for  Non-equilibrium.  We  used  a  graphical  user  interface  (GUI,  Tcaciuc  et  al.14) 
to  calculate  the  extent  of  equilibrium  of  target  PCBs  from  PRC  losses.  This  GUI  is  developed 
based  on  the  mass  transfer  model  presented  in  Fernandez  et  al.6  PRC  losses  in  static  and 
vibrating  systems  were  used  to  determine  the  extent  of  equilibrium  of  target  PCBs  in  the 
corresponding  systems.  PCB  and  PRC  concentration  measurements  were  conducted  in 
triplicates  and  PRC  losses  from  each  PE  strip  were  used  to  calculate  the  extent  of  equilibrium  of 
PCBs  in  the  same  strip. 

Modeling  Uptake  of  Analytes  from  Sediment  Porewter.  Mass  transfer  models  described  in 
Jalalizadeh  and  Ghosh1  were  employed  to  describe  PCB  diffusion  from  sediment  particles  into 
porewater,  and  from  porewater  into  the  polymer  under  static,  vibrating,  and  well-mixed  modes. 
Briefly,  for  static  deployment  mass  transfer  within  the  polymer  and  in  porewater  was  described 
by  Fick’s  2nd  law  of  diffusion.  Instantaneous  equilibrium  was  assumed  at  the  polymer  surface 
boundary  with  porewater.  PCB  desorption  from  sediment  particles  into  porewater  followed  a 
first  order  kinetics,  assuming  all  PCBs  are  associated  with  the  fast  desorbing  pool  in  sediment. 
The  diffusion  model  with  the  assumption  of  instantaneous  equilibrium  between  sediment 
particles  and  porewater  (local  equilibrium  model)  adequately  predicted  the  PE  uptake  in  a  static 
system  in  some  previous  work  (Fernandez  et  al.,  2009).  However,  local  equilibrium  cannot  be 
assumed  in  a  vibrating  system  (Jalalizadeh  and  Ghosh1)  because  mass  transfer  from  sediment 
becomes  limiting  due  to  disruption  of  the  depletion  layer  by  vibration.  In  our  previous  work 
(Jalalizadeh  and  Ghosh1)  we  showed  that  model  predictions  based  on  local  equilibrium  and  fast 
rate  of  desorption  were  identical  in  a  static  system.  However,  the  simulations  deviated  for  the 
vibrating  system.  We  also  showed  that  the  overall  kinetics  of  PAHs  for  the  vibrating  system  were 
driven  by  the  fast  desorption  kinetics  and  not  the  slow  desorption,  likely  because  the  sediments 
are  not  depleted  enough  to  reach  the  slow  regime.  The  fast  desorption  rate  constants  determined 
from  PCB  desorption  study  were  used  for  model  simulation. 

The  same  modeling  structure  for  the  static  deployment  was  employed  for  the  vibrating  system 
when  motors  are  in  the  pause  mode.  We  assumed  that  vibration  mixes  up  and  increases  the 
sediment  concentration  at  the  vicinity  of  the  polymer  to  the  initial  concentration  in  sediment. 

This  assumption  was  validated  by  visual  observation  of  the  vibrating  system  in  our  previous 
work  (Jalalizadeh  and  Ghosh1),  which  indicated  that  the  extent  of  fluidized  mixing  by  vibration 
is  at  least  1  cm  from  the  surface  of  the  device.  Vibrating  motors  in  the  previous  work  had  smaller 
amplitude  of  vibration  (6  G)  and  the  sediment  slurry  had  similar  porosity  and  density  compared 
to  the  current  study. 

For  the  well-mixed  system  mass  transfer  was  assumed  to  be  limited  by  the  polymer  side.  Since 
sediment  is  perfectly  mixed,  porewater  concentration  remains  constant  and  equal  to  the  initial 
value  during  deployment  time.  In  addition,  local  equilibrium  was  assumed  between  sediment 
particles  and  porewater  as  mixing  brings  large  volume  of  sediment  to  the  vicinity  of  the  polymer 
and  kinetics  are  fast  enough  to  appear  close  to  instantaneous  equilibrium. 

Model  equations  were  solved  in  Matlab  using  an  explicit,  finite-difference  numerical  modeling 
technique  (Crank15).  The  equations  and  the  details  of  parameter  estimation,  numerical  solution, 
and  Matlab  codes  are  indicated  in  the  Supporting  Information. 
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6.3-  Results  and  discussion 

PCB  Concentration  in  Sediment.  As  shown  in  Figure  6.2,  hexa-,  hepta-,  and  octachloro-  PCBs 
were  the  most  abundant  homologs  in  the  sediment.  Mono-  and  dichloro-  PCBs  were  not  detected, 
and  the  total  PCB  concentration  was  88  ng/g  in  the  sediment.  This  sediment  sample  was  found  to 
be  suitable  for  the  present  study  because  it  had  a  relatively  low  level  of  PCBs  (pushing  our 
analytical  effort)  and  our  target  homologs  were  dominant  in  sediment  (mimicking  PCDD/Fs). 
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Figure  6.2.  PCB  Homolog  distribution  in  sediment  based  on  the  mean  of  three  measurements 


Effect  of  Vibration  with  High  Frequency  on  PCB  Uptake  into  PE.  The  uptake  profiles  of  six 
PCB  congeners  (each  representing  a  homolog  group)  are  compared  in  static,  vibrating  and  well- 
mixed  systems  in  Figure  6.3.  Since  measurements  of  PE  uptake  in  static  system  from  experiment 
1  and  2  were  not  significantly  different  for  all  PCBs  (Supporting  Information  Table  S6.4),  only 
the  results  from  experiment  1  are  shown  for  the  static  deployment.  The  concentrations  reported 
are  in  ng/g  PE.  PCBs  had  very  slow  uptake  rates  into  PE  in  the  static  exposure  system.  Periodic 
vibration  of  the  PE  sampler  with  high  frequency  resulted  in  much  faster  uptake  compared  to 
static  deployment  for  all  congeners.  For  example,  after  28  days  of  deployment,  the  fractional 
uptake  of  PCB  99,  PCB  132,  and  PCB  (87+182)  were  37%,  26%,  and  22%,  respectively  in  the 
static  system  compared  to  100%,  85%,  and  72%  in  the  vibrating  system.  All  PCBs  reached  to  95 
-  100%  of  equilibrium  after  56  days  of  deployment  in  the  high  frequency  vibrating  system.  This 
is  while,  none  of  the  congeners  reached  to  more  than  50%  of  equilibrium  in  static  deployment  for 
the  same  period.  For  example,  the  uptake  of  PCB  201  after  56  days  was  only  30%  of  equilibrium 
in  static  exposure  compared  to  100%  in  the  vibrating  system.  Even  after  7  days  of  exposure,  the 
vibrating  system  reaches  35%  of  equilibrium  for  PCB  201  compared  to  16%  in  the  static  system. 
Thus,  with  vibration  it  is  apparent  that  in  situ  measurements  of  the  full  range  of  PCB  congeners 
in  sediments  would  approach  very  high  fraction  equilibrium  allowing  more  accurate  estimation 
of  the  existing  porewater  concentration  in  sediment. 
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In  the  well-mixed  system,  no  change  in  the  PE  concentration  was  observed  for  any  of  the 
congeners  after  28  days  (PE  concentrations  after  28  and  56  days  of  exposure  were  not 
significantly  different  as  shown  in  Supporting  Information  Table  S6.2).  Thus,  PE  concentrations 
from  fully  mixed  system  after  56  days  of  exposure  were  chosen  as  equilibrium  concentrations  for 
all  PCB  congeners. 
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Figure  6.3.  Comparison  of  PCBs  uptake  into  PE  in  static,  vibrating,  and  fully  mixed  systems 


Effect  of  Vibration  with  Low  Frequency  on  PCB  Uptake  into  PE.  The  fractional  uptake  of  all 
PCBs  in  the  low  frequency  vibrating  system  improved  within  a  factor  of  1.5  and  2  over  the  static 
system  after  28  days  (Figure  6.2,  and  Supporting  Information  Figures  S6.1-S6.3).  However,  all 
PCBs  reached  to  only  30%  to  60%  of  equilibrium  in  the  vibrating  system,  except  for  PCB  31, 
PCB  70+76,  and  PCB  185  that  reached  to  approximately  70%.  Even  after  56  days  the  fractional 
uptake  of  all  PCBs  were  between  40%  and  72%  (77%  for  PCB  185).  These  results  indicated  that 
the  motors  need  to  vibrate  with  pause  durations  less  than  5  days  in  order  to  improve  the  uptake  of 
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larger  molecular  weight  PCBs  (hexa-,  hepta-,  and  octachloro-congeners)  to  more  than  60%  in  28 
days. 

Estimation  of  PCB  desorption  rate  constants.  Modeling  of  uptake  rates  into  PE  were  conducted 
for  PCB  128  (2,2',3,3',4,4'-hexachlorobiphenyl;  i0g  Kow  =  6.74),  PCB  183  (2, 2', 3, 4, 4', 5', 6- 
heptachlorobiphenyl;  log  Kow  =  7.2),  and  PCB  194  (2,2',3,3',4,4',5,5'-octachlorobiphenyl;  log 
Kow  =  7.8).  Thus,  the  fast  (kf)  and  slow  (ks)  desorption  rate  constants  are  only  reported  for  the 
three  mentioned  PCBs  (Supporting  Information  Table  S6.8)  and  only  Kf  values  were  used  for 
model  simulation. 

The  measured  Kfs  for  three  PCBs  were  approximately  one  order  of  magnitude  larger  compared 
to  the  reported  rates  by  Zimmerman  et  al.16  The  authors  measured  the  desorption  rate  of  PCBs 
from  an  untreated  sediment  collected  from  the  intertidal  region  of  South  Basin  at  Hunters  Point. 
Kf  values  from  this  study  were  also  two  orders  of  magnitude  larger  compared  to  the  rates 
measured  by  Ghosh  et  al.17  in  field-contaminated  sediment  obtained  from  Aluminum 
Corporation  of  America  (at  25  °C). 

Modeling  the  Uptake  of  PCBs.  Experimental  and  modeling  results  for  PCB  128,  PCB  183,  and 
PCB  194  in  static,  vibrating  and  well-mixed  systems  are  indicated  in  Figure  6.4.  As  shown  in  the 
Figure,  modeling  results  had  good  agreement  with  the  uptake  rate  of  three  PCBs  in  static  and 
vibrating  systems  with  low  and  high  frequency  of  vibration. 

The  well-mixed  model  over  predicted  the  earlier  time  point  measurements  (Figure  6.3),  likely 
due  to  slow  desorption  of  PCBs  from  sediment  particles  or  inefficient  mixing  of  the  sediment. 
Slower  kinetics  of  larger  hydrophobic  chemicals  in  well-mixed  system  has  been  also  observed  in 
previous  studies  (Jalalizadeh  and  Ghosh1;  Arp  et  al.18).  Deviations  between  the  model  simulation 
and  experimental  data  were  more  obvious  for  more  hydrophobic  PCBs  (PCB  183  and  PCB  194), 
as  mass  transfer  of  these  compounds  is  mostly  controlled  by  the  sediment-side  (Jalalizadeh  and 
Ghosh1). 
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Figure  6.4.  Fractional  uptake  of  PCB  128  (a),  PCB  183  (b),  and  PCB  194  (c)  in  PE  passive  sampler  in 
four  differently  exposed  systems.  Experimental  data  have  been  shown  by  symbols,  and  fast  desorption 
model  simulations  have  been  shown  by  lines. 


Effect  of  Vibration  on  PRC  Loss  from  PE.  The  fraction  of  PRC  remaining  in  PE  (fpRc)  was 
measured  for  the  four  PRCs  in  static,  vibrating,  and  well  mixed  systems  (Figure  6.5).  The 
fractional  loss  of  lower  molecular  weight  PRCs  such  as  PCB  29  and  PCB  69,  were  55%  and  46% 
after  28  days  of  deployment  in  the  static  system,  respectively.  However,  both  PCBs  were  100% 
dissipated  from  PE  after  the  same  exposure  time  in  the  high  frequency  (2  min-pause)  vibrating 
system.  The  effect  of  vibration  on  PRC  loss  was  more  evident  for  larger  molecular  weight  PRCs. 
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For  example,  the  fractional  loss  was  increased  from  22%  to  87%  for  PCB  155  (hexachloro- 
congener)  and  from  12%  to  63%  for  PCB  192  (heptachloro-congener)  after  28  days.  As  shown  in 
Figure  6.5,  fpRC  for  PCB  29  and  69  were  nearly  identical  in  high  frequency  vibrating  and  well- 
mixed  exposures.  In  fact,  more  than  95%  of  both  PRCs  were  lost  in  the  vibrating  system  only 
after  7  days  of  exposure.  The  differences  between  fpRcs  in  high  frequency  vibrating  and  well 
mixed  systems  were  larger  for  more  hydrophobic  PRCs  (PCB  155  and  PCB  192).  However,  the 
fractional  losses  of  these  PRCs  were  still  reasonably  large  in  the  high  frequency  vibrating  system 
(63%  after  28  days  and  80%  after  56  days  for  PCB  192). 

The  measured  fpRc  values  in  static  system  from  experiment  1  and  2  were  identical  for  all  PRCs  as 
compared  in  the  Supporting  Information  Tables  S6.5  and  S6.6.  Thus,  only  the  results  from 
experiment  1  are  shown  for  the  static  deployment  in  Figure  6.4.  As  indicated  in  the  figure, 
vibration  with  low  frequency  (5  d-pause)  was  not  as  effective  in  improving  the  dissipation  rate  of 
high  molecular  weight  PRCs.  For  example,  the  fractional  losses  of  PCB  29  and  PCB  69  were 
increased  to  approximately  70%  after  28  days  in  the  low  frequency  vibrating  system.  However, 
the  fractional  loss  was  less  than  60%  for  155,  and  only  25%  for  PCB  192.  We  anticipate  that 
higher  fractional  loss  of  PRCs  and  closer  approach  to  equilibrium  for  the  analytes  would  result  in 
greatly  improved  accuracy  of  porewater  measurement.  This  hypothesis  has  been  tested  using 
PRC  data  in  high  and  low  frequency  vibrating  systems  in  the  following  section. 


Figure  6.5.  Comparison  of  PRC  fraction  remaining  in  PE  (fpRc)  in  static,  vibrating,  and  fully  mixed 
systems 
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Prediction  of  Porewater  Concentration  Based  on  Static  vs.  Shaken  Deployment.  PCB 
concentrations  in  PE  measured  using  the  static  and  vibration  deployment  were  corrected  for  non¬ 
equilibrium  based  on  the  fractional  loss  of  PRCs.  Since  four  PRCs  were  available  for  use  in 
correction,  the  method  by  Fernandez  et  al. 6  was  used  as  described  in  the  methods  section.  PE 
equilibrium  concentrations  were  also  determined  from  the  56-day  well-mixed  deployment.  PCB 
concentrations  in  PE  were  converted  to  Cfree  using  the  individual  congener  PE-water  partition 
constants  (Kpew).  Kpew  values  for  PCB  congeners  were  estimated  from  the  correlation  provided  in 
Ghosh  et  al.19  (Supporting  Information  Table  S6.7).  To  compare  the  accuracy  of  the  porewater 
estimations,  the  estimated  CfreeS  based  on  PRC-corrected  concentrations  in  PE  were  compared  to 
the  actual  measurements  of  Cfree  from  the  well-mixed  system.  The  PRC-corrected  CfreeS  measured 
with  28-  and  7-day  static  and  vibrating  passive  samplers  are  plotted  in  Figure  6.6  and  Figure  6.7, 
respectively  and  compared  against  the  true  CfreeS.  In  both  figures,  the  PRC-corrected  CfreeS  that 
are  statistically  different  from  the  true  values  (with  alpha  level  of  5%)  are  indicated  with  a  star. 
Since  the  in  situ  exposures  were  performed  in  lab  mesocosms  with  static  overlying  water,  we 
expect  that  the  porewater  concentrations  in  sediment  would  approach  thermodynamic 
equilibrium  with  the  sediment,  especially  for  the  strongly  hydrophobic  and  recalcitrant 
compounds.  Thus,  porewater  concentrations  measured  in  the  well-mixed  equilibrium  exposures 
should  closely  represent  the  true  porewater  concentration  to  be  expected  in  the  static  and 
vibrating  passive  sampler  deployments. 


Figure  6.6.  Comparison  of  PRC-corrected  Cfree  using  28-day  static  and  vibrating  passive  samplers  with 
measured  Cfree-  The  corrected  Cfree  was  calculated  based  on  the  PRC-corrected  PE  equilibrium  concentration 
in  static  and  vibrating  systems.  The  measured  Circe  was  estimated  based  on  the  measured  equilibrium 
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concentration  in  PE  after  56  days  of  deployment  in  the  fully  mixed  system.  The  corrected  concentrations 
that  are  statistically  different  from  the  measured  Cfree  values  (with  alpha  level  of  5%)  are  indicated  with  a 
star 
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Figure  6.7.  Comparison  of  PRC-corrected  Cfree  using  7-day  static  and  vibrating  passive  samplers  with 
measured  Cfree.  The  corrected  Cfree  was  calculated  based  on  the  PRC-corrected  PE  equilibrium  concentration  in 
static  and  vibrating  systems.  The  measured  Cfree  was  estimated  based  on  the  measured  equilibrium 
concentration  in  PE  after  56  days  of  deployment  in  the  fully  mixed  system.  The  corrected  concentrations  that 
are  statistically  different  from  the  measured  Cfree  values  (with  alpha  level  of  5%)  are  indicated  with  a  star 


As  shown  in  Figure  6.6,  the  tri-,  tetra,  and  pentachloro-congeners  are  generally  well  predicted  by 
both  static  and  vibration  passive  samplers  (within  a  factor  of  2).  For  these  moderate  Kow 
congeners  the  PRC  losses  are  high  in  either  system  and  correction  for  non-equilibrium  is 
relatively  accurate.  The  measurements  using  static  deployment  start  deviating  strongly  for  the 
dominant  higher  chlorinated  PCBs.  There  appears  to  be  a  consistent  positive  bias  for  the 
estimation  from  static  deployment.  The  error  bars  are  also  large,  indicating  poor  precision.  The 
PRC-corrected  CfreeS  from  the  low  frequency  vibration  deployment  are  more  accurate  and  have 
smaller  error  bars  compared  to  measurements  from  the  static  deployment.  However,  these 
estimations  are  statistically  different  from  the  true  values  for  highly  chlorinated  congeners,  as 
shown  in  Figure  6.6.  For  example,  for  the  dominant  heptachloro-congener,  PCB  180,  Cfree 
estimation  based  on  static  and  low  frequency  vibration  deployment  is  nearly  4-fold  and  3-fold 
higher  compared  to  the  true  Cfree,  respectively.  In  the  case  of  the  dominant  octachloro-congener, 
PCB  (203  +  193),  the  PRC-corrected  Cfree  based  on  static  and  low  frequency  vibration 
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deployment  is  nearly  5 -fold  and  4-fold  higher  than  the  true  measurement,  respectively.  In 
comparison,  the  PRC-corrected  CfreeS  from  the  high  frequency  vibration  deployment  are 
indistinguishable  from  the  true  values  for  all  PCBs. 

As  expected,  Cfree  estimations  using  the  7-day  deployment  data  in  static  system  were  less 
accurate  compared  to  the  estimations  using  28-day  deployment  data,  especially  for  higher 
molecular  weight  PCBs  (Figure  6.7).  For  example,  for  hexa-,  hepta-,  and  octachloro-congeners 
deviations  of  Cfree  from  true  values  were  within  a  factor  of  1.4  to  10  using  7-day,  and  within  a 
factor  of  1.1  to  5  when  using  28-day  data.  As  indicated  in  Figure  6.6,  PRC-corrected  Cfree  values 
estimated  with  static  deployment  are  not  statistically  different  from  the  measurements,  as  the 
estimated  concentrations  have  very  poor  precisions.  Higher  precisions  in  estimation  of  Cfree  using 
the  low  frequency  vibrating  deployment  results  in  concnetrations  that  are  statistically  different 
from  the  true  values.  This  is  while,  predicted  CfreeS  in  the  latter  system  has  less  deviations  from 
the  measured  values.  Similar  to  the  28-day  deployment  results,  the  predicted  CfreeS  for  all  PCBs 
were  statistically  indistinguishable  from  the  true  values  when  using  the  high  frequency  vibrating 
passive  sampler. 


6.4-  Implications 

Strongly  hydrophobic  compounds  such  as  large  molecular  weight  PCBs  (hexa-  to  octachloro 
biphenyls)  in  sediment  porewater  are  known  to  be  difficult  to  measure  in  situ  due  to  slow  mass 
transfer.  The  slow  mass  transfer  also  introduces  errors  when  correcting  for  non-equilibrium  using 
the  loss  kinetics  of  PRC  compounds.  In  chapter  5  we  demonstrated  the  development  of  a 
periodically  vibrating  platform  for  the  in  situ  deployment  of  passive  samplers  and  demonstrated 
that  the  technology  enhances  the  mass  transfer  of  PAHs  into  passive  samplers.  In  this  chapter, 
we  extend  the  work  by  presenting  the  applicability  of  the  vibrating  platform  in  sediment 
contaminated  with  strongly  hydrophobic  PCBs.  The  current  platform  is  more  robust  compared  to 
the  one  in  chapter  5  and  enables  deployment  of  larger  passive  samplers  to  target  low 
concentrations  in  sediment.  The  proposed  platform  allowed  accurate  measurement  of  porewater 
concentrations  of  large  molecular  weight  PCBs  in  short  deployment  times  and  reduced  the  errors 
with  the  PRC  correction  method.  Results  from  this  study  indicate  that  for  the  strongly 
hydrophobic  congeners,  the  static  deployment  is  prone  to  a  strong  positive  bias  in  estimation  of 
porewater  concentrations.  The  estimations  were  improved  even  with  a  very  low  frequency  of 
vibration  (with  5  d-pause).  However,  the  predicted  concentrations  were  still  statistically  different 
from  the  measured  values.  High  frequency  vibration  deployment  (with  2  min-pause)  in  the  other 
hand,  provided  accurate  measurement  of  porewater  concentrations  only  after  7  days  of 
deployment.  Shorter  deployment  times  lower  the  risk  of  loss,  destruction,  and  vandalism  of 
deployed  passive  sampler  platforms  in  the  field.  In  general,  the  findings  from  this  study,  pave  the 
way  for  extension  of  the  in  situ  passive  sampling  approach  to  full  range  of  PCBs,  and  possibly 
PCDD/Fs  as  well. 
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7.0  RESEARCH  OBJECTIVE  #3:  Evaluate  a  range  of  polymer  types  for 
use  as  passive  equilibrium  samplers  for  methylmercury. 


7.1  Introduction 

The  most  relevant  form  of  mercury  from  the  standpoint  of  human  health  risk  assessment  is 
monomethylmercury,  i.e.,  CH3Hg+  and  CFTHg-R,  hereafter  collectively  abbreviated  MeHg 
(Sunderland  et  al.  2007).  It  is  a  potent  neurotoxin  that  is  capable  of  crossing  biological 
membranes  and  the  blood-brain  barrier.  MeHg  is  not  typically  introduced  to  the  environment  due 
to  direct  contamination,  but  rather  as  a  result  of  the  methylation  of  inorganic  mercury(II)  (Hgi) 
by  anaerobic  microorganisms  in  sediments  and  saturated  soils  (Gilmour  et  al.  2013). 
Methylmercury  bioaccumulates  in  aquatic  food  webs  through  direct  uptake  by  single  cell 
organisms,  by  benthic  organisms  and  by  efflux  to  the  water  column  (e.g.  Williams  et  al.  2010; 
Taylor  et  al.  2014).  For  risk  assessors  and  decision  makers,  it  is  critical  to  be  able  to  predict  both 
net  MeHg  production  and  its  bioavailability  to  food  webs. 

MeHg  complexation  in  sediments  and  soils.  In  natural  waters,  both  Hgi  and  MeHg  form  a  wide 
variety  of  complexes  with  ligands  such  as  Cl',  OH ",  and  SH  .  Due  to  the  extremely  high  affinity 
of  the  soft  acid-soft  base  associations  involved,  reduced  sulfur  of  all  kinds  tends  to  dominate 
mercury  speciation  in  sediments,  sediment  pore  waters  and  many  surface  waters  (Hsu-Kim  et  al. 
2013).  This  includes  inorganic  sulfides,  simple  organic  thiols  such  as  the  amino  acids  cysteine 
and  glutathione,  and  the  reduced  sulfur  functional  groups  associated  with  dissolved  organic 
matter  (DOM)  (Skyllberg  et  al.  2006;  Haitzer  et  al.  2002).  Generally,  DOM  is  defined 
operationally  as  the  organic  fraction  that  passes  a  0.45  pm  filter  (Aiken  et  al.  201 1).  A  large  part 
of  this  fraction  comprises  humic  substances,  which  occur  as  complex  mixtures  of  aliphatic  and 
aromatic  compounds  with  varying  levels  of  solubility,  and  molecular  weights  ranging  from  a  few 
hundred  Daltons  to  many  thousands  (Liao  et  al.  1982).  The  sulfur  content  of  humic  substances  is 
typically  around  one  percent  and  comprises  a  variety  of  chemical  functionalities,  including 
sulfides,  thiols,  and  thiophenes  under  reducing  conditions  and  sulfonates  and  sulfates  under 
oxidizing  conditions  (Xia  et  al.  1998).  Conditional  stability  constants  for  DOM  complexes  with 
Hgi  and  MeHg  are  very  high  (approaching  those  reported  for  simple  thiols),  suggesting  that  thiol 
functionalities  are  the  primary  sites  of  interaction  with  DOM  molecules  (Schwarzenbach  and 
Schellenberg  1965;  Amirbahman  et  al.  2002;  Haitzer  et.  al  2002).  Because  of  the  considerable 
heterogeneity  of  reduced  sulfur  and  DOM  ligands  in  natural  waters,  Hgi  and  MeHg  species  can 
occupy  a  broad  range  of  molecular  weights,  ionic  charges,  and  hydrophibicities.  Hgi  speciation  is 
further  complicated  by  the  formation  of  non-equilibrium,  mercury-sulfide  polynuclear  clusters 
and  nanoparticles,  whose  dynamic  interactions  are  an  area  of  ongoing  research  (Deonarine  et  al 
2009;  Hsu-Kim  et  al.  2013). 
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In  most  natural  waters,  even  high  chloride  marine  systems,  high  molecular  weight  DOC 
complexes  are  the  dominant  form  of  both  Hg  and  MeHg  (e.g.  Han  et  al.  2005;  Skyllberg  et  al. 
2010).  In  anoxic  sediments  and  soils,  where  sulfide  is  commonly  present,  nanoparticulate  HgS 
may  dominate,  but  it  is  generally  associated  with  DOC  colloids  (Aiken  et  al.  201 1).  The  role  of 
sulfides  and  nanoparticles  in  MeHg  complexation  in  sulfidic  systems  has  received  little  attention. 

MeHg  bioavailability  to  aquatic  food  webs.  For  toxic  trace  metals,  bioavailability  has  often  been 
modeled  based  on  an  estimate  of  concentration  of  the  free  ion  (e.g.  AVS/SEM  model;  Allen  et  al. 
1993;  or  the  Sediment  BLM;  Di  Toro  et  al.,  2005).  However,  free  Hg2+  and  MeHg+  ions  are 
present  in  natural  waters  at  vanishingly  small  concentrations  (commonly  <  10  25  M).  Therefore,  it 
is  ineffective  and  probably  inappropriate  to  model  Hg  and  MeHg  bioavailability  based  on  their 
free  ions. 

Unfortunately,  the  MeHg  uptake  mechanism(s)  by  single  cell  organisms  at  the  bottom  of  aquatic 
food  webs  remain  effectively  unknown,  and  so  bioavailable  species  have  been  evaluated  based 
on  empirical  observations  with  selected  complexes.  An  early  paradigm  was  bioavailability  of 
neutrally  charged  species  (like  MeHgCl  and  HgCh)  via  passive  diffusion  (Mason  et  al.  1996; 
Benoit  et  al.  1999),  and  early  models  for  bioavailability  were  based  on  the  modeled  equilibrium 
concentrations  of  these  species.  Several  studies  in  the  last  decade  support  active,  or  at  least 
facilitated,  uptake  of  MeHg  by  algae  (Pickhardt  and  Fisher  2007;  Le  Faucheur  et  al.  2011,  2014) 
and  Hg(II)  uptake  by  bacteria  (Golding  et  al.  2002;  Schaefer  et  al.  2009).  One  well-known  MeHg 
active  uptake  system  is  the  broad  Hg  transporter  ( merE )  in  the  mercury  detoxification  system  of 
Hg-resistant  bacteria  (Boyd  and  Barkay  2012).  No  other  Hg  or  MeHg  transporters  have  been 
identified. 


Empirical  measurements  show  that  MeHg  salts  (MeHgOH,  MeHgCl,  MeHgNOs)  are  highly 
available  for  uptake  by  microbes  (Benoit  et  al.  1999;  Hsu-Kim  et  al.  2013).  It  is  possible  that 
algae  can  readily  access  MeHg  bound  to  small  thiols,  as  bacteria  can  do  for  Hg-cysteine 
(Schaefer  et  al.  2009;  Gilmour  et  al.  2011;  Graham  et  al.  2012).  In  oxic  natural  waters,  the 
presence  of  dissolved  organic  matter  (DOM)  inhibits  but  does  not  block  MeHg  uptake  by  algae 
(Fuengen  et  al.  2012)  and  bacteria  (Ndu  et  al.  2012)  suggesting  that  these  large  complexes  are 
not  directly  taken  up.  However,  DOC  can  dramatically  enhance  the  bioavailability  of  Hg(II)  to 
microbes  in  sulfidic  sediments  and  soils,  by  slowing  the  precipitation  of  particulate  HgS,  keeping 
Hg(II)  in  the  aqueous  phase  (Graham  et  al.  2012;  Zhang  et  al.  2012).  We  don’t  know  how 
microbes  access  Hg  from  these  high  mw  and/or  nanoparticulate  forms.  Most  likely  they  are  in 
rapid  equilibrium  with  smaller  bioavailable  forms;  alternatively  some  exchange  may  be 
occurring  on  cell  surfaces. 

Benthic  invertebrates  have  many  uptake  routes  for  MeHg,  including  ingestion  of  MeHg 
associated  with  sediment  organic  and  inorganic  matter  (e.g.  Gagnon  et  al.  1997),  grazing  from 
the  water  column,  and  probably  uptake  from  dissolved  and  particulate  MeHg-DOC  complexes 
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(Williams  et  al.  2010;  Taylor  et  al.  2014),  although  DOC  can  reduce  MeHg  bioavailability 
(Lawrence  et  al.  2001).  Laboratory  toxicity  and  exposure  tests  rarely  take  all  of  these  routes  into 
account  (Wang  and  Fisher  et  al.  1999).  Development  of  PS  that  mimics  bioavailable  MeHg 
would  allow  at  least  partial  replacement  of  expensive  exposure  studies,  and  a  wider  spatial  and 
temporal  coverage  for  lower  cost. 

To  summarize,  MeHg  may  be  available  to  aquatic  food  webs  from  colloids  and  nanoparticles  in 
addition  to  salts  and  low  molecular  weight  thiols,  but  bioavailability  from  these  classes  varies, 
and  the  mechanism(s)  of  their  bioavailability  are  not  understood.  All  of  this  implies  a  high 
degree  of  uncertainty  in  modeling  and  quantifying  MeHg  bioavailability,  and  in  developing 
passive  sampling  devices  that  assess  or  mimic  the  bioavailable  pool. 


Measuring  MeHg  bioavailability.  Pore  water  MeHg  concentrations  and/or  sediment/water  MeHg 
partition  coefficients  (Kd)  have  proved  the  most  accurate  indicators  of  MeHg  bioavailability  to 
benthic  organisms  to  date,  but  remain  uncertain  predictors  (Gilmour  et  al.  2013).  “Bioreporters” 

-  microbes  genetically  engineered  with  a  Hg  or  MeHg  transporter  (Scott  et  al.  2001;  Golding  et 
al.  2002;  Dahl  et  al.  2011;  Ndu  et  al.  2012)-  can  provide  information  of  the  concentration  of 
available  forms,  but  only  forms  available  to  those  specific  transporters.  Active  sampling  methods 
such  as  pore  water  withdrawal  via  syringe  (“sippers”)  or  sediment  core  removal  and 
centrifugation  can  be  labor  intensive  and  may  obscure  important  temporal,  spatial,  and  speciation 
information.  They  also  do  not  provide  a  measure  of  the  bioavailable  fraction  of  Hg  or  MeHg  in 
pore  water.  Equilibrium  dialysis  chambers  (“peepers”)  exclude  (or  limit)  higher  molecular 
weight  MeHg  complexes  that  may  have  some  bioavailability. 


Development  of  passive  samplers  for  Hg/MeHg.  The  two  major  strategies  for  passive  sampling 
of  contaminants  are  equilibrium  samplers  and  diffusion  based  samplers,  such  as  DGT  or 
Chemcatchers®.  Equilibrium  samplers  are  mainly  used  for  HOCs,  for  which  the  freely  dissolved 
concentration  in  sediment  pore  waters  is  widely  used  as  a  proxy  measure  of  bioavailability. 
Equilibrium  samplers  make  use  of  thin  polymers  sheets  with  partition  coefficients  for  HOCs 
similar  to  their  binding  constants  to  sediment  organic  matter.  They  effectively  concentrate 
dissolved-phase  HOCs  onto  the  polymers  sheets,  which  can  then  be  extracted  for  analysis. 
Because  of  the  relatively  small  mass  to  surface  area  ratio  of  the  polymer  exposed  to  sediment, 
the  samplers  don't  deplete  pore  water  pools,  or  create  concentration  gradients  that  could  lead  to 
significant  HOC  depletion  from  the  adjacent  sediments.  If  diffusive  boundary  layer  (DBL)  issues 
are  overcome,  and  the  samplers  come  to  equilibrium,  the  equation  for  predicting  aqueous  HOC 
concentrations  is  simple: 

Cpw  =  Csorbent/K(sorbent-water) 

Here,  C  is  the  concentration  of  the  HOC  of  interest,  and  K  is  the  partition  coefficient  between 
water  and  the  polymer  sorbent.  However,  the  bioavailable  fraction  is  MeHg  in  natural  waters  is 
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poorly  defined,  and  is  probably  not  fully  reflected  by  the  “freely  dissolved”  fraction. 

Performance  reference  materials  are  used  to  check  for  equilibrium  and  to  correct  for  non¬ 
equilibrium  when  needed. 

Non-equilibrium  diffusion  samplers  have  been  the  more  common  approach  for  metals.  These 
methods  are  based  on  the  movement  of  metal  ions  through  an  ion-permeable  diffusive  layer  and 
accumulation  in  the  binding  resin  behind.  The  diffusion  layer  limits  the  size  of  the  molecules  that 
penetrate  the  samplers,  generally  limiting  the  sampling  to  low  molecular  weight  compounds. 
Diffusion-based  samplers  for  metals  are  designed  to  estimate  the  “labile”  fraction  of  metals  of 
interest.  This  fraction  is  empirically  defined,  and  depends  on  the  sampler  design  and  chemistry. 

It  includes  free  metal  cations,  metal  salts  and  may  include  some  low  molecular  weight 
compounds  depending  on  the  diffusion  gel.  The  accumulation  rate  and  selectivity  of  the  device 
for  different  metal  complexes  are  regulated  by  the  choice  of  materials  for  the  diffusion  and 
binding  layers.  Diffusions  layers  can  be  made  with  materials  or  different  effective  pore  sizes,  and 
materials  can  be  selected  to  limit  sorption  of  the  target  compound.  The  sorptive  ligands  in  the 
binding  layer  can  also  be  adjusted  for  specific  chemicals. 

Diffusion-based  samplers  are  designed  for  use  as  non-equilibrium,  samplers.  The  binding  layer, 
containing  a  strong  metal-binding  ligand,  and  often  100’s  of  pm  thick,  acts  effectively  as  an 
infinite  sink.  The  “labile”  concentration  of  contaminant  in  the  aqueous  phase  (C dgt)  is  calculated 
based  on  a  lab-calibrated  diffusion  rate  (in  cm  sec-1)  through  the  diffusion  gel,  which  generally 
has  long  path  length  (100’s  of  pm)  in  comparison  to  the  DBL. 

The  estimated  concentration  of  “labile”  metal  (C dgt)  depends  on  both  the  size  and  the  binding 
strength  of  metal  complexes  in  solution.  If  metal-ligand  complexes  that  don’t  penetrate  the 
diffusion  gel  dissociate  rapidly  and  completely  in  solution,  the  mass  of  metal  accumulated  in  the 
DGT  binding  layer  is  the  sum  of  the  fluxes  from  each  species  present.  Effective  fluxes  from 
species  with  higher  binding  strengths  will  be  lower  (fill  eqts  from  Holger  paper).  Work  with 
several  metals  other  than  Hg  suggests  that  DGT  assesses  most  of  the  clkDa  fraction  in  natural, 
oxic  surface  waters  (Davison  and  Zhang  2012). 

By  testing  DGT  metal  accumulation  against  a  range  of  DOC  concentrations,  the  fraction  of  metal 
complexed  to  DOC  can  be  estimated.  Further,  a  comparison  of  Cog  r  with  the  molar  ratio  of  metal 
to  DOC  can  provide  information  on  DOC  binding  site  affinities  (Town  et  al.  2009). 

Additionally,  adjustment  of  the  pore  size  in  the  diffusion  gel  may  allow  some  fraction  of  metal- 
DOC  or  metal  nanoparticles  to  penetrate  the  gel  (Femandez-Gomez  et  al  2014).  Agarose  pore 
sizes  may  be  as  large  as  74  nm,  while  standard  DGT  polyacrylamide  gels  may  be  <2  nm. 

Particles  with  a  diameter  of  5-6  nm  should  have  a  diffusion  coefficient  in  water  approximately  8 
times  less  than  that  of  free  metal  ions  (Davison  and  Zhang  2012).  Thus  readily  exchangeable 
metals  on  DOC,  and  metal-DOC  complexes  may  be  accessible  by  DGT.  However,  the  samplers 
would  need  to  be  calibrated  for  each  site  and  against  a  range  of  DOC  concentrations  in  order  to 
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estimate  what  complexes  comprise  C dgt,  and  in  what  ratio. 

Chemcatchers  (Vrana  et  al.  2005)  also  have  diffusion  and  binding  layers,  made  from 
commercially-available  thin  membranes.  The  diffusive  boundary  layer  in  the  aqueous  phase 
becomes  particularly  important  for  these  thin  diffusion  layers.  For  Chemcatcher-type  samplers, 
Cw  is  estimated  from  an  empirical  “sampling  rate,”  Rs,  with  units  of  volume  sampled  over  time. 
Rs  is  also  derived  from  laboratory  calibration: 

m  =  C  Rst 

Where  m  equals  the  mass  of  analyte  accumulated  in  the  receiving  phase  after  an  exposure  time 
(t)  and  C  is  the  analyte  concentration  in  solution  and  Rs  is  sampling  rate  in  ml/h 

Of  these  approaches,  DGT  has  been  the  most  fully  developed  for  inorganic  Hgll  (hereafter  Hgll) 
and  MeHg  (Docekalova  et  al.  2005;  Divis  et  al.  2005;  Clarisse  et  al.  2009,  2011;  Divis  et  al. 

2010;  Gao  et  al.  2011;  Fernandez-Gomez  et  al.  2011,  2014,  2015;  Pelcova  et  al.  2015;  Noh  et  al. 
2015).  Several  adjustments  to  standard  metals  DGT  have  been  made.  Ideally  a  diffusion  gel 
should  have  diffusion  coefficient  the  same  as  water,  and  ions  should  move  freely  through  the  gel 
to  the  binding  layer.  However,  Docekalov  et  al.  (2005)  measured  a  concentration  factor  of  700 
for  Hg(II)  in  the  polyacrylamine  gel  commonly  used  in  DGT,  relative  to  water  (Hg  binds  amide 
groups  in  the  polyacrylamide).  They  switched  to  agarose  diffusion  gel,  which  showed  only  a 
factor  of  4  concentrations  from  water.  Membrane  properties  have  also  been  adjusted  to  also  limit 
sorption  (e.g.  Clarisse  et  al.  2006;  Fernandez-Gomez  et  al.  2014).  Several  binding  materials  have 
been  tested  to  improve  sorption  and  specificity  (Docekalov_et  al.  2005;  Clarisse  et  al.  2006; 

Divis  et  al.  2010;  Gao  et  al.  2011;  Fernandez-Gomez  et  al.  2014).  All  of  the  most  successful  are 
some  form  of  thiol,  most  effectively  3-mercaptopropyl-functionalized  silica  gel  (Gao  et  al.  2011; 
Clarisse  et  al.  2009,  2011). 

Some  deployments  of  MeHg  DGT  samplers  have  predicted  aqueous  total  MeHg  concentrations 
reasonably  well,  especially  in  saline  surface  waters  where  MeHgCl  may  be  an  important 
complex  (Clarisse  et  al.  2006;  2011;  Noh  et  al.  2015).  However,  the  effective  diffusion  rate  for 
MeHg  must  be  tested  ahead  of  time  in  the  lab  for  water  evaluated  (Noh  et  al.  2015).  In  high  DOC 
waters,  DGT  samplers  under  predicts  Cw,  unless  one  assumes  different  diffusivities  for 
MeHgDOC,  as  calibrates  the  samplers  to  the  sample  water.  In  one  of  only  tests  of  DGT  samplers 
for  prediction  of  MeHg  uptake  by  animals  (Clarisse  et  al.  2009),  DGT  samplers  predicted  MeHg 
uptake  by  Macoma  balthica  under  controlled  conditions  in  the  lab,  across  a  wide  range  of 
salinities;  but  under  predicted  uptake  at  very  low  MeHg  concentrations,  or  when  DOC  was  added 
to  the  tests.  DGT  samplers  specific  for  MeHg  also  predicted  MeHg  uptake  by  rice  after  the 
samplers  were  calibrated  for  the  rice  paddy  soil  pore  waters  (Liu  et  al.  2012).  In  sulfidic  systems, 
DGT  models  do  particularly  poorly  in  predicting  Cw  (Clarisse  et.  al  2009).  Laboratory  studies 
show  the  presence  of  HgS  nanoparticles  in  sulfidic  samples  deposit  on  the  surface  of  samplers, 
slowing  the  migration  of  dissolved  species  into  the  samplers  (Pham  et  al.  2015). 
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To  summarize,  the  most  developed  passive  sampling  approach  for  MeHg,  DGT,  remains 
complex  to  use,  but  can  provide  estimates  of  bioavailable  MeHg  if  the  samplers  are  calibrated 
based  on  a  specific  measured  diffusion  rates  for  the  study  water.  They  provide  uncertain  results 
in  with  regard  to  the  precise  species  of  mercury  sampled,  although  a  range  of  deployment  times 
and  gel  thicknesses  can  be  used  to  help  assess  MeHg  size,  speciation  and  its  relationship  to 
bioavailability.  DGT  samplers  require  site-specific  calibration  (Pelcova  et  al.  2014;  Noh  et  al. 
2015).  DGT  does  not  directly  access  higher  molecular  weight  complexes  like  MeHg-DOC 
colloids  HgS  nanoparticles  that  are  potentially  bioavailable  forms  of  MeHg. 


MeHg  passive  sampler  development  strategy.  There  remains  a  need  for  a  simple,  passive 
sampling  system  designed  specifically  to  quantify  the  time-weighted  average  bioavailable  MeHg 
concentration  in  a  reproducible  fashion.  However,  we  still  must  define,  at  least  empirically,  the 
MeHg  complexes  and  phases  that  are  available  to  organisms.  We  believe  that  one  of  the  most 
important  aspects  of  MeHg  passive  sampler  development  will  be  the  use  of  a  range  of  sampler 
types  and  chemistries  to  begin  to  understand  and  define  the  bioavailable  pool  of  MeHg  in 
sediments  and  soils. 


While  our  main  goal  is  to  develop  a  passive  sampler  that  provides  an  estimate  of  the  bioavailable 
fraction  of  MeHg,  our  near-term  objective  is  to  use  a  range  of  passive  sampler  chemistries  and 
designs  to  assess  the  chemical  and  physical  forms  of  MeHg  available  to  organisms.  This  would 
be  an  improvement  on  current  methods  for  MeHg,  which  assume  the  “freely  dissolved”  fraction 
is  equivalent  to  the  bioavailable  pool. 

Our  general  approach  is  to  develop  an  equilibrium  method  in  which  the  binding  material  mimics 
the  chemistry  and  partition  coefficient  of  natural  MeHg  ligands  on  cell  surfaces.  We  chose  an 
equilibrium  method  to  try  and  avoid  the  limitations  of  methods  that  rely  on  kinetic  calibrations 
for  each  potential  MeHg  complex.  We  are  testing  both  diffusion-based  and  single-layer  polymers 
gels  in  order  to  evaluate  the  bioavailability  MeHg  in  colloids  or  nano-particles. 

In  this  study,  we  began  development  of  MeHg  samplers  by  evaluating  potential  materials.  Our 
goal  in  this  research  was  to  identify  and  prepare  materials  with  affinity  for  MeHg  and  capable  of 
accumulating  detectable  and  reproducible  amounts  in  proportion  to  aqueous  concentrations. 

More  specifically,  we  are  working  toward  a  small  set  of  materials  for  testing  against  animal 
bioavailability.  These  materials  should  have  a  range  of  properties  including: 

•  Partition  coefficients  close  to  the  measured  K  for  MeHg  in  natural  sediments  (10A3  to 
10A5). 

•  Linear  uptake  isotherms  across  the  concentration  range  of  interest,  i.e.  from  the  pg/L 
range  commonly  found  in  natural  uncontaminated  sediment  pore  waters  to  the  high  ng/L 
concentrations  found  in  the  most  contaminated  pore  waters. 
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•  Capable  of  accumulating  detectable  and  reproducible  masses  of  MeHg  in  thin  layer 
deployments  (ultimately  10’s  of  microns)  at  natural  MeHg  concentrations. 

o  The  absolute  MeHg  DLs  in  our  lab  are  about  1-5  pg  per  sample.  A  good  target  for 
easy  detection  for  labs  with  less  sophisticated  analytical  methods  (i.e.  CVAf) 
would  be  ~0.1  ng  per  sample. 

•  Reasonably  rapid  equilibrium  from  complex  aquatic  matrices,  including  dissolved 
organic  matter. 

•  Ability  to  manufacture  in  consistent  thin  films,  to  limit  depletion  of  MeHg  complexes  in 
pore  waters 


A  range  of  polymer  and  polymer  composites  were  prepared  to  test  sorption  capacity  and 
linearity.  Target  sorbents  included  activated  carbon  and  a  suite  of  thiol  ligands  and  thiol- 
functionalized  sorbents  in  the  solid  and  aqueous  phase.  Polymers  were  chosen  for  pore  sizes, 
sorbent  capacity  of  MeHg,  and  ability  to  potentially  include  sorbents.  A  wide  variety  of  polymers 
with  and  without  sorbent  inclusions  were  synthesized  and  tested.  The  binding  capacity  of 
polymers  and  polymer/sorbents  mixtures  were  measured.  MeHg  partitioning  coefficients  were 
measured  in  14-d  isotherms  in  3  ppt  saline  solutions  using  simple  MeHg  salts.  Several  sets  of 
isotherms  were  done  as  we  learned  about  the  properties  of  our  materials  and  began  to  build 
mixtures  that  best  approximated  MeHg  equilibrium  with  sediments.  Finally,  materials  showing 
promise  in  isotherms  were  tested  under  more  environmentally  realistic  experimental  conditions 
with  the  inclusion  of  DOM. 


7.2  Materials  and  methods 

Development  of  passive  sampling  polymers  for  methylmercury.  Many  polymers  with  and 
without  sorbent  inclusions  were  tested  in  this  study  (Table  7.1).  The  first  materials  tested  were 
simple  hydrophobic  polymers.  Some,  including  polyethylene  (PE),  polyoxymethylene  (POM), 
and  polydimethylsiloxane  (PDMS),  are  widely  used  for  passive  sampling  of  organic 
contaminants.  However,  these  standard  polymers  were  not  sufficiently  sorptive  for  MeHg,  which 
is  much  less  lipophilic  (see  Results  and  Discussion).  To  address  this,  attempts  were  made  to 
achieve  greater  chemical  affinity  by  modifying  existing  polymer  types.  To  this  end,  two  separate 
strategies  were  investigated  in  parallel.  The  first  was  incorporation  of  activated  carbon  in  the 
sampling  material.  Selection  of  ACs  for  trial  was  based  on  Gomez-Eyles  et  al.  (2013),  in  study  of 
Hgi  and  MeHg  isotherms  for  a  variety  of  commercially  available  and  lab-synthesized  ACs  and 
biochars.  All  exhibited  linear  partitioning  across  a  range  of  environmentally  relevant  water 
concentrations,  with  log  Kd  values  from  3.8  to  5.5.  It  may  be  possible  to  exploit  this  behavior  for 
a  passive  sampling  application  by  immobilizing  AC  in  a  polymeric  matrix.  Choi  (2010)  reported 
a  method  for  embedding  powdered  activated  carbon  in  a  poly(vinylidene  fluoride)  (PVDF) 
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polymer.  This  procedure  was  followed  in  the  present  work  and  also  adapted  to  produce  an 
assortment  of  PVDF-ACs  with  varying  carbon  contents  and  physical  attributes  (Figure  7.1).  For 
example,  it  is  possible  to  adjust  the  cross-sectional  pore  structure,  porosity,  and  permeability  of 
PVDF  by  manipulating  variables  such  as  solvent  proportions  during  synthesis  (Sukitpaneenit  and 
Chung  2009;  Wang  et  al.  2008;  Wang  et  al.  2013).  Agarose  was  chosen  as  another  possible 
support  material  for  AC  because  it  has  proven  effective  as  the  diffusive  gel  in  DGT  sampling  for 
MeHg  (Gao  et  al.  2014).  It  is  also  possible  to  manipulate  the  structure  of  agarose  gels  by 
incorporating  dopants  during  gel  preparation  to  promote  either  micro-  or  macroporosity 
(Charlionet  et  al.  1996).  For  this  project,  agarose  gels  containing  suspended  AC  were  prepared 
and  found  to  be  stable  in  water  for  at  least  28  d,  the  typical  deployment  time  for  our  lab’s  passive 
sampler  deployments  in  the  field. 


Figure  7.1.  Assortment  of  agarose  polymers.  Clockwise  from  upper  left:  with  embedded  GAC; 
with  embedded  GAC  and  coated  on  a  glass  fiber  filter;  with  embedded  thiol-SAMMS  and  coated 
on  a  glass  fiber  filter;  with  embedded  thiol-SAMMS. 


The  second  type  of  polymer  modification  investigated  was  the  incorporation  of  reduced  sulfur 
chemical  functionality  in  the  sampling  material.  This  approach  not  only  promises  to  enhance 
affinity  for  MeHg,  but  may  also  represent  a  step  toward  the  goal  of  biomimesis  by  drawing  upon 
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the  chemistry  underlying  many  of  mercury’s  biological  interactions.  Biomedical  researchers 
have  developed  an  assortment  of  thiolated  polymers  for  use  in  pharmaceuticals,  where  their 
ability  to  form  sulfide  bonds  confers  mucoadhesive  properties  that  facilitate  controlled  drug 
delivery  (Bernkop-Schniirch  2005).  Several  of  these  were  synthesized  for  the  present  work, 
including  a  chitosan-thioglycolic  acid  conjugate  (Kast  and  Bernkop-Schnurch  2001),  and 
alginate  and  xyloglucans  functionalized  with  L-cysteine  (Bernkop-Schnurch  et  al.  2001; 
Bhalekar  et  al.  2013).  Generally,  the  physical  and  mechanical  properties  of  these  materials  as 
prepared  were  unsuitable  for  deployment  in  sediment,  so  suspensions  of  the  materials  in  agarose 
were  prepared  to  obtain  greater  stability.  Also  included  in  preliminary  isotherms  were 
diatomaceous  earth  particles  that  were  covalently  thiolated  according  to  Yu  et  al.  (2012)  and 
suspended  in  agarose.  Other  tested  materials  included  polyethylene  terephthalate  sheets  with 
covalently  immobilized  L-cysteine,  PVDF  doped  with  L-cysteine,  and  PVDF  and  agarose 
suspensions  of  thiol-SAMMS  (Self-Assembled  Monolayers  on  Mesoporous  Supports,  an 
engineered  sorbent  for  aqueous  metal  cations  (Feng  et  al.  1997)).  Photographs  of  some  of  the 
PVDF-based  polymer  composites  are  shown  in  Figure  7.2. 


Figure  7.2.  Assortment  of  magnified  photos  of  lab-synthesized  PVDF  polymers.  Clockwise  from 
upper  left:  unmodified;  with  embedded  cysteine;  with  embedded  thiol-SAMMS  and  coated  on  a 
glass  fiber  filter;  with  embedded  GAC  and  coated  on  a  glass  fiber  filter. 
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MeHg  isotherm  studies  with  polymers.  In  initial  isotherm  experiments,  a  piece  of  test  polymer 
was  added  to  50  mL  deionized  water  adjusted  to  3  ppt  salinity  with  Instant  Ocean  (Spectrum 
Brands,  Blacksburg,  VA)  and  spiked  with  MeHg  at  one  of  a  series  of  concentrations  typically 
spanning  two  orders  of  magnitude.  Samples  were  contained  in  60-mL  polyethylene  terephthalate 
glycol  copolymer  (PETG)  bottles.  Incubations  were  carried  out  at  4  °C  with  orbital  shaking  at 
120  rpm.  Following  incubation,  samples  were  passed  through  0.45  pm  glass  microfiber  filters 
using  disposable  plastic  syringes.  Polymers  were  saved  and  frozen,  and  filtered  water  samples 
were  distilled  and  analyzed  according  to  methods  described  by  Mitchell  and  Gilmour  (2008). 
Briefly,  an  isotopic  dilution  spike  was  added,  the  samples  were  buffered  with  citrate,  derivatized 
with  sodium  tetraethylborate  to  facilitate  volatilization,  purged  and  concentrated  on  a  Tenax  trap, 
thermally  desorbed,  separated  on  an  OV-3/Chromasorb  column,  and  introduced  into  an  ICP-MS 
for  detection.  Kpoiymer  for  each  sample  was  calculated  as  follows: 

Kpolymer  —  ~  (1) 

where  Cw  is  the  measured  final  concentration  in  water,  and  Cs  is  the  sorbent  concentration.  The 
mass  of  MeHg  on  each  sorbent  was  inferred  from  the  difference  between  spiked  and  “recovered” 
MeHg  (the  final  measured  concentration  multiplied  by  the  sample  volume).  A  significant  amount 
of  sorption  to  PETG  bottle  walls  occurred  in  blank  spikes,  and  this  could  have  introduced 
substantial  error  in  the  Kd  measurements  of  very  weak  sorbents.  However,  this  phenomenon  is  a 
much  smaller  source  of  relative  error  for  sorbents  stronger  than  PETG  by  one  or  more  orders  of 
magnitude  (i.e.,  any  polymers  worthy  of  further  consideration).  For  this  reason,  no  blank 
correction  was  applied  to  calculations. 

For  this  preliminary  screening  phase,  a  target  range  of  sampler  partitioning  was  set  at  between 
103  °  and  104  5  L  kg"1.  This  range  was  chosen  as  a  first  attempt  to  approximate  benthic 
bioaccumulation  and  so  that  the  sampler  can  compete  with  sediment  for  MeHg.  For  comparison, 
MeHg  log  Kd  values  in  the  Chesapeake  Bay  have  been  measured  in  the  range  of  about  2.5  to  4 
(Hollweg  et  al.  2009).  It  is  also  important  to  avoid  oversampling.  Materials  with  excessively  high 
affinity  can  act  as  sinks  for  MeHg  in  pore  water,  inducing  desorption  and  resupply  from 
sediment  stores  and  resulting  in  overestimation  of  equilibrium  water  concentrations  (Claris  se  et 
al.  2009). 


Mercury  and  Methylmercury  Analyses.  Total  Hg  and  MeHg  analyses  were  carried  out  by  isotope 
dilution  ICP-MS,  after  sample  preparation  by  digestion,  distillation,  extraction  or  other  methods 
as  described  in  detail  in  our  cited  publications  (Gilmour  et  al.  2013;  Graham  et  al.  2012;  Mitchell 
et  al.  2012;  Hollweg  et  al.  2010).  These  methods  are  adaptations  of  EPA  Methods  1630  and 
1631.  Isotope  dilution  ICP-MS  significantly  improves  the  accuracy  and  precision  of  these 
multistep  methods.  QA/QC  will  include  blanks,  replicates,  spikes  and  suitable  certified  reference 
materials. 
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Modeling  sorption  in  polymers.  To  assess  diffusion  of  MeHg  into  the  test  polymers,  measured  K 
values  were  compared  with  those  that  would  be  predicted  strictly  on  a  mass  basis.  For  example, 
given  the  composition  of  an  AC-containing  polymer  and  the  Kd  values  of  the  AC  measured  by 
Gomez-Eyles  et  al.  (2013;  105  4  L  kg  *)  and  the  unmodified  polymers  (this  work),  one  could 
make  the  following  approximation: 


K 


=  f, 


d, predicted  1polymerKpolymer  T  lAC^AC 


+  fArK/ 


(2) 


where  f  is  the  fractional  contribution  of  each  constituent  to  the  total  mass  of  the  material. 


Most  of  the  modified  polymers  performed  well  (log  K  in  the  target  range).  To  narrow  the 
selection,  a  subset  of  the  materials  with  the  most  desirable  partitioning  behavior  and  physical 
properties  was  chosen  to  proceed  to  a  second  phase  of  testing.  This  phase  was  designed  to 
investigate  sampler  partitioning  of  MeHg  in  the  presence  of  dissolved  organic  matter  (DOM),  a 
crucial  environmental  ligand.  Here,  MeHg  isotherms  were  set  up  similarly,  but  with  the  addition 
of  Suwannee  River  Humic  Acid  II  (hereafter  SRHA;  International  Humic  Substances  Society),  a 
well-characterized,  organic -rich  isolate  with  a  high  degree  of  aromaticity  that  tends  to  make  it 
reactive  toward  mercury  (Graham  et  al.  2013).  SRHA  was  added  to  achieve  a  DOM: MeHg  ratio 
of  106  in  each  sample.  This  was  done  because,  according  to  Haitzer  et  al.  (2002),  affinity 
experiments  at  DOM:Hg  mass  ratios  at  or  below  105  are  misleading  because  they  give  the 
impression  that  mercury  associates  primarily  with  the  relatively  plentiful  hydroxyl  groups  on 
DOM.  In  fact,  speciation  is  controlled  by  a  small  fraction  of  reactive  thiol  functional  groups 
under  typical  environmental  conditions.  MeHg  spikes  were  chosen  in  consideration  of  the 
detection  range  of  the  analytical  method  and  the  anticipated  partitioning  coefficients  of  the 
materials  in  the  presence  of  DOM,  which  as  a  first  guess  were  anticipated  to  be  one  order  of 
magnitude  lower  than  those  measured  in  the  DOM-free  solutions.  DOM  and  MeHg  were  allowed 
to  equilibrate  overnight  before  the  addition  of  sorbents  (Figure  7.3).  The  samples  incubated  for 
14  d  and  were  processed  and  analyzed  for  MeHg  as  before. 


Figure  7.3.  Sample  bottles  for  MeHg-DOM/polymer  isotherm  experiment. 
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7.3  Results  and  discussion 

The  first  polymers  tested,  including  PE,  POM,  and  PDMS,  have  sufficiently  high  affinity  for 
PCBs  and  polyaromatic  hydrocarbons  (PAHs)  to  serve  as  passive  sampling  materials  due  to  the 
hydrophobicity  of  the  compounds  (log  Kow  between  3  and  8  (Hawker  and  Connell  1988;  de 
Maagd  et  al.  1998)).  MeHg  salts,  by  contrast,  are  only  weakly  hydrophobic,  if  at  all  (log  Kow 
between  -1  and  2  (Mason  2002)).  It  is  therefore  not  surprising  that  the  measured  partitioning 
coefficients  in  this  first  experiment  were  one  to  two  orders  of  magnitude  below  our  target  range 
(see  Figure  7.4  below).  This  experiment  also  included  early  attempts  at  AC-enhanced  PVDF  and 
agarose.  The  initial  agarose+AC  performed  poorly,  leading  to  adjustments  in  the  preparation 
procedure  for  subsequent  experiments.  However,  the  measured  Kd  of  the  PVDF+AC  was  roughly 
an  order  of  magnitude  higher  than  those  of  the  unmodified  polymers,  providing  the  first  evidence 
that  the  modification  strategy  could  be  effective.  Overall,  the  key  findings  of  the  first  set  of 
isotherm  studies  were  that  the  pure  polymers  had  weak  sorption  for  MeHg  and  inclusion  of  AC 
particles  greatly  enhanced  sorption  and  maintained  the  observed  linear  sorption  behavior  of 
MeHg  in  AC. 
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Figure  7.4.  Results  of  first  MeHg  isotherm  study  focusing  on  unmodified  hydrophobic 
polymers.  See  Table  7.1  for  descriptions  of  all  materials  tested. 
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A  second  set  of  isotherms  provided  further  support  for  the  polymer-AC  approach  (see  Figure  7.5 
below).  An  assortment  of  PVDF+ACs  prepared  using  different  ACs  and  polymerizing 
nonsolvents  performed  similarly,  and  all  were  an  order  of  magnitude  more  sorptive  than 
unmodified  PVDF.  A  second  iteration  of  agarose+AC  partitioned  much  more  strongly  than  the 
first,  and  more  than  unmodified  agarose.  The  thiolated  polymers  tested  in  this  set  leached  sulfur 
into  solution,  causing  analytical  interferences  and  confounding  the  results.  For  subsequent 
preparations,  greater  care  was  taken  to  pre-clean  these  materials,  leading  to  fewer  analytical 
difficulties. 
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Figure  7.5.  Results  of  second  isotherm  study,  including  AC-  and  thiol-based  polymers.  See 
Table  7.1  for  descriptions  of  all  materials  tested. 


The  third  set  of  isotherms  focused  on  a  selection  of  thiolated  materials  (Figure  7.6).  PVDF  with 
suspended  cysteine  performed  very  well,  while  an  agarose  suspension  of  cysteine  did  not.  This 
may  have  been  because  the  cysteine  was  not  stable  in  the  more  hydrophilic  and  porous  agarose 
matrix  and  may  have  partially  leached  into  solution.  By  contrast,  the  agarose  suspension  of  thiol- 
SAMMS  was  the  strongest  sorbent  yet  tested,  and  the  only  one  to  date  with  a  measured  log  Kd 
greater  than  5.  The  biomedical  thiomers  performed  reasonably  well,  but  given  the  complexities 


58 


involved  in  their  preparation  and  the  fact  that  they  were  no  more  sorptive  than  polymers 
containing  SAMMS  or  AC,  they  were  excluded  from  further  consideration.  All  polymers  tested 
with  inclusions  of  thiolated  ligands  demonstrated  strong  and  linear  sorption  behavior  as 
demonstrated  in  Figure  7.6. 


Figure  7.6.  Results  of  third  isotherm  study,  including  thiomers  and  other  sulfur-based  polymers. 
See  Table  7.1  for  descriptions  of  all  materials  tested. 


In  the  fourth  set  of  screening  isotherms  (Figure  7.8),  the  mercapto-functionalized  diatomaceous 
earth  proved  to  be  an  exceedingly  strong  sorbent  and,  importantly,  was  as  effective  as  predicted 
when  suspended  in  agarose  (Kmeasured:KPredicted  =  1).  A  second  preparation  of  cysteine- 
functionalized  polyethylene  terephthalate  (Figure  7.7)  also  performed  well.  Suspensions  of  AC 
and  SAMMS  in  an  agarose  doped  with  polyethylene  glycol  (which  promotes  macroporosity) 
performed  similarly  to  their  undoped  analogs.  This  result  is  not  surprising  in  the  simple, 
moderately  saline  test  solutions,  but  also  leaves  open  the  possibility  that  PEG  doping  could 
confer  different  sorption  behavior  in  more  realistic  matrices  with  bulkier  ligands. 
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Figure  7.7.  Polyethylene  terephthalate  covalently  modified  with  cysteine. 


Figure  7.8.  Results  of  fourth  isotherm  study.  See  Table  7.1  for  descriptions  of  all  materials 
tested. 
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For  completeness,  thiol-SAMMS  and  unmodified  PET  were  tested  in  a  final  screening  isotherm. 
The  results  were  used  in  calculations  of  Kpredicted  and  are  shown  in  Figure  7.9.  As  expected, 
partitioning  by  thiol-SAMMS  was  very  strong,  and  comparable  to  that  measured  for  DE- 
MPTMS,  another  thiolated,  highly  porous  silica  material.  Refer  to  the  attached  table  for  a 
complete  summary  of  screening  isotherm  results. 
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Figure  7.9.  Results  of  fifth  and  final  screening  isotherm.  See  Table  7.1  for  descriptions  of  all 
materials  tested. 


The  most  promising  polymers  in  the  screening  phase  were  included  in  MeHg/DOM  isotherms. 
Given  the  composition  of  the  sample  solutions  in  this  experiment,  the  chemical  speciation 
software  MINEQL+  (v.  4.6;  Environmental  Research  Software)  predicted  virtually  complete 
complexation  of  MeHg  by  DOM  prior  to  addition  of  the  sorbent  materials.  As  such,  the 
experiment  can  be  said  to  have  measured  the  ability  of  the  materials  to  partition  these  complexes. 


61 


All  of  the  tested  materials  exhibited  linear  partitioning  across  the  concentration  range  (Figure 
7.10).  However,  partitioning  was  depressed  by  an  average  of  78%  compared  to  the  screening 
isotherms,  in  which  MeHgCl  was  the  dominant  species  in  solution  according  to  MINEQL+.  The 
materials  performed  remarkably  similarly,  with  an  average  log  Kd  of  3.1.  As  shown  in  Figure 
7.10  below,  for  each  of  the  sorption  isotherms,  the  correlation  coefficient  (in  log-log  scale)  was 
high  with  most  values  at  0.99  with  concentration  values  in  water  spanning  nearly  3  orders  of 
magnitude.  It  is  important  to  note  that  the  aqueous  concentrations  were  in  the  typical  range  of 
environmental  values  of  a  few  to  hundreds  of  ng/L. 
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Figure  7.10.  Results  of  MeHg-DOM  isotherms.  See  Table  7.1  for  descriptions  of  all  materials 
tested. 
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7.4  Implications 

In  this  research  we  demonstrate  for  the  first  time  that  a  composite  polymer  can  be  synthesized 
that  responds  linearly  to  aqueous  concentration  of  MeHg  in  water  at  equilibrium,  in  the  presence 
of  complexing  ligands  chloride  and  DOM.  Thus,  conceptually,  we  can  develop  an  equilibrium 
passive  sampler  for  MeHg  in  sediment  porewater  that  is  able  to  respond  to  porewater 
concentration  of  MeHg. 

A  material  with  a  log  Kd  of  3  with  respect  to  MeHg,  exposed  to  an  aqueous  MeHg  concentration 
of  5  ng/L  (a  moderate-to-low  value  for  contaminated  sediment  pore  water),  would  be  capable  of 
accumulating  500  pg  of  MeHg  (a  readily  detectable  amount)  in  a  100  mg  mass  of  polymer.  For 
most  of  the  materials  under  consideration,  this  is  roughly  the  size  of  a  penny.  Several  of  these 
materials  achieved  this  degree  of  partitioning  of  MeHg  complexed  to  DOM,  one  of  the  most 
likely  ligands  to  dominate  speciation  in  sediments  and  also  one  of  the  least  diffusive.  This 
constitutes  a  promising  proof  of  concept  for  an  equilibrium-based,  polymeric  passive  sampling 
system  for  MeHg  in  sediment  pore  water. 

Before  that  goal  is  reached,  a  few  important  questions  remain  to  be  answered.  In  many  cases, 
Kmeasured: Kpredicted  was  significantly  less  than  1.  This  could  be  because  MeHg  did  not  have 
sufficient  time  to  diffuse  throughout  the  polymers  in  14  d  (i.e.,  a  kinetic  limitation),  or  because 
sorbent  particles  suspended  in  the  polymers  were  permanently  inaccessible  to  MeHg  due  to 
limitations  in  the  effective  porosity  of  the  materials  as  prepared  (implying  that  surface  adsorption 
was  the  primary  mechanism  of  accumulation).  To  address  these  questions,  the  kinetics  of  MeHg 
accumulation  by  polymers  should  be  studied  in  a  time  series  experiment.  It  may  also  be  possible 
to  devise  an  in  situ  calibration  method  for  polymer  equilibration,  perhaps  by  pre-loading  the 
sampler  with  isotopically-labelled  MeHg  as  a  performance  reference  compound.  Several  other 
practical  aspects  of  deployment  in  an  actual  sediment  will  require  investigation  and  refinement, 
including  a  method  for  eluting  MeHg  from  the  sampler  and  the  physical  configuration  of  the 
final  device.  Finally,  the  utility  of  the  sampler  should  be  validated  by  comparing  its 
measurements  to  bioaccumulation  by  a  relevant  test  organism  in  a  side-by-side  deployment. 
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Table  7.1.  Summary  of  results  of  MeHg  partitioning  screening  isotherms. 


Material 

Abbreviation 

Reference 

l°gKd 

(Kd  in  L/kg) 

log  Cs/log 
Cw  slope 

log  Cs/log 
Cw  int. 

log  Cs/log 

cw  R2 

f  polymer 

*i sorbent 

log  Kpredict 

^measured 

^predicted 

polyethylene 

PE 

n/a 

2.35 

0.78 

2.9 

0.92 

polytetrafluoroethylene  (Teflon) 

PTFE 

n/a 

2.19 

0.67 

2.6 

0.97 

polydimethylsiloxane 

PDMS 

n/a 

2.14 

0.57 

2.6 

0.99 

polyoxymethylene,  38  pm  thick 

POM38 

n/a 

2.05 

0.63 

3.0 

0.99 

cellulose  acetate 

CA 

n/a 

2.58 

0.54 

2.6 

0.93 

cellulose  nitrate 

CN 

n/a 

2.95 

0.69 

3.2 

0.91 

polyethersulfone 

PES 

n/a 

2.72 

0.76 

2.8 

0.72 

heterogeneous  mixture  of  alkanes  (paraffin) 

Parafilm 

n/a 

2.06 

1.1 

1.9 

0.96 

PAC  suspended  in  PVDF;  preparation  was  not  optimized 

PVDF+PAC  [old] 

Choi  2010 

3.43 

0.75 

3.8 

0.99 

0.89 

0.11 

4.5 

8E-02 

PAC  suspended  in  agarose;  preparation  was  not  optimized 

ag+PAC  [old] 

Gao,  et  al.  2011 

2.60 

0.55 

3.4 

0.96 

0.95 

0.05 

4.2 

3E-02 

powdered  activated  carbon 

PAC 

n/a 

5.48 

0.84 

5.6 

0.97 

agarose 

ag 

Gao,  et  al.  2011 

2.59 

0.53 

3.2 

0.96 

powdered  activated  carbon  suspended  in  agarose 

ag+PAC 

Gao,  et  al.  2011 

3.45 

0.86 

3.7 

1.0 

0.96 

0.04 

4.1 

2E-01 

poly(vinylidene fluoride)  prepared  with  water  nonsolvent 

PVDF  Dl 

Choi  2010 

2.29 

0.12 

3.3 

0.13 

coal-based  granular  AC  suspended  in  poly(vinylidene  fluoride) 

PVDF+GAC 

Choi  2010 

3.39 

1.0 

3.3 

0.84 

0.89 

0.11 

3.9 

3E-01 

PAC  suspended  in  PVDF  prepared  with  water  nonsolvent 

PVDF+PAC  Dl 

Choi  2010 

3.24 

0.92 

3.4 

0.98 

0.89 

0.11 

4.5 

5E-02 

PAC  suspended  in  PVDF  prepared  with  1:1  methanol :water  nonsolvent 

PVDF+PAC  meOH 

Choi  2010 

3.58 

0.89 

3.7 

0.98 

0.89 

0.11 

4.5 

IE-01 

mercapto-functionalized  cellulose  acetate 

CA+MA 

Aoki,  et  al.  2007 

2.36 

1.0 

2.2 

0.99 

L-cystei ne-functi ona  1  i zed  pol yuretha ne 

PU+cys 

Duan  and  Lewis  2002 

1.92 

0.56 

2.9 

0.46 

L-cystei  ne-functi  ona  1  i  zed  pol  yethyl  ene  terephtha  late 

PET+cys 

Duan  and  Lewis  2002 

1.11 

n/a 

n/a 

n/a 

L-cystei  ne-functi  ona  li  zed  alginate  suspended  in  agarose 

alg+cys  ag 

Bernkop-Schnurch,  et  al.  2001 

3.59 

0.72 

4.1 

0.94 

mercapto-functionalized  chitosan  suspended  in  agarose 

chit+MAA  ag 

Kastand  Bernkop-Schnurch  2001 

3.21 

0.78 

3.8 

0.98 

L-cystei  ne-functi  ona  li  zed  xyloglucans  suspended  in  agarose 

xylo+cys  ag 

Bhalekar,  et  al.  2013 

3.45 

0.73 

4.0 

0.95 

L-cysteine  dissolved  in  agarose 

ag+cys 

Gao,  et  al.  2011 

2.94 

1.0 

2.8 

0.98 

0.98 

0.02 

2.6 

2E+00 

L-cystei  ne  dissolved  in  PVDF 

PVDF+cys 

Choi  2010 

4.36 

0.75 

4.8 

0.99 

0.94 

0.06 

2.4 

1E+02 

thiol-SAMMS  suspended  in  agarose 

ag+SAMMS 

Gao,  et  al.  2011 

5.11 

0.54 

5.2 

0.84 

0.98 

0.02 

3.8 

2E+01 

thiol-SAMMS  suspended  in  PVDF 

PVDF+SAMMS 

Choi  2010 

4.11 

0.82 

4.4 

0.98 

0.947 

0.053 

4.2 

7E-01 

diatomaceous  earth 

DE 

n/a 

2.96 

0.66 

3.2 

0.93 

mercapto-functionalized  diatomaceous  earth 

DE-MPTMS 

Yu,  eta  1.2012 

5.65 

1.2 

5.5 

0.98 

mercapto-functionalized  diatomaceous  earth  suspended  in  agarose 

ag+DE-MPTMS 

Yu,  eta  1.2012 

3.82 

0.75 

4.2 

0.98 

0.99 

0.01 

3.7 

1E+00 

L-cystei  ne-functi  ona  1  i  zed  pol  yethyl  ene  terephtha  late  (second  attempt) 

PET+cys  2 

Duan  and  Lewis  2002 

3.73 

0.62 

4.4 

0.94 

agarose  doped  with  polyethyleneglycol  and  glycerol 

PEGag 

Chari ionet,  et  al.  1996 

2.68 

0.67 

2.9 

0.98 

granular  activated  carbon  suspended  in  PEG-doped  agarose 

PEGag+GAC 

Chari ionet,  et  al.  1996 

3.55 

0.71 

4.1 

0.98 

0.97 

0.03 

3.5 

1E+00 

hiol-SAMMS  suspended  in  PEG-doped  agarose 

PEGag+SAMMS 

Chari  ionet,  et  al.  1996 

3.93 

0.96 

4.0 

1.0 

0.98 

0.02 

3.8 

1E+00 

pol  yethyl  ene  terephtha  late 

PET 

n/a 

3.64 

0.76 

4.0 

0.99 

thiol-SAMMS 

SAM  MS 

n/a 

5.51 

1.3 

5.0 

0.97 
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8.0  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE 
RESEARCH: 


This  1-year  SEED  project  advanced  the  field  of  passive  sampling  on  two  major  fronts:  1) 
Developed  an  engineering  approach  to  extend  in-situ  passive  sampling  to  high  molecular  weight 
compounds  in  sediment  porewater,  and  2)  Explored  an  equilibrium  passive  sampling  approach 
for  MeHg  in  sediment  porewater. 

In  this  research  we  demonstrate  through  laboratory  experiments  and  numerical  modeling  that 
short  periodic  shaking  of  a  passive  sampler  deployed  in  static  sediment  greatly  enhances  the  rate 
of  mass  transfer  and  reduces  the  difference  in  the  extent  of  equilibrium  achieved  compared  to  a 
well  stirred  laboratory  equilibrium.  The  improvement  over  static  sediment  deployment  is 
especially  evident  for  the  high  molecular  weight  PAH  compounds  such  as  benzo(a)pyrene.  We 
also  demonstrate  this  method  for  strongly  hydrophobic  chlorinated  organics  using  PCB 
congeners  in  the  log  Kow  range  of  6-8.  The  accuracy  of  measurement  of  porewater  PCB 
concentrations  is  greatly  enhanced  in  the  vibrating  passive  sampler.  Further  work  is  needed  to 
operationalize  the  vibrating  passive  sampler  concept.  Recommendations  for  future  work  on  the 
vibrating  passive  sampler  include: 

1)  Further  optimization  of  vibration  frequency  through  laboratory  experiments.  In  this  SEED 
project  we  introduced  frequent  vibrations  (once  every  few  minutes).  However,  our  modeling 
results  indicate  that  there  is  an  optimal  frequency  of  vibration  that  can  be  much  lower  (e.g.  once 
per  day).  Further  experimental  work  is  needed  to  develop  a  better  understanding  of  how  the 
optimal  frequency  is  impacted  by  compound  hydrophobicity  and  sediment  characteristics. 

2)  Testing  on  strongly  hydrophobic  compounds  such  as  dioxins  and  furans.  In  this  SEED  project 
we  tested  the  vibrating  sampling  platform  using  field  sediments  with  PAHs  and  highly 
chlorinated  PCBs.  Ultimately,  this  needs  to  be  tested  in  the  laboratory  for  sediments  impacted 
with  dioxins  and  furans.  There  is  very  little  in  the  current  literature  on  the  use  of  passive 
sampling  for  measuring  freely  dissolved  porewater  concentrations  of  these  compounds. 

3)  Confirmation  through  laboratory  experiments  that  the  concept  of  using  freely  dissolved 

concentrations  in  sediment  porewater  to  predict  biouptake  can  be  extended  to  dioxins  and  furans. 

Several  studies  have  correlated  the  freely  dissolved  concentration  of  PAHs  and  PCBs  to  a 
biological  endpoint,  but  such  empirical  data  on  dioxins  and  furans  are  lacking. 

4)  Field  testing  of  the  sampler  platform  along  with  traditional  deployment  of  passive  samplers. 
The  prototype  platforms  developed  in  this  research  are  ready  for  field  deployment.  These  field 
studies  should  ideally  target  strongly  hydrophobic  compounds  such  as  highly  chlorinated  PCBs 
or  dioxins/furans. 
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This  research  also  developed  and  tested  polymer  composites  for  use  as  equilibrium  passive 
samplers  for  MeHg  in  sediment  porewater.  Several  polymer  composites  demonstrated  strong 
sorption  characteristics  and  linear  sorption  behavior  that  we  were  seeking  for  use  in  passive 
sampling.  These  tests  were  also  performed  in  the  presence  of  two  common  complexing  ligands 
in  sediment  porewater:  chloride  and  dissolved  organic  matter.  We  believe  that  initial  results 
presented  in  this  report  show  great  promise  for  pursuing  the  concept  of  equilibrium  passive 
sampling  for  Hg  and  MeHg  in  sediment  porewater. 

Our  study  also  took  the  first  step  toward  demonstrating  that  equilibrium  passive  samplers  for 
MeHg  are  feasible  using  our  proposed  approach,  in  which  materials  with  similar  partition 
coefficients  to  MeHg  in  sediments  are  deployed.  Specifically,  we  developed  several  polymers 
with  sorbent  inclusions  that  meet  our  requirements  for  this  type  of  passive  samplers:  linear 
uptake  isotherms,  and  appropriate  partition  coefficients.  Several  of  these  materials  should  be  able 
to  be  produced  in  thin  films  with  sufficient  capacity  for  MeHg  analysis,  but  this  will  need 
additional  testing. 

Further  work  is  needed  to  advance  the  technology  development  for  passive  sampling  of  MeHg 
and  are  listed  below: 

1)  Complete  testing  of  a  suite  of  selected  polymers  prior  to  testing  against  benthic 
organisms,  specifically: 

a.  Complete  kinetic  testing  over  time  with  simple  salts  (MeHgCl)  and  with 
MeHgDOM  complexes. 

b.  Test  equilibrium  partitioning  into  5-8  selected  passive  sampler  designs,  for  key 
classes  of  synthesized  MeHg  complexes,  including  MeHg-thiols,  MeHg-DOM 
complexes  (with  a  suite  of  DOC  compounds  across  a  range  of  aromaticity  and  S 
content,  see  Graham  et  al.  2012,  2013),  and  colloids  that  contain  both  DOC  and 
sulfide. 

2)  Empirically  determine  which  passive  sampler  sorbents  and  pore  size  predict  MeHg 
uptake  by  organisms. 

•  We  propose  to  use  Leptocheirus  plumulosus  as  a  test  organism  because  it  is  a  mixed- 
type  of  feeder  and  because  of  our  experience  in  developing  MeHg  exposure  systems 
for  this  organism  in  collaboration  with  ERDC. 

•  Directly  test  key  classes  of  MeHg  complexes  for  uptake  by  both  the  animals  and  the 
suite  of  passive  samplers,  using  enriched  stable  isotope  MeHg  spikes  in  artificial 
sediment  test  matrices.  We  have  significant  experience  synthesizing  these  complexes, 
included  HgS  nanoparticulates  and  colloids. 

•  Use  this  information  to  narrow  the  set  of  passive  samplers  to  2-3. 

•  Extend  testing  to  natural  sediments  with  the  select  set  of  passive  samplers.  Include 
measurement  of  total  pore  water  MeHg,  and  variables  needed  for  equilibrium 
modeling  of  MeHg  complexation. 

3)  Develop  enriched  stable  isotope  spikes  as  performance  reference  compounds  in  these 
samplers 
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9.0  RESPONSE  TO  ACTION  ITEM 


Action  Items  Provided  at  Last  In-Progress  Review: 

Please  plan  on  considering  the  following  comments  from  the  Sediment  Review  Panel  over  the 
course  of  this  project.  Recognizing  that  this  project  is  a  SEED  project,  and  is  just  getting 
underway,  we  are  not  expecting  answers  per  se,  but  rather,  will  be  interested  in  seeing  that  these 
items  are  addressed  and  included  in  the  Final  Report. 

1.  The  presentation  and  proposal  is  based  on  mathematical  modeling  that  suggests  that  uptake 
will  be  facilitated  and  shortened  based  on  the  vibratory  disruption  of  the  aqueous  boundary 
layer.  Empirical  evidence  will  be  needed  to  confirm  that  this  is  in  fact  the  case. 

Response:  Yes,  a  major  focus  of  this  research  was  to  perform  laboratory  experiments  to  confirm 
the  theoretical  expectation.  Field  sediments  were  used  in  laboratory  mecosoms  to  test  the 
performance  of  vibrating  samplers  and  compare  with  static  and  well-mixed  deployments. 

Chapter  5  presents  the  results  in  the  form  of  a  peer-review  publication  based  on  experimental  and 
modeling  results  for  PAHs  in  sediments.  Chapter  6  presents  experimental  results  for  highly 
chlorinated  PCBs  that  would  fall  in  the  range  of  most  dioxins  and  furans. 

2.  Multiple  engineering  issues  will  need  to  be  considered  in  the  development  of  this  sampler. 
These  would  include  the  size  and  construction  of  the  sampler,  effectively  vibrating  the  sampling 
media  across  depth  and  width  of  a  sampler  (e.g.,  a  10  x  10  cm  sampler),  whether  the  vibrations 
will  cause  the  sampler  to  come  out  of  the  sediment,  weight  of  the  sampler,  this  to  stay  in  place, 
power  source  (size  and  longevity  of  the  battery).  We  are  confident  that  you  are  aware  of  most  of 
these  issues,  and  are  not  expecting  those  to  be  answered  here,  just  thought  about  and  considered 
in  the  Final  Report. 

Response:  The  engineering  aspects  of  the  development  of  a 
passive  sampler  were  addressed  in  this  research.  As 
demonstrated  in  chapters  5  and  6,  these  samplers  were  able  to 
function  effectively  for  56  h  deployments  in  laboratory 
sediment  mesocosms.  We  went  through  several  iterations  of 
design  and  settled  on  a  final  prototype  shown  on  the  right  that 
was  tested  successfully  in  sediments.  We  were  able  to 
comfortably  measure  low  aqueous  concentrations  of  highly 
chlorinated  PCBs  with  only  0.1  ppm  total  PCBs  in  sediment. 

We  also  calculated  our  power  requirement  and  demonstrate  in 
chapter  5  that  the  motors  need  to  be  active  very  infrequently 
(as  low  as  once  a  day)  and  a  regular  AA  size  battery  is 
sufficient  for  a  1 -month  deployment. 
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10.0  APPENDICES 


10.1  Supporting  information  associated  with  section  5 

Local  equilibrium  model  for  a  static  system 

A  one  dimensional  diffusion  model,  based  on  Fernandez  et  al. 1  was  used  to  describe  the  mass 
transfer  in  a  system  containing  a  PE  strip  with  the  thickness  of  21P  and  sediment  with  thickness  of 
lw  on  both  sides.  Following  the  Fick’s  second  law  of  diffusion,  for  a  PE  strip  with  concentration 
Cpe  and  at  point  x  and  time  t: 

^  =  -lp<x<  lp  (S-l) 

Where  Dpe  is  chemical  diffusivity  in  PE  (cm2/s) 

Assuming  instantaneous  equilibrium  between  pore  water  and  sediment  particles,  transport 
equation  in  pore  water  with  concentration  of  Cw  at  point  x  and  time  t  will  be  as  follows: 

=  D  d~£r  lp  <x  <  lp+lw  and  -  lp  <  x  <  -lp-lw  (S-2) 

Where  p  is  sediment  bulk  density  (g/cm3),  R  is  retardation  factor  and  D  is  the  effective 
diffusivity  (cnr/s) 

R  and  D  are  calculated  as  follows: 


R  =  8/(1  +  P-Kd) 


D  = 


Du 


l-in£2 


(S-3) 

(S-4) 


Where  Dw  is  chemical  diffusivity  in  water  (cnr/s),  Kd  is  sediment-water  partition  coefficient 
(cm3/g)  and  £  is  porosity  (cm3/cm3) 

The  polymer  is  initially  clean  and  sediment  particles  have  initial  concentration  of  SO. 
Cwo=so/Kd  lp  <  x  <  lp+lw  and  -  lp  <  x  < -lp-lw  (S-5) 

CPE  —  0  —lp  <  x  <lp  (S-6) 

At  the  boundary  of  PE-sediment,  Cpe  can  be  determined  using  the  PE-water  partitioning 
coefficient  Kpew  and  pore  water  concentration: 


CpE  —  KPEwCw  x  —  lp  and  x  —  lp  (S-V) 

Continuity  of  flux  is  also  assumed  at  the  boundary: 

D',£S  =  0S  x  =  l„andx  =  -lr  (S-8) 

Due  to  symmetry  the  flux  will  be  zero  at  the  center  of  PE.  Pore  water  concentration  is  equal  to 
the  initial  concentration  far  away  from  the  polymer  at  x=lP+lw  and  does  not  change  by  time: 
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x  =  0 


x  =  L+lw  and  x  =  —lp—lw  (S-10) 


Details  of  model  simulation 


Numerical  solution 


Explicit,  finite-difference  numerical  modeling  techniques,  as  described  by  Crank  2  were  used  to 
solve  the  model’s  equations.  Briefly,  the  explicit  method  is  implemented  by  dividing  the  x  region 
with  the  length  of  L,  into  N  intervals  each  of  width  Ax  such  that  Xi=iAx;  i=l,  2...N  and  NAx=L. 

In  this  model,  polymer  region  with  the  length  of  1P  (half  of  the  polymer  thickness)  was  divided 
into  NP  intervals.  Pore  water  and  sediment  regions  had  the  same  length  of  lw  and  were  divided 
into  Nw  intervals. 

Cl  is  denoted  for  concentration  in  location  idx  and  at  time  jdt,  or  in  other  words  C  (iAx,  j At);  j=l, 
2,3... 

The  numerical  solution  to  model  equations  in  pore  water,  polymer  and  sediment  regions  will  be 
as  Equations  S-ll,  S-12  and  S-13,  respectively: 


Pore  water  region:  l<i<N\. 

11) 


Polymer  region:  Nw+l<i<N 


c/+1-c/ 


_  n  i- 1  i  ^L;+1 

—  Upp  .  ? 

A  Xpp 


(S-12) 


Sediment  region:  N+l<i<N+Nv 
13) 


c!+1-d 


L  =  K KdCi-N  -  Cd 


Where  N=NW+NP,  and  Axpe  and  Axw  are  resolutions  in  PE  and  pore  water  region,  respectively 


We  added  an  extra  point  at  the  polymer-water  boundary  with  the  concentration  of  Cg  in  order  to 
implement  the  finite  difference  approximation  of  second  order  derivatives: 


for  i=Nv 


(S- 14) 


C-+1—C- 


KPEwcg  2C/+l+C/+2 


for  i=Nw 


(S-15) 
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The  concentration  of  the  extra  point  can  be  obtained  using  our  boundary  layer  condition 
(Equation  8): 


CNW+1  KPEwCg  _  Cg  CNW 


A  XPE 


Axu 


From  Equation  S- 16,  Cg  is  calculated  as: 

r 
“ a 


Cn  —  A*  CN  +1  +  B  *  CN 


A  = 


(— ) 
'  dx-n  J 


£*L+RPE-KpEw 
\xw  AxPE 


B 


f— ) 

_ _ dxw ’ _ 

^l+£pjlKpEw 

Axw  A  xPE 


(S-16) 


(S-17) 

(S-18) 

(S-19) 


Substituting  Equation  (S-17)  in  Equations  (S-14)  and  (S-15),  concentrations  at  the  water-polymer 
boundary  can  be  estimated  as  follows: 


lnw  lnw 


At 


Axu 


■  iA  *  CNW+ 1  +  (5  -  2)  *  CNw  +  CWw_i) 


20) 


*  CNW+N 


(S- 


ri+1  —  ('  i 

LJVW  +  1  CJVW+1 

At 


Axpl;2  *  KpEw  *  Qw  +  (A  *  Kppw 


2)  *  C1 


Ww+1 


+  C 


Nw+ 2> 


(S-21) 


Parameter  estimation 


In  order  to  calculate  the  slow  and  fast  desorption  rate  constants  of  pyrene,  chrysene,  and 
benzo(a)pyrene.  Equation  S-22  (Ghosh  et  al.3)  was  fitted  with  desorption  kinetics  data  of 
mentioned  PAHs.  The  desorption  data  were  presented  for  manufactured  gas  plant  site  sediments 
in  Ghosh  et  al.4  (Table  SI). 

-  =  fe^  +  (1  -  f)e~kst  (S-22) 

s0 

Where  f  is  fraction  of  fast  desorbing  pool  in  sediment,  kf  is  first-order  rate  constant  for  fast 
component,  and  ks  is  first-order  rate  constant  for  slow  component. 


Table  S5.1.  Fractional  loss  of  PAHs  from  Harbor  point  sediment  measured  by  Ghosh  et  al.4. 


Time 

(day) 

Pyrene 

Chrysene 

Benzo(a)pyrene 

0 

0 

0 

0 

0.5 

0.22 

0.06 

0.02 

2 

0.37 

0.13 

0.04 

10 

0.49 

0.23 

0.08 

77 


The  equilibrium  concentration  of  PAHs  in  PE  (CpE.eq)  was  determined  by  extracting  the  PE 
samplers  after  56  days  of  deployment  in  the  fully  mixed  PE-sediment  system.  Another  time  point 
measurement  of  PE  concentration  after  77  days  of  deployment  confirmed  that  equilibrium  had 
reached  in  PE  for  all  PAHs  in  56  days  (The  77  data  have  been  compared  to  56  data  in  Table  S3). 
The  pore  water  initial  (bulk)  concentration  (Cwo)  and  sediment-water  partitioning  coefficient  (Kd) 
were  calculated  using  the  measured  sediment  concentration  (Csed)  and  reported  values  for  Kpew  3. 

Cwo  =  ^  (S-23) 

H-PEw 

Kd  =  Cf±  (S-24) 

LW0 

Table  S5.2  summarizes  the  parameter  values  used  in  model  simulations  for  three  PAHs. 


Table  S5.2.  Parameter  values  used  in  model  simulations 


Parameter 

Pyrene 

Chrysene 

Benzo(a)pyrene 

Log  Dpe3 

-9.53 

-10.28 

-10.59 

(cm2/s) 

Log  Dwb  (cnr/s) 

-5.3 

-5.39 

-5.4 

Log  Kpewc 

4.9 

5.6 

6.2 

(L/kg) 

Log  Kd  (L/kg) 

3.3 

4.2 

4.9 

Log  kf  d(s_1) 

-4.63 

-4.76 

-4.63 

Log  Ks  d(s_1) 

-6.53 

-6.73 

-7.2 

s  (cm3/  cm3) 

0.6 

p  (g/  cm3) 

0.6 

a)  Estimated  from  the  correlation  provided  by  Gschwend  5 

b)  Gustafson  and  Dickhut  6 

c)  Lohmann 7 

d)  Measured  by  Ghosh  et  al.  4 


Table  S5.3.  Comparison  of  PAH  concentration  in  PE  in  the  fully  mixed  system  after  56  days 
(C56)  and  77  days  (C77)  of  exposure.  P-value  is  more  than  0.05  for  all  pairs.  This  indicates  that 
the  difference  between  the  77  and  56  data  is  not  significant  at  alpha  level  of  5%  for  all 
congeners. 


Compound 

C56  (ng/mg) 

C77  (ng/mg) 

p-value 

Phenanthrene 

44.5  ±  2.7 

42.0+1.4 

0.23 

Anthracene 

16.2  +  0.8 

15.3  +  0.8 

0.26 

Fluoranthene 

51.5+3.0 

45.3  +  4.8 

0.12 

Pyrene 

57.6  +  3.4 

51.4  +  5.5 

0.17 

Benz(a)anthracene 

31.1  +  1.3 

30.4  +  3.0 

0.72 

Chrysene 

26.1  +  1.4 

26.3  +  2.5 

0.86 
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Benzo(b)fluoranthene 

20.1  ±  1.4 

19.6+1.8 

0.71 

Benzo(k)fluoranthene 

23.2  +  1.4 

21.1  +  1.7 

0.17 

Benzo(a)pyrene 

48.7  ±  1.5 

45.7  +  4.0 

0.25 

Indeno(l,2,3,- 

42.2  ±  0.7 

30.1  +  16 

0.26 

cd)pyrene 

Dibenz(a,h)anthracene 

4.3  ±0.3 

4.5  +  0.5 

0.76 

Quantification  of  difference  between  low  frequency  (5  min  pause)  and  high  frequency  (2min 
pause)  vibrating  data 

PE  concentration  of  pyrene,  chrysene,  and  benzo(a)pyrene  in  high  frequency  and  low  frequency 
vibrating  systems  at  different  exposure  times  have  been  shown  in  Table  S4.  Due  to  loss  of  some 
PE  strips  or  low  recoveries,  duplicate  measurements  were  not  available  for  7  day  and  28  day 
data.  Thus,  only  14  day  and  56  day  data  were  compared  for  significance  of  difference.  As 
indicated  is  Table  S5,  except  for  the  pyrene  concentration  in  day  14,  PE  concentrations  in  the 
two  systems  were  not  significantly  different  at  alpha  level  of  0.05  (p-values  were  more  than 
0.05).  The  maximum  difference  between  two  models’  predictions  is  for  pyrene  in  day  7,  which  is 
only  2%  (Table  S6). 


Table  S5.4.  PAH  concentration  in  PE  in  high  frequency  and  low  frequency  vibrating  systems  at 
different  exposure  times 


Compound 


PE  concentration  in  low  frequency  vibrating  PE  concentration  in  high  frequency  vibrating  system 
system  (pg/g)  (pg/g) 


7  days 

14  days 

28  days 

56  days 

7  days 

14  days 

28  days 

56  days 

Pyrene 

243 

279  +  2.8 

279 

371+0.1 

303 

311+5 

326 

395  +  24 

Chrysene 

69.1 

75  +  19 

95 

113  +  8 

99 

95  +  3 

105 

139  +  8 

Benzo(a)pyrene 

99.8 

137+11 

123 

180  +  2 

152 

142  +  2 

124 

201  +  12 

Table  S5.5.  P-values  calculated  for  the  experimental  data  from  high  and  low  frequency  vibrating 


systems. 


Compound 

P-value 

14  days 

56  days 

Pyrene 

0.014 

0.28 

Chrysene 

0.28 

0.08 

Benzo(a)pyrene 

0.62 

0.13 

Model  prediction: 
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Table  S5.6.  Percent  difference  between  high  and  low  frequency  model  predictions  at  different 
exposure  times. 


Exposure  time 

(day) 

Pyrene 

Chrysene 

Benzo(a)pyrene 

7 

2 

0.69 

0.31 

14 

0.5 

0.48 

0.18 

28 

0.1 

0.2 

0.12 

56 

0 

0 

0 

Quantification  of  difference  between  2  min  and  6  hour  pause  models 

The  percent  differences  between  2-min  and  6-hour  models  have  been  indicated  in  Table  S7  for 
different  exposure  times.  The  maximum  difference  is  approximately  6%. 

Table  S5.7.  Percent  difference  between  2-min  and  6-hour  model  predictions  for  chrysene 


Exposure  time 

(day) 

%  difference 

7 

6.38 

14 

3.95 

28 

1.24 

56 

0.1 

Matlab  codes 


Fast  desorption  model  in  static  system 


outdi 
outf  i 
Lw  = 
Lp  = 
rho  = 
kpew 
kd  = 
phi  = 
Dwate 
Dw  = 
Dp  = 
kf  = 
CwO  = 


r  =  ' Output_Static '  ; 
le  =  strcat (outdir,  ' /  1 
0.1; 

0.00125; 

0.6; 

=  370000; 

10 A4 .2 /rho; 

0.6; 

r  =  4  *  1 0  A  -  6 ; 


' Output_Static . txt ' )  ; 


output  file 


%  half-length  of  water  region  (cm) 

%  half-length  of  polymer  region  (cm) 

%  sediment  bulk  density  (g/cm3) 

%  PE-water  partition  coefficient  for  chrysene (cm3/cm3) 
%  sediment-water  partitioning  coefficient  (cm3/g) 

%  porosity (cm3/cm3) 

%  diffusivity  in  water  (cm2/s) 

Dwater/ ( 1-log (phi A2 )) ;  %  diffusivity  in  pore  water  (cm2/s) 

5*10A-11;  %  diffusivity  in  PE  (cm2/s) 

1.7e-5;  %  fast  desorption  rate  constant  (s-1) 

0.4;  %initial  pore  water  concentration  (ng/cm3) 


Np  =  5 ; 

Nw  =  200; 
dxp  =  Lp/Np; 
dxw  =  Lw/Nw; 


%  discretization  number  on  lp 
%  discretization  number  on  lw 
%  resolution  of  water/sediment  region 
%  resolution  of  water  or  sediment  region 


A  =  (Dp/dxp) / ( (Dw/dxw) + (Dp*kpew/dxp) ) ;  %parameter 
B  =  (Dw/dxw) /( (Dw/dxw) + (Dp*kpew/dxp) ) ;  ^parameter 
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%  initial  condition 

SO  =  CwO*kd;  %initial  sediment  concentration  (ng/g) 

uO  =  [CwO*ones (Nw, 1 ); zeros (Np, 1 ); S0*ones (Nw, 1 )] ;  %  Initial  condition  matrix 

%  defining  the  second-order  differentiation  operator 
N  =  Np+Nw; 

dfl  =  ones (Nw, 1 ) * [ 1 , -2 , 1 ] ; 

Dfl  =  (Dw/ (dxwA2) ) *spdiags (dfl, [-1  0  1] , Nw, Nw) ; 
df2  =  ones (Nw, 1) * [0,  (-rho*kf *kd/phi)  ,  0]  ; 

Df2  =  spdiags (df2 , [-1  0  l],Nw,Nw); 

DfTl  =  Df 1+Df 2 ; 

df2_bbb  =  ones (Np, 1 ) * [ 1 , -2 , 1 ] ; 

Df2bbbb  =  spdiags (df2_bbb, [-1  0  l],Np,Np); 

DfT2  =  [DfTl , zeros (Nw, Np) ; zeros (Np, Nw) , Df2bbbb]  ; 

DfT2(l,:)  =  zeros (1,N) ; 

DfT2 (Nw, : )  =  [zeros(l,Nw-2),Dw/dxwA2,Dw/dxwA2*(B-2)~ 
rho*kf *kd/ phi , Dw/ dxwA2  *A, zeros ( 1 , Np-1 ) ] ; 

DfT2(Nw+l,:)  =  [  zeros ( 1 , Nw-1 ), B*kpew,  kpew*A-2,  1 , zeros ( 1 , Np-2 )] ; 

DfT2 (end,  : )  =  [  zeros ( 1 , Nw+Np-2 ) , 2 ,  -2]; 

Df T2 (Nw+1 : end, : )  =  DfT2 (Nw+1 : end, : ) / (dxpA2 ) *Dp; 

dfsl  =  ones (Nw, 1) * [0, rho*kf/phi, 0] ; 

Dfsl  =  spdiags (dfsl,  [-1  0  1] , Nw,Nw) ; 

Dfsl  =  [Df si ; zeros (Np, Nw) ] ; 

Dfsl(l,:)  =  zeros (l,Nw); 

dfs2  =  ones (Nw, 1 ) * [ 0 , kf *kd, 0 ] ; 

Dfs2  =  spdiags (dfs2 ,  [-1  0  l],Nw,Nw); 

Dfs2(l,:)  =  zeros (l,Nw); 

Dfs2  =  [Df s2 , zeros (Nw, Np) ]  ; 

dfs3  =  ones (Nw, 1 ) * [ 0 , -kf , 0 ] ; 

Dfs3  =  spdiags (dfs3,  [-1  0  l],Nw,Nw); 

Dfs3(l,:)  =  zeros (l,Nw); 

Df T  =  [DfT2,Dfsl;Dfs2,Dfs3]  ; 

DfT  =  sparse (DfT) ; 

Mp  =  @(u_current) ( sum (u_current (Nw+1 : N) ))-( (u_current (Nw+1 ) +u_current (N) ) /2 ) 
%mass  per  unit  area  of  PE  (ng/cm2) 

u_current  (:,1)  =  uO; 

u_inf  =  expm (10A14*DfT) *u0;  %  Equilibrium  concentration  (ng/cm3) 

M_inf  =  Mp(u_inf);  %  Equilibrium  mass  per  unit  area  of  PE  (ng/cm2) 

T_end  =  56*24*3600;  %  deployment  time  (s) 

tvec  =  10000 : 10000 : T_end+1 0000 ; 

for  jtt  =  1 : length (tvec) 

T_current  =  tvec (jtt); 

u_current  (:,jtt)  =  expm (T_current*DfT) *u0; 

M(jtt)  =  Mp  (u_current  (:, jtt) ) ; 


end 


fid  =  f open (' Static . txt w ') ; 
fprintf(fid, ' %2 . 2e\n ' , M/M_inf ) ; 
f close (fid) ; 


Fast  desorption  model  in  high  frequency  vibrating  system 


outdi 
outf  i 
Lw  = 
Lp  = 
rho  = 
kpew 
kd  = 
phi  = 
Dwate 
Dw  = 
Dp  = 
kf  = 
CwO  = 


r  =  ' Output_Vibration '  ; 
le  =  strcat (outdir Output_Vibration . txt ')  ; 

0.1; 

0.00125; 

0.6; 

=  370000; 

10 A4 .2 /rho; 

0.6; 

r  =  4  *  1 0  A  -  6  ; 


%  output  file 

half-length  of  water  region  (cm) 
half-length  of  polymer  region  (cm) 
sediment  bulk  density  (g/cm3) 

PE-water  partition  coefficient  for  chrysene (cm3/  cm3) 
sediment-water  partitioning  coefficient  (cm3/g) 
porosity (cm3/cm3) 
diffusivity  in  water  (cm2/s) 


Dwater/ ( 1-log (phi A2 ) )  ; 


diffusivity  in  pore  water  (cm2/s) 


5*10A-11; 
1 . 4e-5 ; 
0.4; 


%  diffusivity  in  PE  (cm2/s) 

%  fast  desorption  rate  constant  (s-1) 
%initial  pore  water  concentration  (ng/cm3) 


Np  =  5 ; 

Nw  =  200; 
dxp  =  Lp/Np; 
dxw  =  Lw/Nw; 


%  discretization  number  on  lp 
%  discretization  number  on  lw 
%  resolution  of  water/sediment  region 
%  resolution  of  water  or  sediment  region 


A  =  (Dp/dxp) / ( (Dw/dxw) + (Dp*kpew/dxp) ) ;  %parameter 
B  =  (Dw/dxw) /( (Dw/dxw) + (Dp*kpew/dxp) ) ;  %parameter 


%  initial  condition 

SO  =  Cw0*kd;  %initial  sediment  concentration  (ng/g) 

uO  =  [Cw0*ones (Nw, 1) ; zeros (Np, 1) ; S0*ones (Nw, 1) ] ;  %  Initial  condition  matrix 


%  defining  the  second-order  differentiation  operator 
N  =  Np+Nw; 

dfl  =  ones (Nw, 1 ) * [ 1, -2 , 1 ] ; 

Dfl  =  (Dw/ (dxwA2) ) *spdiags (dfl, [-1  0  1] ,Nw,Nw) ; 
df2  =  ones (Nw, 1) * [0, (-rho*kf *kd/phi) , 0] ; 

Df2  =  spdiags (df2 ,  [-1  0  l],Nw,Nw); 

DfTl  =  Df 1+Df 2 ; 

df2_bbb  =  ones (Np, 1 ) * [ 1 , -2 , 1 ] ; 

Df2bbbb  =  spdiags (df2_bbb, [-1  0  l],Np,Np); 

DfT2  =  [DfTl , zeros (Nw, Np) ; zeros (Np, Nw) , Df2bbbb]  ; 

DfT2(l,:)  =  zeros (1,N) ; 

DfT2 (Nw, : )  =  [ zeros ( 1 , Nw-2 ), Dw/dxwA2 , Dw/dxwA2 * (B-2 ) - 
rho*kf *kd/ phi , Dw/ dxwA2 *A, zeros ( 1 , Np-1 ) ]  ; 

DfT2(Nw+l,:)  =  [ zeros ( 1 , Nw-1 ), B*kpew,  kpew*A-2,  1 , zeros ( 1 , Np-2 )] ; 
DfT2 (end,  :)  =  [ zeros ( 1 , Nw+Np-2 ) , 2 ,  —2]; 

Df T2 (Nw+1 : end, : )  =  DfT2 (Nw+1 ; end, : ) / (dxpA2 ) *Dp; 

dfsl  =  ones (Nw, 1) * [0, rho*kf/phi, 0] ; 

Dfsl  =  spdiags (dfsl, [-1  0  1] ,Nw,Nw) ; 

Dfsl  =  [Dfsl; zeros (Np, Nw) ] ; 
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Dfsl(l,:)  =  zeros ( 1 , Nw) ; 

dfs2  =  ones (Nw, 1 ) * [ 0 , kf *kd, 0 ] ; 

Dfs2  =  spdiags (df s2 ,  [-1  0  l],Nw,Nw); 

Dfs2(l,:)  =  zeros (l,Nw); 

Dfs2  =  [Df s2 , zeros (Nw, Np) ]  ; 

dfs3  =  ones (Nw, 1 ) * [ 0 , -kf ,  0 ]  ; 

Dfs3  =  spdiags (dfs3, [-1  0  l],Nw,Nw); 

Dfs3(l,:)  =  zeros (l,Nw); 

Df T  =  [DfT2, Dfsl; Dfs2, Dfs3] ; 

DfT  =  sparse (DfT) ; 

Mp  =  @(u_current) ( sum (u_current (Nw+1 : N) ) ) - ( (u_current (Nw+1 ) +u_current (N) ) /2 ) ; 
%mass  per  unit  area  of  PE  (ng/cm2) 

u_current  (:,1)  =  uO; 

u_inf  =  expm ( 1 0 A 14 *Df T) *u0 ;  %  Equilibrium  concentration  (ng/cm3) 

M_inf  =  Mp(u_inf);  %  Equilibrium  mass  per  unit  area  of  PE  (ng/cm2) 

T_end  =  120;  %  pause  cycle  duration (s) 

tvec  =  60:60:T_end; 
n  =  2 ; 

M(l)  =  0; 

for  k=l : (7*24*3600) /T_end 

for  jtt  =  1 : length (tvec) 

T_current  =  tvec(jtt); 

u_current  (:,jtt)  =  expm (T_current*DfT) *u0; 

M(n)  =  Mp(u_current  (:,jtt)); 
n  =  n+1; 

end 

uO  =  [Cw0*ones (Nw, 1 ) ; u_current (Nw+1 : N, length (tvec) ) ; S0*ones (Nw,  1 ) ]  ; 

%  resetting  the  initial  condition  after  each  pause  cycle 

end 

fid  =  f open (' Vibration . txt w ') ; 
fprintf ( f id, ' %2 . 2e\n ' , M/M_inf ) ; 
f close (fid) ; 


Fully  mixed  model 


function  mixed 
clc; 

outdir  =  'Output  ; 

outfile  =  strcat (outdir,  Output_Mix ') ;  %  output  file 

m  =  50;  %  discretization  number  in  time 

N  =  250;  %  discretization  number  on  x 

L  =  0.0025;  %  Full-length  of  polymer  region  (cm) 

rho  =0.6;  %  sediment  bulk  density  (g/cm3) 
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kpew  =  370000; 

kd  =  10A4.2/rho;  % 

dx  =  L/ (N-l ) ; 

D  =  5*10A-11;  % 

T  =  56*24*3600; 

dt  =  length (T) /m;  % 

%initial  condition 
CO  =  0.4; 

Cl=  kpew*C0; 

CpO  =  [Cl ; zeros (N-2 , 1 ) 


PE-water  partition  coefficient  for  chrysene (cm3/  cm3) 

sediment-water  partitioning  coefficient  (cm3/g) 

resolution  of  polymer  region 

pyrene  diffusivity  in  PE  (cm2/s) 

deployment  time  (sec) 

time  interval 


%  initial  water  concentration  (ng/cm3) 

%  concentration  at  polymer-water  boundary  (ng/cm3) 
Cl];  %  initial  condition  in  polymer  region 


%  defining  the  second-order  differentiation  operator 
Df2  =  D/dxA2 *Df 2_fd_gen (N) ; 

Df2  =  sparse (Df2); 

Cp  =  CpO; 
for  jtt  =  2 :m 

Cp(:,jtt)  =  Cp (1 :N, end) +dt*Df2*  Cp(l:N,end); 

M(jtt)  =  (sum (Cp ( :  ,  jtt) ) - (  Cp(l,jtt)+  Cp (N, jtt) ) / 2)  * (L/ (N-l) ) ; 
end 

fid  =  f open ( ' Output_Mix . txt ' , ' w ' ) ; 
fprintf(fid, '%2.2e\n',M) ; 
f close (fid) ; 
end 

function  Df2  =  Df2_fd_gen (N) 
d  =  [-1;  0;  1]; 

B  =  ones  (N,  1)  *  [1,  -2, 1]  ; 

B (N-l , 1 )  =  0; 

B  (2, 3)  =  0; 

B ( 1 , 2  )  =  0; 

B  (N,  2  )  =  0; 

Df2  =  spdiags (B, d, N, N)  ; 
end 
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Figures 


Figure  S5.1-  Vibrating  system  set  up  before  insertion  into  sediment,  a)  Two  motors  were 
connected  in  parallel  to  a  timer  and  a  power  supply  (2  rechargeable  batteries  1.2V  each).  The 
timers  were  powered  by  a  12  V  battery,  b)  PE  strip  (6mmx20mm)  closed  up  in  stainless  steel 
mesh,  c)  &  d)  Four  meshes  were  attached  to  each  motor  like  fins. 
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8  PE  strips,  5  cm  apart  (static  deployment) 


Vibrating 

motors 


Vibrating 

motors 


Figure  S5.2-  Vibrating  system  set  up  after  insertion  into  sediment.  The  sediment  slurry  was 
placed  in  a  large  glass  tray  (25cmx35cmx6cm).  Four  motors  were  placed  inside  the  sediment. 
Eight  additional  PE  samplers,  enclosed  in  stainless  steel  meshes  without  motors,  were  placed 
inside  the  same  tray  in  a  static  mode  far  from  the  vibrating  motors 
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Figure  S5.3-  Motor  vibration  in  sediment.  Based  on  our  visual  observation  of  vibration  in 
sediment,  the  extent  of  mixing  appears  to  extend  at  least  1  cm  from  the  surface  of  the  fins.  As 
shown  in  the  figure  (taken  with  the  platform  slightly  lifted  from  the  sediment  surface  while 
vibrating),  the  vibration  causes  pressure  waves  to  extend  several  cm  away  from  the  device,  with 
actual  fluidized  mixing  limited  to  about  1  cm  around  the  sampler  fins.  The  vibration  effect  is 
well  beyond  the  thickness  of  the  depletion  layer  at  the  end  of  each  pause  mode,  which  is  only  a 
few  microns. 
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Figure  S5.4-  Experimental  data  for  fractional  uptake  of  fluoranthene  and  four  heavy  PAHs  in  PE 
passive  sampler  in  four  different  systems. 


88 


M/Minit 


1.0  ■ 

0.8  - 


a.Phenanthrene-dlO 


0.4 


0.2 


a  ■  i 


10  20  30  40 


b.Pyrene-dlO 


% ; 


t 


if 


10 


20  30 

Time  (days) 


50 


60 


♦  Fully  mixed 
■  Static 

A  Low  freq  vibration 
x  High  freq  vibration 


Figure  S5.5.  Loss  profile  of  a)  deuterated  phenanthrene  and  b)  deuterated  pyrene  from  PE 
exposed  to  sediment  under  different  mixing  conditions. 
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1.2 


Figure  S5.6-  Comparison  of  three  modeling  approaches  for  the  vibrating  system  for  pyrene. 


Figure  S5.7.  Fractional  uptake  of  three  PAHs  into  PE  versus  energy  consumption  of  one  motor 
after  15  days  of  deployment 
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10.2  Supporting  information  associated  with  section  6 

Parameter  estimation 


PCB  concentrations  in  PE  after  56  days  of  deployment  in  the  fully  mixed  PE-sediment  system 
were  used  as  the  equilibrium  concentration  of  PCBs  in  PE  (CpE,eq).  Since  no  change  was 
observed  in  the  PE  concentrations  after  28  days  of  exposure  in  this  system  (The  28  data  have 
been  compared  to  56  data  in  Table  S5).  The  pore  water  initial  (bulk)  concentration  (Cwo)  and 
sediment-water  partitioning  coefficient  (Kd)  were  calculated  using  the  measured  sediment 
concentration  (Csed)  and  reported  values  for  Kpew.7 

Cwo  =  (S-l) 

kPEw 

Kd  =  Cf±  (S-2) 

LW0 

A  two  compartment  model2  was  used  to  describe  the  desorption  kinetics  of  PCBs  from  sediment 
particles  into  porewater: 

-  =  fe~kft  +  (1  -  f)e~kst  (S-3) 

So 

Where  S  is  PCB  concentration  in  sediment,  So  is  PCB  initial  concentration  in  sediment,  f  is 
fraction  of  fast  desorbing  pool  in  sediment,  kf  is  first-order  rate  constant  for  fast  component,  and 
ks  is  first-order  rate  constant  for  slow  component. 

The  model  was  fitted  to  the  normalized  desorption  data  for  PCB  128,  183,  and  194  in  order  to 
determine  in  f,  kf,  and  ksin  Equation  S-3  (Figure  S6.6).  These  parameters  are  reported  in  Table 
S6.7. 


Table  S6.1  summarizes  the  parameter  values  used  in  model  simulations  for  three  PCBs. 


Model  equations 

Static  system 

In  a  system  containing  a  PE  strip  with  the  thickness  of  21P  and  sediment/porewater  with  thickness 
of  lw  on  both  sides,  sediment  concentration  (S)  changes  following  first  order  kinetics.4'5 

^  =  k(KdCw  -  S)  lp  <x  <  lp+lw  and  -  lp  <  x  <  -lp-lw  (S-4) 

Where  t  is  time  (s),  S  is  chemical  concentration  in  sediment  (ng/g),  Kd  is  sediment-water 
partition  coefficient  (cm3/g),  k  is  first  order  desorption  rate  constant  (s'1)  and  Cw  is  chemical 
concentration  in  water  (ng/cm3) 

For  a  PE  strip  with  concentration  Cpe  and  at  point  x  and  time  t: 

d-rr  =  D^dJSr  -iP<x<ip  (s-5) 

Where  Cpe  is  chemical  concentration  in  PE  (ng/cm3)  and  Dpe  is  chemical  diffusivity  in  PE 
(cm2/s) 
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The  transport  equation  in  porewater  with  concentration  of  Cw  at  point  x  and  time  t  will  be  as 
follows: 


dCw  ^  d2Cw  ( P\  dS 

d  t  dx 2 


L  <  x  <  L+lw  and  —  L  <  x  <  —lv—lw  (S-6) 


£J  gt  -p  ^  A  ^  Lp  1  LW  *-p  '•>  '•>  *-p  <-W 

3\  nnrl  r*  10  t f  \  1  ( r*tvi  3\ 


Where  p  is  sediment  bulk  density  (g/cnr  )  and  s  is  porosity  (cm  '  cmJ) 
D  is  the  diffusivity  in  water  (Dw)  after  correction  for  tortuosity: 


D  = 


Du 


1  -Ins2 


(S-7) 


Substituting  Equation  (1)  in  Equation  (3),  the  transport  equation  in  porewater  can  be  re-written 
as: 


9CW 

dt 


=  D?^-(7 )Ww-S) 


(S-8) 


-lp  <x  <  Ip  (S-9) 

Zp  <  x  <  lp+lw  and  —  lp  <  x  < —lp—lw  (S-10) 


Initial  conditions 

The  polymer  was  initially  clean  and  porewater  was  assumed  to  be  in  equilibrium  with  sediment. 

CPE  =  0 

Q/o  =  ^0  /Kd 

Where  So  is  initial  chemical  concentration  in  sediment  (ng/g)  and  Cwo  is  initial  chemical 
concentration  in  water  (ng/cm3) 

Boundary  conditions 

Continuity  of  flux  and  equilibrium  condition  was  assumed  at  the  PE-water  boundary  as  done  in 
previous  work.8 


n  dcPE  _  n  dCw 
UPE  ~  —  u  ~ 


dxPE 


dx,, 


x  —  Zp  and  x  —  —lp,  t  >  0  (S-ll) 


x  =  Zp  and  x  =  —lp  (S-12) 


CpE  KppyyCyy  A  t-p  LAI  l- LA  A  Lp 

Where  Kpew  is  the  PE-water  partition  coefficient  (cm3/cm3) 

Due  to  symmetry  the  flux  will  be  zero  at  the  center  of  the  PE  sheet.  Porewater  concentration  is 
equal  to  the  initial  concentration  far  away  from  the  polymer  at  x=lP+lw  and  does  not  change  ov 


over 


time: 

9  CpE  _ 


dx 


=  0 


x  =  0 


(S-12) 


dCw  _dS  _ 
dt  dt 


X 


Zp+Zw  and  x  —  —lp—lw  (S-13) 
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Vibrating  system 

In  the  vibrating  system,  when  the  motor  is  in  pause  mode,  the  mass  transfer  is  similar  to  the  static 
mode  and  sediment  and  porewater  concentration  in  the  vicinity  of  the  polymer  depletes  with 
time.  Every  time  the  motor  vibrates,  the  sediment  and  porewater  at  the  vicinity  of  the  polymer  is 
mixed  up.  We  assumed  that  this  mixing  is  enough  to  increase  sediment  concentration  to  the 
initial  concentration  in  sediment  (So).  Porewater  concentration  right  after  each  vibration  pulse 
will  also  increase  to  the  initial  concentration: 

Cwq  t  S  Sq  Ip  X  CLTld  Ip  ^  X  ^  Ip  (S-14) 

Fully  mixed  system 

Diffusion  in  a  polymer  with  thickness  of  21P  follows  Fick’s  second  law: 

d-jT  =  DPEeJ£r  -lp<x<lp  (S-15) 

The  passive  sampler  was  assumed  to  be  initially  clean.  Since  the  system  is  perfectly  mixed,  the 
porewater  concentration  remains  constant  and  equal  to  the  initial  value  (So/Kd)  during  the 
deployment  time.  The  boundary  condition  is  defined  as: 

Cp£  —  KppwCw 0  x  =  Ip  cmd  x  =  Ip  (S-16) 


Details  of  model  simulation 


Numerical  solution 

Explicit,  finite-difference  numerical  modeling  techniques,  as  described  by  Crank  2  were  used  to 
solve  the  model’s  equations.  Briefly,  the  explicit  method  is  implemented  by  dividing  the  x  region 
with  the  length  of  L,  into  N  intervals  each  of  width  Ax  such  that  Xi=iAx;  i=l,  2...N  and  NAx=L. 

In  this  model,  polymer  region  with  the  length  of  1P  (half  of  the  polymer  thickness)  was  divided 
into  NP  intervals.  Porewater  and  sediment  regions  had  the  same  length  of  lw  and  were  divided 
into  Nw  intervals. 


C(  is  denoted  for  concentration  in  location  idx  and  at  time  jdt,  or  in  other  words,  C  (iAx,  j At); 
j=l,  2,  3... 

The  numerical  solution  to  model  equations  in  pore  water,  polymer  and  sediment  regions  will  be 
as  Equations  S-ll,  S-12  and  S-13,  respectively: 


Porewater  region:  l<i<Nw 
Polymer  region:  Nw+l<i<N 


c/+1— c. 
At 


D  -  (f)  HKtCt  -  CUN)  (S-17) 


c(+1-c] 


At 


l-  =  D 


rJ  _2 r-^+r^ 
Li- 1 


PE  A xPE2 


(S-18) 
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Sediment  region:  N+l<i<N+Nw  — - L  =  k(KdCi_N  —  C) )  (S- 

19) 

Where  N=NW+NP,  and  Axpe  and  Axw  are  resolutions  in  PE  and  pore  water  region,  respectively 


We  added  an  extra  point  at  the  polymer-water  boundary  with  the  concentration  of  Cg  in  order  to 
implement  the  finite  difference  approximation  of  second  order  derivatives: 


c/+1-c/ 


At 


C-+1—C- 


-  —  D 


C /  , -2 C/+Cn 


At 


D 


PE  ' 


Axw 

KpEw  (•/+ 1 


~  W  k(KdCi  —  Ci+N) 


hXpE 


for  i=Nv 

for  i=Nw 


(S-20) 

(S-21) 


The  concentration  of  the  extra  point  can  be  obtained  using  our  boundary  layer  condition 
(Equation  8): 


CNW+1  KpEwCg  _  p 


~CNU 


Axps 


Ax,, 


(S-22) 


From  Equation  S-16,  Cg  is  calculated  as: 


-  A  *  CWw+1  +  B  *  CNw 

(S-23) 

A  _  W 

Axw  A  xPE  PEw 

(S-24) 

(—) 
p  _  dxwJ 

~DW  DpE 

Axw  AxPE  PEw 

(S-25) 

Substituting  Equation  (S-17)  in  Equations  (S-14)  and  (S-15),  concentrations  at  the  water-polymer 
boundary  can  be  estimated  as  follows: 


rt+l_r7 
LNw  LNw 


At 


{a  *  CNw+1  +  (B  -  2)  .  CNw  +  CNw.,)  -  (^)  .  CNw  +  (fj  ,  C, 


Nw+N 


(S- 


26) 


cj+1  —pi 

CJVW+1  CJVW  +  1 

At 


■  (A  *  KPEw  —  2)  *  CNw+1  +  CNw+2 ) 


(S-27) 
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Matlab  codes 


Fast  desorption  model  in  static  system  for  PCB  128 

outdir  =  ' Output_Static '  ; 

outfile  =  strcat (outdir, 'Output_Static.txt ') ;  %  output  file 

Lw  =0.1;  %  half-length  of  water  region  (cm) 

Lp  =  0.00125;  %  half-length  of  polymer  region  (cm) 

rho  =0.9;  %  sediment  bulk  density  (g/cm3) 

kpew  =  10A6.6*0.93;  %  PE-water  partition  coefficient  for  chrysene (cm3/cm3) 

kd  =  10A5.47;  %  sediment-water  partitioning  coefficient  (cm3/g) 

phi  =0.7;  %  porosity (cm3/cm3) 

Dwater  =  10A-5.41;  %  diffusivity  in  water  (cm2/s) 

Dw  =  Dwater/ ( 1-log (phi A2 ))  ;  %  diffusivity  in  pore  water  (cm2/s) 

Dp  =  10A-9.35;  %  diffusivity  in  PE  (cm2/s) 

kf  =  3.33e-5;  %  fast  desorption  rate  constant  (s-1) 

CwO  =  0.4;  %initial  pore  water  concentration  (ng/cm3) 


Np  =  5 ; 

Nw  =  200; 
dxp  =  Lp/Np; 
dxw  =  Lw/Nw; 


%  discretization  number  on  lp 
%  discretization  number  on  lw 
%  resolution  of  water/sediment  region 
%  resolution  of  water  or  sediment  region 


A  =  (Dp/dxp) / ( (Dw/dxw) + (Dp*kpew/dxp) ) ;  %parameter 
B  =  (Dw/dxw) /( (Dw/dxw) + (Dp*kpew/dxp) ) ;  %parameter 


%  initial  condition 

SO  =  Cw0*kd;  %initial  sediment  concentration  (ng/g) 

uO  =  [Cw0*ones (Nw, 1) ; zeros (Np, 1) ; S0*ones (Nw, 1) ] ;  %  Initial  condition  matrix 

%  defining  the  second-order  differentiation  operator 
N  =  Np+Nw; 

dfl  =  ones (Nw, 1 ) * [ 1, -2 , 1 ] ; 

Dfl  =  (Dw/ (dxwA2) ) *spdiags (dfl, [-1  0  1] ,Nw,Nw) ; 
df2  =  ones (Nw, 1) * [0, (-rho*kf *kd/phi) , 0] ; 

Df2  =  spdiags (df2 ,  [-1  0  l],Nw,Nw); 

DfTl  =  Df 1+Df 2 ; 

df2_bbb  =  ones (Np, 1 ) * [ 1 , -2 , 1 ] ; 

Df2bbbb  =  spdiags (df2_bbb, [-1  0  l],Np,Np); 

DfT2  =  [DfTl , zeros (Nw, Np) ; zeros (Np, Nw) , Df2bbbb] ; 

DfT2(l,:)  =  zeros (1,N) ; 

DfT2 (Nw, : )  =  [ zeros ( 1 , Nw-2 ), Dw/dxwA2 , Dw/dxwA2 * (B-2 ) - 
rho*kf *kd/ phi , Dw/ dxwA2  *A, zeros ( 1 , Np-1 ) ] ; 

DfT2(Nw+l,:)  =  [ zeros ( 1 , Nw-1 ), B*kpew,  kpew*A-2,  1 , zeros ( 1 , Np-2 )] ; 

DfT2 (end,  : )  =  [  zeros ( 1 , Nw+Np-2 ) , 2 ,  -2]; 

Df T2 (Nw+1 : end, : )  =  DfT2 (Nw+1 : end, :)/ (dxpA2 ) *Dp; 

dfsl  =  ones (Nw, 1) * [0, rho*kf/phi, 0] ; 

Dfsl  =  spdiags (dfsl, [-1  0  1] ,Nw,Nw) ; 

Dfsl  =  [Df si ; zeros (Np, Nw) ] ; 

Dfsl(l,:)  =  zeros (l,Nw) ; 


dfs2  =  ones (Nw, 1 ) * [ 0 , kf *kd, 0 ] ; 

Dfs2  =  spdiags (dfs2 ,  [-1  0  l],Nw,Nw); 
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Dfs2(l,:)  =  zeros (l,Nw); 

Dfs2  =  [Df s2 , zeros (Nw, Np) ] ; 

dfs3  =  ones (Nw, 1 ) * [ 0 , -kf , 0 ] ; 

Dfs3  =  spdiags (dfs3, [-1  0  l],Nw,Nw); 

Dfs3(l,:)  =  zeros (l,Nw); 

Df T  =  [Df T2 , Df si ; Df s2 , Df s3 ] ; 

DfT  =  sparse (DfT); 

Mp  =  @(u_current)  ( sum (u_current (Nw+1 : N) ) ) - ( (u_current (Nw+1 ) +u_current (N) ) /2 )  ; 
%mass  per  unit  area  of  PE  (ng/cm2) 

u_current  (:,1)  =  uO; 

u_inf  =  expm ( 1 0 A 14 *Df T) *u0 ;  %  Equilibrium  concentration  (ng/cm3) 

M_inf  =  Mp(u_inf);  %  Equilibrium  mass  per  unit  area  of  PE  (ng/cm2) 

T_end  =  56*24*3600;  %  deployment  time  (s) 

tvec  =  10000 : 10000 : T_end+10000  ; 

for  jtt  =  1 : length (tvec) 

T_current  =  tvec (jtt); 

u_current  (:,jtt)  =  expm (T_current*DfT) *u0; 

M(jtt)  =  Mp  (u_current  (:,jtt)); 

end 

fid  =  f open (' Static . txt w ') ; 
fprintf ( f id,  ' %2 . 2e\n ' , M/M_inf ) ; 
f close (fid) ; 


Fast  desorption  model  in  two  min-pause  vibrating  system  for  PCB  128 

outdir  =  ' Output_Vibration ' ; 

outfile  =  strcat (outdir Output_Vibration . txt ') ;  %  output  file 

Lw  =0.1;  %  half-length  of  water  region  (cm) 

Lp  =  0.00125;  %  half-length  of  polymer  region  (cm) 

rho  =0.6;  %  sediment  bulk  density  (g/cm3) 

kpew  =  10^6.6*0.93;  %  PE-water  partition  coefficient  for  chrysene (cm3/cm3) 

kd  =  10A5.47;  %  sediment-water  partitioning  coefficient  (cm3/g) 

phi  =0.7;  %  porosity (cm3/cm3) 

Dwater  =  10A-5.41;  %  diffusivity  in  water  (cm2/s) 

Dw  =  Dwater/ ( 1-log (phi A2 ))  ;  %  diffusivity  in  pore  water  (cm2/s) 

Dp  =  10A-9.35;  %  diffusivity  in  PE  (cm2/s) 

kf  =  3.33e-5;  %  fast  desorption  rate  constant  (s-1) 

CwO  =  0.4;  %initial  pore  water  concentration  (ng/cm3) 


Np  =  5; 

Nw  =  200; 
dxp  =  Lp/Np; 
dxw  =  Lw/Nw; 


%  discretization  number  on  lp 
%  discretization  number  on  lw 
%  resolution  of  water/sediment  region 
%  resolution  of  water  or  sediment  region 


A  =  (Dp/dxp) / ( (Dw/dxw) + (Dp*kpew/dxp) ) ;  %parameter 
B  =  (Dw/dxw) /( (Dw/dxw) + (Dp*kpew/dxp) ) ;  %parameter 
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%  initial  condition 

SO  =  CwO*kd;  %initial  sediment  concentration  (ng/g) 

uO  =  [CwO*ones (Nw, 1 ); zeros (Np, 1 ); S0*ones (Nw, 1 )] ;  %  Initial  condition  matrix 

%  defining  the  second-order  differentiation  operator 
N  =  Np+Nw; 

dfl  =  ones (Nw, 1 ) * [ 1 , -2 , 1 ] ; 

Dfl  =  (Dw/ (dxwA2) ) *spdiags (dfl, [-1  0  1] , Nw,Nw) ; 
df2  =  ones (Nw, 1) * [0, (-rho*kf *kd/phi) , 0] ; 

Df2  =  spdiags (df2 ,  [-1  0  l],Nw,Nw); 

DfTl  =  Df 1+Df 2 ; 

df2_bbb  =  ones (Np, 1 ) * [ 1 , -2 , 1 ] ; 

Df2bbbb  =  spdiags (df2_bbb, [-1  0  l],Np,Np); 

DfT2  =  [DfTl , zeros (Nw, Np) ; zeros (Np, Nw) , Df2bbbb]  ; 

DfT2(l,:)  =  zeros (1,N) ; 

DfT2  (Nw,  :  )  =  [zeros(l,Nw-2),Dw/dxwA2,Dw/dxwA2*(B-2)~ 
rho*kf *kd/ phi , Dw/ dxwA2  *A, zeros ( 1 , Np-1 ) ] ; 

DfT2(Nw+l,:)  =  [ zeros ( 1 , Nw-1 ), B*kpew,  kpew*A-2,  1 , zeros ( 1 , Np-2 )] ; 

DfT2 (end,  : )  =  [ zeros ( 1 , Nw+Np-2 ) , 2 ,  -2]; 

Df T2 (Nw+1 : end, : )  =  DfT2 (Nw+1 :end, : ) / (dxpA2) *Dp; 

dfsl  =  ones (Nw, 1) * [0, rho*kf/phi, 0] ; 

Dfsl  =  spdiags (dfsl ,  [-1  0  l],Nw,Nw); 

Dfsl  =  [Df si ; zeros (Np, Nw) ] ; 

Dfsl(l,:)  =  zeros ( 1 , Nw) ; 

dfs2  =  ones (Nw, 1 ) * [ 0 , kf *kd, 0 ] ; 

Dfs2  =  spdiags (dfs2 ,  [-1  0  l],Nw,Nw); 

Dfs2(l,:)  =  zeros (l,Nw); 

Dfs2  =  [Df s2 , zeros (Nw, Np) ]  ; 

dfs3  =  ones (Nw, 1 ) * [ 0 , -kf , 0 ] ; 

Dfs3  =  spdiags (dfs3,  [-1  0  l],Nw,Nw); 

Dfs3(l,:)  =  zeros (l,Nw); 

Df T  =  [DfT2, Dfsl; Dfs2, Dfs3] ; 

DfT  =  sparse (DfT); 

Mp  =  @(u_current)  ( sum  (u_current (Nw+1 : N) ))-( (u_current (Nw+1 ) +u_current (N)  )  /2  ) 
%mass  per  unit  area  of  PE  (ng/cm2) 

u_current  (:,1)  =  uO; 

u_inf  =  expm (10A14*DfT) *u0;  %  Equilibrium  concentration  (ng/cm3) 

M_inf  =  Mp(u_inf);  %  Equilibrium  mass  per  unit  area  of  PE  (ng/cm2) 

T_end  =  120;  %  pause  cycle  duration (s) 

tvec  =  60:60:T_end; 
n  =  2 ; 

M(l)  =  0; 

for  k=l : (7*24*3600) /T_end 

for  jtt  =  1 : length (tvec) 

T_current  =  tvec(jtt); 

u_current  ( : ,  jtt)  =  expm (T_current*DfT) *u0; 


M(n)  =  Mp(u_current  (:,jtt)); 
n  =  n+1; 

end 

uO  =  [CwO*ones (Nw, 1 ) ; u_current (Nw+1 : N, length (tvec) ) ; S0*ones (Nw,  1 ) ]  ; 
%  resetting  the  initial  condition  after  each  pause  cycle 

end 

fid  =  f open (' Vibration . txt w ') ; 
fprintf(fid, ' %2 . 2e\n ' , M/M_inf ) ; 
f close (fid) ; 


Well-mixed  model  for  PCB  128 


function  mixed 
clc; 

outdir  =  'Output1; 
outfile  =  strcat (outdir, 
m  =  50; 

N  =  250; 

L  =  0.0025; 
rho  =  0.6; 
kpew  =  10A6. 6*0. 93; 
kd  =  10A5.47;  % 

dx  =  L/ (N-l) ; 

D  =  10--9.35; 

T  =  56*24*3600; 
dt  =  length (T)/m; 


/ '  ,  ' Output_Mix ' ) ;  %  output  file 

%  discretization  number  in  time 
%  discretization  number  on  x 
%  Full-length  of  polymer  region  (cm) 

%  sediment  bulk  density  (g/cm3) 

%  PE-water  partition  coefficient  for  chrysene (cm3/cm3) 
sediment-water  partitioning  coefficient  (cm3/g) 

%  resolution  of  polymer  region 
%  pyrene  diffusivity  in  PE  (cm2/s) 

%  deployment  time  (sec) 

%  time  interval 


%initial  condition 
CO  =  0.4; 

Cl—  kpew*C0; 

CpO  =  [Cl; zeros (N-2 , 1 ) ; Cl] ; 


initial  water  concentration  (ng/cm3) 
concentration  at  polymer-water  boundary  (ng/cm3) 
%  initial  condition  in  polymer  region 


%  defining  the  second-order  differentiation  operator 
Df2  =  D/dxA2 *Df 2_fd_gen (N) ; 

Df2  =  sparse (Df2); 

Cp  =  CpO; 
for  jtt  =  2 :m 

Cp(:,jtt)  =  Cp (1 :N, end) +dt*Df2*  Cp(l:N,end); 

M(jtt)  =  (sum (Cp ( : , jtt) ) - (  Cp(l,jtt)+  Cp (N, jtt) ) /2) * (L/ (N-l) ) ; 
end 

fid  =  f open (' Output_Mix . txt w ') ; 
fprintf (fid, ' %2 . 2e\n ' , M) ; 
f close (fid) ; 
end 

function  Df2  =  Df2_fd_gen (N) 


d  = 

[-1; 

0; 

B  = 

ones 

(N, 

B  (N- 

-1,1) 

= 

B  (2 , 

3)  = 

0; 

B  ( 1 , 

2)  = 

0; 

B(N, 

2)  = 

0; 
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Df  2 
end 


spdiags  (B,  d,  N,  N)  ; 


Tables 

Table  S6.1.  Parameter  values  used  in  model  simulations 


Parameter 

PCB  128 

PCB  183 

PCB  194 

Log  Dpe3 
(cm2/s) 

-9.35 

-9.53 

-9.7 

Log  Dwb  (cm2/s) 

-5.41 

-5.44 

-5.47 

Log  Kpew3 
(L/kg) 

6.6 

7.1 

7.3 

Log  Kd  (L/kg) 

5.47 

6.13 

6.75 

kf  (s'1) 
s  (cm3/cm3) 
p  (g/cm3) 

3.33xl0"5 

0.7 

0.9 

3. 06x1 0'5 

2.78xl0'5 

e)  Lohmann  1 

f)  Gustafson  and  Dickhut  3 


Table  S6.2.  Comparison  of  PCB  concentration  in  PE  in  the  fully  mixed  system  after  28  days 
(C28)  and  56  days  (C56)  of  exposure.  P-value  is  more  than  0.05  for  all  pairs.  This  indicates  that 
the  difference  between  the  28  and  56  data  is  not  significant  at  alpha  level  of  5%  for  all 
congeners. 


Compound 

C28  (ng/g) 

C56  (ng/g) 

P-values 

PCB  31 

60.2  ±  8.3 

58.4+10.9 

0.71 

PCB  70+  76 

56.7  +  10.1 

34.9  +  6.01 

0.52 

PCB  99 

23.91+2.8 

20.5+4.5 

0.36 

PCB  146 

61.5+23.8 

50.8  +  8.9 

0.51 

PCB  153 

208.6  +  84.6 

178.7  +  35.1 

0.60 

PCB  132 

52.6  ±  27.9 

38.0  +  8.4 

0.47 

PCB  163+PCB  138 

294.3  +  156.3 

221.1  +46.5 

0.52 

PCB  158 

40.1  ±31.2 

26.5  +  7.4 

0.54 

PCB  178+PCB  129 

48.8  +  19.9 

42.7  +  8 

0.66 

PCB  187+PCB  182 

109.9+15.3 

118.3  +  18 

0.57 

100 


PCB  183 

63.0+15.9 

62.3  +  12.3 

0.95 

PCB  128 

26.2  ±  20.8 

16.3  +  5.3 

0.50 

PCB  185 

11.8  +  1.8 

12.4  +  2.2 

0.74 

PCB  201 

115.3  +  17 

136.6  +  22.1 

0.25 

PCB  203+PCB  196 

129.5  +  22 

155.8  +  27.7 

0.26 

PCB  208+PCB  195 

26.0  +  5.4 

30.8  +  5.4 

0.34 

PCB  194 

54.3  +  9.2 

64.9  +  14.6 

0.36 

Table  S6.3.  PCB  concentration  in  static  PE  at  different  exposure  times  from  experiment  1  and  2. 


Compound 

Experiment  1 

Experiment  2 

7  d 

14  d 

28  d 

56  d 

7  d 

14  d 

28  d 

56  d 

PCB  31 

19.9  +  2.3 

18.5  +  2.4 

26.1+5.6 

25.7  +  4.6 

15.7+1.9 

16.8  +  1.8 

18.5+2.2 

18.1  +2.1 

PCB  70+  76 

12.1  +  1.9 

11.7+1.5 

16.1  +2.4 

17.7  +  3.1 

10.3  +  1.2 

11.1  +  1.3 

12.4+1.4 

12.3  +  1.4 

PCB  99 

3.5  +0.3 

3.7  +  0.6 

8.1  +  1.7 

7.1  +  1.5 

3.6  +  0.6 

4.6 +  0.5 

5.4  +  0.03 

6.2  +  0.5 

PCB  146 

15.4  +  0.6 

17.3  +  1.6 

22.0  +  0.4 

26.5  +  2.7 

14.9  +  0.4 

22.3  +  1.3 

21.8+2.2 

21.5  +  1.1 

PCB  153 

31.5  +  6.5 

36.0  +  3 

45.8+0.9 

54.0  +  5.8 

25.7+1.4 

28.2  +  3.7 

42.1  +8.3 

42.3  +  4.6 

PCB  132 

6.7+  1.4 

9.0 +  0.9 

11.3+0.3 

13.6+1.1 

4.5  +  0.4 

7.2  +  1 

10.4+1.9 

10.9+1.1 

PCB  163+PCB  138 

38.6  +  0.2 

44.7  +  4.6 

56.7  +  1 

66.4  +  6.4 

34.4  +  2.6 

29.2  +  4.8 

49.9  +  9.9 

52.2  +  5.4 

PCB  158 

3.1  +0.4 

3.9 +  0.3 

5.1  +0.1 

5.9 +  0.7 

2.6  +  0.4 

2.9  +  0.9 

4.6+  1.2 

4.8  +0.4 

PCB  178+PCB  129 

5.2 +  0.3 

8.6 +  0.9 

11.0  +  0.8 

12.7+1.7 

4.3  +  0.6 

7.9 +  5.3 

8.8  +  1.4 

9.8  +  1.3 

PCB  187+PCB  182 

25.4+1.7 

27.0  +  2.7 

33.2  +  1.9 

36.6  +  4.5 

22.7+1.7 

25.3  +  2.6 

31.8+4.8 

33.8  +  4.8 

PCB  183 

8.3 +  0.3 

10.6  +  1 

17.5  +  1 

19.1+2.4 

7.2  +  0.5 

9.2+  1.3 

13.9  +  2.1 

14.6  +  2.2 

PCB  128 

2.5  +  0.1 

2.9  +  0.2 

3.8 +  0.1 

4.7  +  0.5 

1.5 +  0.1 

1.9 +  0.4 

3.6 +  0.8 

3.7 +  0.4 

PCB  185 

3.8 +  0.1 

4.0 +  0.4 

4.7  +  0.4 

5.1  +0.7 

2.1  +0.3 

2.7 +  0.4 

3.8 +  0.7 

3.9  +  0.6 

PCB  201 

24.4+1.6 

26.3  +  3.2 

41.2  +  3.1 

43.8  +  5.7 

21.5  +  2.1 

21.3  +  3.5 

35.1+6.5 

36.5  +  5.9 

PCB  203+PCB  196 

26.4  +  0.5 

28.2  +  3 

45.5  +  3.7 

48.6  +  5.9 

23.5  +  2.4 

26.9  +  4 

39.4  +  6.7 

40.5  +  6.5 

PCB  208+PCB  195 

5.7  +  0.7 

6.9 +  0.9 

8.1  +0.6 

8.8  +  1.1 

3.5  +  0.3 

3.7 +  0.6 

6.8  +  1.2 

7.4+  1.1 

PCB  194 

12.0  +  0.8 

18.1  +2 

21.0  +  2.3 

21.6  +  2.2 

10.0+1.5 

8.5 +  0.7 

17.9  +  3.4 

18.0  +  2.9 
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Table  S6.4.  P-values  calculated  for  the  PCB  concentrations  in  static  PE  from  experiment  1  and 
experiment  2.  P-value  is  more  than  0.05  for  all  pairs.  This  indicates  that  the  difference  between  PCB 
uptake  in  two  experiments  is  not  significant  at  alpha  level  of  5%  for  all  congeners. 


Compound 

P-value 

7  d 

14  d 

28  d 

56  d 

PCB  31 

0.1 

0.45 

0.17 

0.12 

PCB  70+  76 

0.38 

0.64 

0.2 

0.1 

PCB  99 

0.7 

0.18 

0.13 

0.51 

PCB  146 

0.19 

0.07 

0.87 

0.24 

PCB  153 

0.2 

0.08 

0.51 

0.14 

PCB  132 

0.27 

0.11 

0.5 

0.10 

PCB  163+  138 

0.1 

0.07 

0.35 

0.12 

PCB  158 

0.34 

0.18 

0.56 

0.26 

PCB  178+  129 

0.11 

0.83 

0.1 

0.18 

PCB  187+  182 

0.11 

0.07 

0.68 

0.57 

PCB  183 

0.06 

0.27 

0.08 

0.17 

PCB  128 

0.06 

0.05 

0.73 

0.15 

PCB  185 

0.08 

0.068 

0.14 

0.18 

PCB  201 

0.09 

0.11 

0.26 

0.30 

PCB  203+  196 

0.08 

0.35 

0.28 

0.29 

PCB  208+  195 

0.15 

0.05 

0.2 

0.27 

PCB  194 

0.07 

0.14 

0.31 

0.21 

Table  S6.5.  PRC  fraction  remaining  in  PE  (fpRc)  deployed  in  static  system  after  different  exposure  times. 


Compound 

Experiment  1 

Experiment  2 

7  d 

14  d 

28  d 

56  d 

7  d 

14  d 

28  d 

56  d 

PCB  29 

0.7  +  0.04 

0.6  +  0.03 

0.45  +  0.03 

0.32  +  0.03 

0.63  +  0.02 

0.60  +  0.02 

0.46  +  0.02 

0.32  +  0.02 

PCB  69 

0.75  +  0.03 

0.67  +  0.03 

0.54  +  0.02 

0.39  +  0.02 

0.73+0.01 

0.71+0.01 

0.56  +  0.00 

0.39  +  0.04 
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PCB  155 

0.89  +  0.02  0.83  +  0.04 

0.78  +  0.07  0.64  +  0.00  0.82  +  0.03  0.81  +0.01  0.76  +  0.09  0.64  +  0.06 

PCB  192 

0.93  +  0.05  0.9  +  0.05 

0.88  +  0.06  0.74  +  0.01  0.94  +  0.01  0.86  +  0.04  0.90  +  0.04  0.74  +  0.01 

Table  S6.6.  P-values  calculated  for  fpRc  in  static  system  from  experiment  1  and  experiment  2.  P-value  is 
more  than  0.05  for  all  pairs.  This  indicates  that  the  difference  between  PRC  loss  in  two  experiments  is 
not  significant  at  alpha  level  of  5%  for  all  congeners. 

Compound 

P-value 

7  d 

14  d  28  d  56  d 

PCB  29 

0.3 

0.84  0.90  0.22 

PCB  69 

0.56 

0.94  0.21  0.30 

PCB  155 

0.13 

0.54  0.41  0.48 

PCB  192 

0.85 

0.32  0.17  0.21 

Table  S6.7.  logPEw  of  PCB  congeners.  For  PCB  (70+76),  PCB  (82+151),  PCB  (163+138),  PCB 
(182+187)  the  average  log  Kpew  of  the  two  PCBs  was  used 

PCB  congener  logKowa 

10gKPEwb 

31 

5.67 

5.43 

70+76 

6.2 

6.06 

76 

6.13 

5.97 

99 

6.39 

6.28 

82+151 

6.2 

6.06 

118 

6.74 

6.69 

128 

6.74 

6.69 

132 

6.58 

6.50 

138 

6.99 

6.99 

151 

6.32 

6.20 

153 

6.92 

6.91 

158 

7.02 

7.02 

163+138 

6.83 

6.80 

177 

7.08 

7.09 

180 

7.36 

7.42 

183 

7.2 

7.24 

185 

7.11 

7.13 

182+187 

7.2 

7.24 

187 

7.17 

7.20 

194 

7.8 

7.94 
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195 

7.56 

7.66 

196 

7.65 

7.77 

201 

7.62 

7.73 

203 

7.65 

7.77 

208 

7.71 

7.84 

a)  Hawker  and  Connell  6 

b)  From  the  correlation  provided  by  Ghosh  et  al.  8 


104 


Figures 


20 

15 

10 

5 

0 


240 


§  180 


-g  120 

<5 

o 

§  60 
O 


PCB  128 


i 


\ 


i  i 


0  10  20  30  40 

PCB153 


*  * 


0  10  20  30  40 


PC 

B183 

80 

I 

> 

60 

T  t  ' 

► 

40 

lit 

20 

*•  *  j 

i 

0  1 

i - 

10  20  30  40 

Time  (d) 


* 


50 


60 


50 


60 


* 


50 


60 


PCB  146 


IS 

c 


c 

o 


c 

o 

o 

c 

o 

O 


60 

40 

20 

0  I 

30 


1 


it 


t 

A 

i 


10 


20 


30 


40 


50 


■■5  20 


10 


PCB  158  (Hexa) 


n 

o 


i  i 


k 


10  20  30  40 

Time  (d) 


60 


•  Static 

•  Two  min  vib 
a  5  d  vib 

•  Fully  mix 


50 


60 


Figure  S6.1.  Comparison  of  uptake  rate  of  hexachloro-congeners  into  PE  in  static,  vibrating,  and  fully 
mixed  systems 
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Figure  S6.2.  Comparison  of  uptake  rate  of  heptachloro-congeners  into  PE  in  static,  vibrating,  and  fully 
mixed  systems 


106 


80 


60 


40 


20 


0 

6D 

c 


C  200 
u 
o 
C 

o 

o  150 


100 


I 


50 


\ 


PCB  194 


10 


20 


30 


40 


PCB203+196 


▲ 

i 


10  20  30  40 

Time  (d) 


i 

i 


50  60 


t 


50  60 


PCB208+195 


"S 

c 


c 

o 


c 

o 

O 


40 


30 


20 


10 


t 


ill 


10 


20 


30  40 

Time  (d) 


I 

t 


■  Static 

♦  Two  min  vib 
a  5  d  vib 

•  Fully  mix 


50 


60 


Figure  S6.3.  Comparison  of  uptake  rate  of  octachloro-congeners  into  PE  in  static,  vibrating,  and  fully 
mixed  systems 


H  Figure  S6.4.  Prototype  of  in  situ  vibrating  passive  sampler  built  with 

copper  casing 
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Fraction  remained  in  sediment 


1.2 


Figure  S6.5.  Normalized  desorption  data  for  three  PCBs  and  dual  compartment  model  fit 


Table  S6.8.  Parameters  estimated  by  fitting  dual  compartment  model  (Equation  S-3)  to  normalized 
desorption  data  for  three  PCBs 


Compound 

f 

kf(s_1) 

ks  (s'1) 

PCB  128 

0.48 

3.33xl0-5 

5.56xl0'6 

PCB  183 

0.8 

3. 06x1  O'5 

7xl0"7 

PCB  194 

0.48 

2.78xl0'5 

7xl0'7 

Reference 


108 


1-  Lohmann,  R.  Critical  review  of  low-density  polyethylene’s  partitioning  and  diffusion 
coefficients  for  trace  organic  contaminants  and  implications  for  its  use  as  a  passive  sampler. 
Environ.  Sci.  Technol.  2011,  46,  606-618. 

2-  Ghosh,  U.;  Talley,  J.  W.;  Luthy,  R.  G.  Particle-Scale  Investigation  of  PAH  Desorption 
Kinetics  and  Thermodynamics  from  Sediment.  Environ.  Sci.  Technol.  2001,  35,  3468-3475. 

3-  Gustafson,  K.  E.;  Dickhut,  R.  M..  Molecular  diffusivitiy  of  polycyclic  aromatic  hydrocarbons 
in  aqueous  solution.  J.  Chem.  Eng.  Data.  1994,  39,  281-285. 

4-  Lampert,  D.  An  assessment  of  the  design  of  in  situ  management  approaches  for  contaminated 
sediments.  Ph.D.  Thesis,  The  University  of  Texas  at  Austin,  May  2010. 

5-  Wemer,  D.;  Ghosh,  U.;  Luthy,  R.  G.  Modeling  polychlorinated  biphenyl  mass  transfer  after 
amendment  of  contaminated  sediment  with  activated  carbon.  Environ.  Sci.  Technol.  2006,  40, 
4211-4218. 

6-  Hawker,  D.  W.;  Connell,  D.  W.  Octanol-water  partition  coefficients  of  polychlorinated 
biphenyl  congeners.  Environ.  Sci.  Technol.  1988,  22,  382-387. 

7-  Sacan,  M.T.;  Ozkul,  M.;  Erdem,  S.  S.  Physico-chemical  properties  of  PCDD/Fs  and  phthalate 
esters.  SAR  QSAR  Environ.  Res.  2005,  16,  443-459. 

8-  Ghosh,  U.;  Kane  Driscoll,  S.;  Burgess,  R.  M.;  Jonker,  M.  T.  O.;  Reible,  D.;  Gobas,  F.;  Choi, 
Y.;  Apitz,  S.  E.;  Maruya,  K.  A.;  Gala,  W.  R.;  Mortimer,  M.;  Beegan,  C.  Passive  sampling 
methods  for  contaminated  sediments:  Practical  guidance  for  selection,  calibration,  and 
implementation.  Integr.  Environ.  Assess.  Manage.  2014,  10,  210-223. 


109 


10.3  Data  associated  with  figures  in  the  report. 


Data  for  Figure  5.2 

Compound  sediment  cone 


Naphthalene 

Acenaphthylene 

Acenaphthene 

Fluorene 

Phenanthrene 

Anthracene 

Fluoranthene 

Pyrene 

Benz(a)anthracene 

Chrysene 

Benzo(b)fluoranthene 

Benzo(k)fluoranthene 

Benzo(a)pyrene 

lndeno(1 ,2,3,-cd)pyrene 

Dibenz(a,h)anthracene 

Benzo(g,h,i)perylene 

stdev 

Porewater  cone  (ug/l) 

stdev 

0.21 

17.59 

1.24 

0.18 

8.45 

0.77 

0.30 

101.17 

1.33 

0.29 

27.55 

0.10 

0.16 

52.53 

0.80 

0.13 

16.02 

0.09 

0.93 

7.49 

0.25 

0.45 

7.87 

0.01 

0.51 

0.60 

0.06 

0.11 

0.51 

0.05 

0.03 

0.06 

0.00 

1.19 

0.07 

0.00 

0.84 

0.17 

0.01 

0.02 

0.09 

0.01 

0.28 

0.01 

0.00 

0.39 

0.03 

0.01 

(ug/g) 

1. 

2. 

3. 

2. 

19. 

7. 

16. 

16. 

7. 

7. 

4. 

5. 

14. 

13. 

2. 

5. 
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PE  fractional  uptake  of  PAHs  into  PE  (C/Ceq)  at  different  exposure  times  (Experimental  data) 
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Data  for  Fig  5.4  and  S5.4 
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Data  for  Fig  6. 1 


PCDD/Fs 

Log  Kow  PCBs 

Log  Kow 

2,3,7,8-TCDF 

6.41  PCB  128 

6.74 

1,2,3,7,8-PeCDF 

6.74  PCB  132 

6.58 

2,3,4,7,8-PeCDF 

6.8  PCB  146 

6.89 

2,3,7,8-TCDD 

7.05  PCB  153 

6.92 

1,2,3,7,8-PeCDD 

7.06  PCB  163 

6.99 

2,3,4,6,7,8-HxCDF 

7.43  PCB  158 

7.02 

1,2,3,7,8,9-FlxCDF  7 

7.44  PCB  185 

7.11 

1,2,3,4,7,8-HxCDF 

7.46  PCB  178 

7.14 

1,2,3,6,7,8-HxCDF 

7.56  PCB  187 

7.17 

1,2,3,4,6,7,8-HpCDF 

7.81  PCB  183 

7.2 

1,2,3,4,7,8,9-HpCDF 

7.92  PCB  201 

7.62 

1,2,3,4,7,8-HxCDD 

7.93  PCB  208 

7.71 

1,2,3,6,7,8-HxCDD 

7.93  PCB  194 

7.8 

Data  for  Fig  6.2 


PCB  homolog 

Cone  (ng/g) 

STDEV 

Mono 

0.00 

0 

Di 

0.00 

0.00 

Tri 

3.97 

0.34 

Tetra 

8.63 

0.33 

Penta 

7.84 

1.14 

Hexa 

26.84 

1.55 

Hepta 

26.07 

0.37 

Octa 

13.85 

0.14 

Nona 

0.60 

0.01 

Deca 

0.52 

0.06 
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Data  for  Fig  6.3;  6.4;  S  6. 1 ;  S  6.2.  S  6.3 


Cone  in  PE  in  Static  system  (ng/g)  STDEV 


PC  congener 

0  7  d 

14  d 

28  d 

56  d 

0  7  d 

14  d 

28  d 

56  d 

r(31 ) 

0.00 

19.01 

18.51 

26.05 

25.69 

0.00 

2.30 

2.36 

5.56 

1.00 

(70+76) 

0.00 

12.06 

11.77 

16.08 

17.74 

0.00 

1.89 

1.50 

2.43 

3.08 

"(99) 

0.00 

3.49 

3.71 

8.09 

7.05 

0.00 

0.30 

0.63 

1.74 

1.48 

"(146) 

0.00 

15.31 

17.28 

22.02 

26.49 

0.00 

0.00 

0.58 

0.15 

0.98 

"(153) 

0.00 

31.85 

35.96 

45.84 

53.95 

0.00 

0.00 

1.72 

0.49 

3.25 

"(132) 

0.00 

7.65 

9.03 

11.34 

13.65 

0.00 

0.00 

0.60 

0.18 

0.78 

(163+138) 

0.00 

38.71 

44.73 

56.74 

66.41 

0.00 

0.00 

2.72 

0.57 

3.75 

"(158) 

0.00 

3.37 

3.91 

5.07 

5.93 

0.00 

0.00 

0.19 

0.08 

0.39 

(178+129) 

0.00 

5.16 

8.64 

11.04 

12.66 

0.00 

0.00 

0.47 

0.39 

0.88 

(187+182) 

0.00 

25.75 

27.05 

33.16 

36.57 

0.00 

0.00 

1.59 

1.13 

2.62 

"(183) 

0.00 

8.26 

10.56 

11.28 

19.06 

0.00 

0.00 

0.73 

0.51 

1.25 

"(128) 

0.00 

2.56 

2.92 

3.81 

4.66 

0.00 

0.00 

0.14 

0.04 

0.31 

"(185) 

0.00 

3.90 

4.03 

4.73 

5.07 

0.00 

0.00 

0.21 

0.21 

0.35 

"(201) 

0.00 

25.50 

23.24 

41.20 

43.85 

0.00 

0.00 

1.86 

1.43 

2.61 

(203+196) 

0.00 

26.82 

26.07 

45.54 

48.56 

0.00 

0.00 

1.98 

1.79 

2.90 

(208+195) 

0.00 

5.68 

6.91 

8.11 

8.83 

0.00 

0.00 

0.48 

0.31 

0.55 

"(194) 

0.00 

12.18 

18.11 

20.96 

21.60 

0.00 

0.00 

0.89 

1.00 

0.97 

Cone  in  PE  in  2  min  pause  vibration  system 

(ng/g) 

STDEV 

PC  congene 

0  7d 

14  d 

28  d  56  d 

0  7d 

14  d 

28  d 

56  d 

"(31) 

0.00 

45.00 

48.00 

51.36 

54.65 

0.00 

1.53 

2.21 

6.73 

8.67 

(70+76) 

0.00 

25.43 

27.87 

38.54 

42.12 

0.00 

0.68 

4.18 

2.55 

8.59 

"(99) 

0.00 

16.06 

16.90 

22.09 

30.33 

0.00 

3.12 

0.92 

2.95 

4.68 

"(146) 

0.00 

43.17 

45.98 

61.47 

54.94 

0.00 

1.37 

1.45 

23.77 

7.39 

"(153) 

0.00 

132.31 

139.68 

208.62 

194.92 

0.00 

1.57 

20.47 

84.58 

29.63 

"(132) 

0.00 

28.05 

28.75 

36.57 

42.53 

0.00 

6.43 

2.26 

2.66 

4.38 

(163+138) 

0.00 

165.35 

170.18 

204.10 

244.19 

0.00 

17.18 

13.47 

3.37 

33.70 

"(158) 

0.00 

18.70 

19.10 

22.15 

30.31 

0.00 

1.91 

1.06 

2.26 

4.94 

(178+129) 

0.00 

25.74 

28.18 

48.76 

46.60 

0.00 

0.16 

4.28 

19.92 

6.22 

(187+182) 

0.00 

76.30 

81.75 

109.91 

126.97 

0.00 

5.17 

11.67 

15.28 

13.87 

"(183) 

0.00 

32.87 

35.33 

63.03 

67.67 

0.00 

1.74 

6.40 

15.91 

11.35 

"(128) 

0.00 

11.30 

12.13 

14.21 

16.30 

0.00 

1.41 

0.53 

1.76 

2.30 

"(185) 

0.00 

7.21 

7.50 

11.77 

13.41 

0.00 

0.13 

0.85 

1.82 

1.72 

"(201) 

0.00 

51.62 

70.10 

115.30 

147.35 

0.00 

5.21 

8.66 

17.03 

16.72 

(203+196) 

0.00 

58.58 

77.19 

129.46 

168.37 

0.00 

5.40 

9.16 

21.97 

24.28 

(208+195) 

0.00 

11.02 

14.63 

26.03 

33.33 

0.00 

0.96 

1.70 

5.35 

4.41 

"(194) 

0.00 

23.25 

30.03 

54.31 

71.63 

0.00 

1.13 

2.69 

9.16 

12.54 
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Cone  in  PE  in  5  d  pause  vibration  system  (ng/g) 

STDEV 

PC  congene 

0  7d 

14  d 

28  d 

56  d 

0  7d 

14  d 

28  d 

56  d 

r(31) 

0.00 

32.50 

35.53 

41.98 

37.73 

0.00 

2.42 

0.88 

0.93 

4.02 

(70+76) 

0.00 

23.67 

23.40 

27.70 

25.53 

0.00 

1.26 

0.17 

0.71 

2.27 

"(99) 

0.00 

7.46 

7.61 

10.14 

11.83 

0.00 

0.78 

0.03 

1.21 

1.04 

"(146) 

0.00 

21.56 

24.08 

32.76 

42.29 

0.00 

0.88 

0.05 

0.69 

1.99 

"(153) 

0.00 

45.87 

52.79 

73.19 

89.44 

0.00 

2.58 

0.75 

1.13 

6.66 

"(132) 

0.00 

10.36 

12.65 

16.05 

20.83 

0.00 

0.51 

0.50 

1.28 

1.25 

(163+138) 

0.00 

56.49 

64.16 

90.50 

108.32 

0.00 

3.63 

1.91 

1.56 

8.00 

"(158) 

0.00 

4.94 

5.20 

7.61 

9.77 

0.00 

0.45 

0.03 

0.13 

0.57 

(178+129) 

0.00 

8.82 

14.76 

20.48 

21.30 

0.00 

1.46 

1.49 

2.04 

1.14 

(187+182) 

0.00 

45.26 

47.04 

60.39 

65.86 

0.00 

2.70 

0.61 

2.01 

5.14 

"(183) 

0.00 

8.64 

18.87 

24.94 

37.09 

0.00 

0.92 

0.17 

0.16 

2.14 

"(128) 

0.00 

3.99 

4.90 

8.48 

9.60 

0.00 

0.92 

0.90 

2.04 

0.36 

"(185) 

0.00 

7.00 

7.50 

9.14 

10.34 

0.00 

0.00 

0.00 

0.10 

0.00 

"(201) 

0.00 

45.47 

45.97 

73.76 

82.07 

0.00 

1.80 

0.07 

0.18 

4.48 

(203+196) 

0.00 

46.34 

49.58 

82.82 

91.72 

0.00 

0.42 

0.06 

0.07 

5.58 

(208+195) 

0.00 

14.60 

13.84 

25.43 

15.19 

0.00 

0.00 

1.40 

0.00 

0.00 

r(194) 

0.00 

13.85 

21.70 

34.06 

41.34 

0.00 

0.00 

0.07 

0.60 

2.21 

Cone  in  PE  in  well-mixed  system  (ng/g)  STDEV 


PC  congenc 

0  2  d  7  d 

14  d 

28  d 

56  d 

0 

2d 

7  d 

14  d 

28  d 

56  d 

’(31) 

0.00 

50.28 

50.70 

59.80 

60.20 

58.39 

0.00 

12.36 

2.67 

1.75 

8.36 

10.97 

(70+76) 

0.00 

36.14 

35.80 

38.49 

56.73 

34.88 

0.00 

7.54 

1.52 

0.59 

10.10 

6.01 

’(99) 

0.00 

16.86 

16.12 

18.68 

23.91 

20.53 

0.00 

4.26 

2.98 

0.48 

2.79 

4.51 

’(146) 

0.00 

38.90 

55.03 

57.67 

54.02 

63.44 

0.00 

8.00 

0.96 

1.98 

7.52 

5.82 

’(153) 

0.00 

110.86 

181.26 

203.62 

190.92 

216.10 

0.00 

26.86 

3.98 

5.59 

25.90 

25.75 

’(132) 

0.00 

26.20 

32.74 

31.09 

38.42 

43.46 

0.00 

6.22 

1.62 

1.22 

0.77 

4.67 

(163+138) 

0.00 

142.79 

216.44 

224.85 

260.00 

252.94 

0.00 

32.60 

5.68 

5.65 

16.16 

29.60 

’(158) 

0.00 

13.26 

24.92 

25.32 

25.30 

29.35 

0.00 

3.48 

1.06 

1.56 

3.84 

3.79 

(178+129) 

0.00 

19.67 

34.64 

43.14 

46.98 

54.75 

0.00 

4.68 

2.51 

1.28 

7.07 

10.58 

(187+182) 

0.00 

65.23 

98.22 

140.61 

155.39 

151.87 

0.00 

14.47 

1.23 

4.09 

11.52 

28.71 

’(183) 

0.00 

29.90 

46.06 

67.60 

80.39 

81.67 

0.00 

7.11 

0.51 

1.07 

2.97 

13.81 

’(128) 

0.00 

10.08 

15.93 

14.78 

18.30 

16.19 

0.00 

2.32 

0.31 

0.67 

3.80 

1.03 

’(185) 

0.00 

7.68 

8.65 

12.44 

13.62 

13.06 

0.00 

1.94 

0.40 

0.19 

0.67 

3.43 

’(201) 

0.00 

44.00 

53.80 

92.78 

147.56 

143.77 

0.00 

2.33 

1.57 

2.61 

18.07 

31.99 

(203+196) 

0.00 

52.12 

58.18 

103.69 

162.00 

165.21 

0.00 

13.69 

1.87 

2.41 

51.24 

35.62 

(208+195) 

0.00 

9.89 

11.50 

19.77 

31.00 

32.99 

0.00 

2.58 

0.69 

0.48 

5.57 

6.23 

’(194) 

0.00 

22.05 

26.23 

38.42 

65.11 

69.79 

0.00 

5.79 

0.49 

0.90 

8.35 

15.22 

114 


Data  for  Fig  6.5 


Static  system 

2  min  Vibration  system 

5  d  Vibration  system 

Mix  system 

Initial  7  d 

14 d  28 d 

56  d  7  d 

14 d  28 d 

56  d  7  d 

14  d  28  d 

56  d  7  d 

14  d  28  d 

56  d 

3)  PRC-PCB29 

1.00  0.70 

0.60  0.45 

0.35  0.06 

0.04  0.00 

0.00  0.59 

0.45  0.29 

0.12  0.00 

0.00  0.00 

0.00 

4)  PRC-PCB69 

1.00  0.75 

0.67  0.54 

0.46  0.08 

0.05  0.00 

0.00  0.62 

0.41  0.32 

0.10  0.00 

0.00  0.00 

0.00 

7)  PRC-PCB155 

1.00  0.89 

0.83  0.78 

0.68  0.43 

0.25  0.13 

0.04  0.69 

0.60  0.42 

0.32  0.00 

0.00  0.00 

0.00 

9)  PRC-PCB192 

1.00  0.93 

0.90  0.88 

0.82  0.68 

0.54  0.37 

0.21  0.92 

0.78  0.74 

0.66  0.24 

0.24  0.12 

0.04 

Data  for  Fig  6.6  and  6.7 


Cfree  using  static  data  (pg/L) 

Cfree  using  2  i 

nin-pause  vibration  data 

Cfree  using  5  d-pause  vibration  data 

Measured  Cfree 

7  d  stdev 

28  d  stdev 

7  d  stdev 

28  d  stdev 

7  d  stdev 

28  d  stdev 

56  d  stdev 

31 

244.76 

20.34 

190.38 

55.22 

209.09 

3.79 

197.25 

18.47 

254.34 

8.04 

210.20 

0.90 

216.64 

40.70 

70+76 

70.75 

14.69 

41.93 

4.73 

30.31 

1.42 

37.64 

2.49 

74.04 

3.74 

43.83 

1.99 

33.72 

13.54 

99 

13.41 

5.46 

14.24 

3.26 

11.34 

2.65 

12.59 

1.68 

14.84 

1.68 

9.65 

1.55 

10.77 

2.36 

82+151 

75.54 

20.30 

48.50 

1.30 

43.75 

3.66 

55.20 

1.38 

66.17 

2.53 

41.74 

4.14 

46.48 

8.94 

118 

15.32 

5.83 

7.78 

1.04 

N.D  N.D 

7.30 

2.61 

15.69 

0.84 

7.95 

0.41 

6.25 

2.48 

128 

4.67 

1.71 

2.97 

0.18 

3.68 

0.99 

3.62 

0.77 

3.31 

1.13 

2.90 

1.01 

3.28 

0.21 

132 

22.33 

7.11 

14.81 

0.53 

13.22 

4.81 

13.89 

2.29 

14.92 

0.30 

9.98 

1.55 

12.82 

1.46 

153 

65.55 

31.02 

33.98 

4.59 

25.62 

0.35 

26.10 

0.20 

35.76 

0.78 

22.55 

3.22 

26.86 

3.20 

158 

5.26 

2.33 

3.11 

0.51 

N.D  N.D 

2.84 

0.50 

3.14 

0.23 

1.88 

0.32 

2.78 

0.36 

163+138 

71.37 

30.47 

41.83 

5.42 

33.64 

4.32 

33.76 

2.63 

43.87 

1.46 

28.03 

3.96 

32.30 

3.78 

177 

18.23 

9.18 

8.90 

1.65 

6.14 

0.51 

6.72 

0.35 

10.57 

0.37 

5.76 

1.19 

5.40 

0.93 

180 

74.42 

48.57 

27.98 

8.03 

10.94 

0.98 

13.00 

1.40 

35.53 

1.02 

18.52 

3.55 

9.69 

2.23 

183 

17.20 

8.35 

8.16 

1.91 

3.41  ~ 

0.55 

4.55 

0.32 

10.46 

0.20 

5.53 

1.08 

3.65 

0.80 

185 

5.57 

2.61 

2.58 

0.63 

0.88 

0.13 

1.13 

0.09 

4.24 

0.10 

1.92 

0.27 

0.90 

0.25 

182+187 

36.29 

20.68 

15.88 

3.68 

8.64 

0.76 

8.62 

0.18 

21.98 

0.02 

10.88 

1.37 

9.19 

1.74 

194 

7.89 

5.71 

3.83 

1.59 

0.86 

0.08 

1.26 

0.28 

3.96 

0.49 

2.27 

0.74 

0.79 

0.17 

201 

19.63 

13.07 

9.97 

3.45 

2.45 

0.60 

3.85 

0.64 

10.40 

0.86 

5.91 

1.83 

2.67 

0.59 

203+196 

19.20 

13.90 

10.51 

3.75 

2.67 

0.58 

3.44 

0.27 

10.13 

1.66 

6.27 

1.97 

2.83 

0.61 

208+195 

4.48 

2.87 

1.93 

0.67 

0.73 

0.19 

1.02 

0.34 

3.27 

0.21 

1.56 

0.00 

0.59 

0.11 

Data  for  Fig  S5.5 


Fraction  of  PRC  remained  in  PE 

Static  system 

2  min  Vibration  system 

5  d  Vibration  system  Mix  system 

Compound  Initial  7  d  14  d 

28  d 

56  d 

7  d  14  d  28  d  56  d 

7  d 

14  d 

28  d  56  d  7  d  14  d 

28  d 

56  d 

phenanthrene_d-10  1.00  0.19  0.16 

0.00 

0.00 

0.00  0.00  0.00  0.00 

0.00 

0.00 

0.00  0.00  0.00  0.00 

0.00 

0.00 

Pyrene_dl0  1.00  0.42  0.38 

0.21 

0.08 

0.01  0.00  0.00  0.00 

0.01 

0.00 

0.01  0.00  0.01  0.00 

0.01 

0.00 
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Data  for  Fig  S5.7 


Pause  time 

C/Ceq 

Energy  usage  (mWh) 

2min 

0.98 

5849 

5min 

0.98 

2339 

6hour 

0.97 

32 

1  day 

0.93 

8 

5  days 

0.72 

1 

Static 

0.48 

0 

Data  for  Fig.  S6.5 


C/Ceq  (Experimental  data) 

PCB  compound 

2d 

4  d 

8  d 

12  d 

48  d 

128 

1 

0.86 

0.77 

0.64 

0.52 

0.19 

183 

1 

0.85 

0.70 

0.56 

0.41 

0.22 

194 

1 

0.92 

0.83 

0.74 

0.64 

0.47 

C/Ceq  (Model  data) 

PCB  compound 

2d 

4  d 

8  d 

12  d 

48  d 

128 

1 

0.88 

0.78 

0.63 

0.52 

0.20 

183 

1 

0.84 

0.71 

0.53 

0.41 

0.18 

194 

1 

0.91 

0.84 

0.73 

0.65 

0.47 
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Fig.  7.4 

material 

log  Cs 

log  Cw 

log  Kd 

log  avg  Kd 

PTFE 

2.75 

0.34 

2.41 

PTFE 

3.20 

0.89 

2.32 

PTFE 

3.64 

1.38 

2.26 

PTFE 

3.76 

1.96 

1.80 

PTFE 

4.25 

2.44 

1.82 

2.19 

PDMS 

2.55 

0.08 

2.47 

PDMS 

3.01 

0.69 

2.32 

PDMS 

3.31 

1.33 

1.98 

PDMS 

3.70 

1.90 

1.80 

PDMS 

3.87 

2.42 

1.45 

2.14 

POM38 

2.69 

0.34 

2.35 

POM38 

3.03 

0.91 

2.12 

POM38 

3.48 

1.40 

2.08 

POM38 

3.73 

1.96 

1.77 

POM38 

3.76 

2.46 

1.29 

2.05 

CA 

2.97 

0.11 

2.86 

CA 

3.45 

0.72 

2.73 

CA 

3.77 

1.30 

2.47 

CA 

4.22 

1.86 

2.35 

CA 

4.40 

2.41 

1.99 

2.58 

CN 

3.09 

0.04 

3.05 

CN 

3.61 

0.55 

3.05 

CN 

4.09 

1.06 

3.03 

CN 

4.58 

1.60 

2.98 

CN 

4.63 

2.37 

2.26 

2.95 

PES 

2.87 

-0.26 

3.13 

PES 

2.86 

0.88 

1.99 

PES 

3.48 

1.31 

2.17 

PES 

4.32 

1.45 

2.87 

PES 

4.62 

2.21 

2.42 

2.72 

Parafilm 

2.21 

0.41 

1.80 

Parafilm 

3.06 

0.83 

2.22 

Parafilm 

3.43 

1.35 

2.08 

Parafilm 

3.81 

1.92 

1.90 

Parafilm 

4.50 

2.33 

2.17 

2.06 

PVDF+AC 

4.05 

0.47 

3.59 

PVDF+AC 

4.66 

1.10 

3.56 

PVDF+AC 

5.04 

1.54 

3.50 

PVDF+AC 

5.57 

2.29 

3.28 

PVDF+AC 

6.01 

3.05 

2.96 

3.43 

ag+AC 

3.91 

0.98 

2.94 

ag+AC 

4.18 

1.51 

2.67 

ag+AC 

4.63 

2.04 

2.59 

ag+AC 

5.02 

2.68 

2.34 

ag+AC 

5.17 

3.38 

1.78 

2.60 
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Fig.  7.5 

material 

log  Cs  log  Cw 

log  Kd 

log  avg  Kd 

PAC 

5.14  #NUM! 

#DIV/0! 

PAC 

5.64  0.14 

5.497795251 

PAC 

6.11  0.57 

5.543732198 

PAC 

6.69  1.11 

5.57996598 

PAC 

7.17  1.96 

5.206737092 

5.48 

ag 

3.02  #NUM! 

#DIV/0! 

ag 

3.37  0.47 

2.896577369 

ag 

4.03  1.37 

2.658172396 

ag 

4.13  1.80 

2.324221724 

ag 

4.40  2.37 

2.021792924 

2.59 

ag+PAC 

3.95  #NUM! 

#DIV/0! 

ag+PAC 

4.37  0.85 

3.516700727 

ag+PAC 

4.85  1.35 

3.508901728 

ag+PAC 

5.45  1.99 

3.457803605 

ag+PAC 

5.93  2.67 

3.267268584 

3.45 

PVDF  Dl 

2.95  #NUM! 

#DIV/0! 

PVDF  Dl 

3.19  0.64 

2.548915092 

PVDF  Dl 

3.70  1.18 

2.516823822 

PVDF  Dl 

3.79  1.89 

1.90674801 

PVDF  Dl 

4.58  2.46 

2.122157866 

2.35 

PVDF+GAC 

3.82  #NUM! 

#DIV/0! 

PVDF+GAC 

4.55  1.28 

3.272048906 

PVDF+GAC 

4.81  1.56 

3.247481195 

PVDF+GAC 

5.16  1.54 

3.610583669 

PVDF+GAC 

5.49  2.15 

3.343249605 

3.39 

PVDF+PAC  Dl 

3.94  #NUM! 

#DIV/0! 

PVDF+PAC  Dl 

4.53  1.18 

3.343787251 

PVDF+PAC  Dl 

4.92  1.80 

3.117051272 

PVDF+PAC  Dl 

5.55  2.24 

3.30875689 

PVDF+PAC  Dl 

6.04  2.88 

3.152196072 

3.24 

PVDF+PAC  meOH 

3.73  #NUM! 

#DIV/0! 

PVDF+PAC  meOH 

4.49  0.91 

3.58099527 

PVDF+PAC  meOH 

5.03  1.54 

3.483708501 

PVDF+PAC  meOH 

5.62  1.84 

3.783443739 

PVDF+PAC  meOH 

5.87  2.53 

3.346898687 

3.58 

CA+MA 

4.17  #NUM! 

#DIV/0! 

CA+MA 

4.24  1.86 

2.378467198 

CA+MA 

4.50  2.23 

2.270623596 

CA+MA 

5.12  2.80 

2.315642576 

CA+MA 

5.72  3.28 

2.442739087 

2.36 

PU+Cys 

#NUM!  1.93 
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PU+Cys 

3.75  1.92 

1.834592946 

PU+Cys 

4.24  2.42 

1.824019589 

PU+Cys 

5.04  2.78 

2.264592329 

PU+Cys 

4.55  3.46 

1.090452481 

1.92 

PET+Cys 

#NUM!  1.51 

#NUM! 

PET+Cys 

3.09  1.98 

1.106768435 

PET+Cys 

#NUM!  2.48 

#NUM! 

PET+Cys 

#NUM!  3.02 

#NUM! 

PET+Cys 

#NUM!  3.52 

#NUM! 

1.11 
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Fig.  7.6 

material 

log  Cs 

log  Cw 

log  Kd 

log  avg  Kd 
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1.80 

3.29 

chit+MAA  ag 

5.66 

2.30 

3.36 
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5.95 

2.86 

3.09 

chit+MAA  ag 

6.35 

3.35 

3.00 

3.21 

xylo+Cys  ag 

4.65 

1.07 

3.58 

xylo+Cys  ag 

5.19 

1.67 

3.52 

xylo+Cys  ag 

5.85 

2.18 

3.66 

xylo+Cys  ag 

5.96 

2.85 

3.11 

xylo+Cys  ag 

6.36 

3.35 

3.01 

3.45 

ag+Cys 

4.34 

1.37 

2.97 

ag+Cys 

4.75 

1.91 

2.84 

ag+Cys 

5.32 

2.40 

2.93 

ag+Cys 

5.68 

2.91 

2.77 

ag+Cys 

6.46 

3.36 

3.10 

2.94 

ag+SAMMS 

5.12 

0.00 

5.12 

ag+SAMMS 

5.50 

0.67 

4.83 

ag+SAMMS 

6.02 

1.17 

4.85 

ag+SAMMS 

6.48 

1.66 

4.82 

ag+SAMMS 

7.12 

1.92 

5.20 

ag+SAMMS 

7.12 

1.86 

5.26 

ag+SAMMS 

7.12 

1.89 

5.23 

5.01 

PVDF+Cys 

4.91 

-0.85 

5.76 

PVDF+Cys 

5.97 

1.44 

4.54 

PVDF+Cys 

5.83 

1.42 

4.41 

PVDF+Cys 

6.52 

2.23 

4.30 

PVDF+Cys 

6.98 

2.88 

4.10 

4.36 

PVDF+SAMMS 

4.96 

0.62 

4.34 

PVDF+SAMMS 

5.51 

1.57 

3.94 

PVDF+SAMMS 

5.96 

1.76 

4.20 

PVDF+SAMMS 

6.45 

2.43 

4.02 

PVDF+SAMMS 

6.88 

3.04 

3.84 

4.11 
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Fig.  7.8 

material  log  Cs 

log  Cw 

log  Kd 

log  avg  Kd 

DE  3.07 
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1.68 
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DE  4.79 
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2.53 
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DE-MPTMS  5.15 

-0.39 

5.54 

DE-MPTMS  5.48 

0.08 

5.40 

DE-MPTMS  6.10 

0.54 
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DE-MPTMS  6.50 

0.76 
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DE-MPTMS  7.00 

1.16 

5.84 

5.65 

ag+DE-MPTMS  4.49 

0.46 

4.03 

ag+DE-MPTMS  5.14 

1.19 

3.95 

ag+DE-MPTMS  5.36 

1.56 
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ag+DE-MPTMS  6.15 

2.43 
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2.98 
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PEGag  2.66 

-0.42 

3.07 

PEGag  3.54 

-0.53 
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PEGag  3.54 

1.09 

2.45 

PEGag  4.12 
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2.52 

PEGag  4.40 
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2.14 
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PEGag+GAC  4.96 
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3.91 

PEGag+GAC  5.15 

1.60 

3.55 

PEGag+GAC  5.64 

2.21 

3.43 

PEGag+GAC  6.02 

2.80 

3.23 

PEGag+GAC  6.49 

3.20 

3.29 

3.55 

PEGag+SAMMS  4.79 

0.84 

3.95 

PEGag+SAMMS  5.51 

1.51 

4.00 

PEGag+SAMMS  5.84 

1.94 

3.90 

PEGag+SAMMS  6.49 

2.61 

3.88 

PEGag+SAMMS  6.92 

3.03 

3.89 

PEGag+SAMMS  6.92 

3.03 

3.89 

PEGag+SAMMS  6.92 

3.03 

3.89 

3.93 
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Fig.  7.9 

material 
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Fig.  7.10 

material 
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2.98 

PVDF+GAC 
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3.62 
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0.9 

3.19 
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4.4 

1.4 

3.03 

PVDF+GAC 
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1.6 

2.83 
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